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J^ , Abstract. This is a comprehensive exposition of the classical moment problem using 

^— N ' methods from the theory of finite difference operators. Among the advantages of this 

f^ , approach is that the Nevanlinna functions appear as elements of a transfer matrix and 

\^ • convergence of Pade approximants appears as the strong resolvent convergence of finite 

^^ , matrix approximations to a Jacobi matrix. As a bonus of this, we obtain new results on 

^^ • the convergence of certain Pade approximants for series of Hamburger. 

C^ '. §1. Introduction 



a 
> 



The classical moment problem was central to the development of analysis in the pe- 
riod from 1894 (when Stieltjes wrote his famous memoir [38]) until the 1950's (when 



/\ I Krein completed his series on the subject [15, 16, 17]). The notion of measure (Stielt- 

j^ ■ jes integrals), Pade approximants, orthogonal polynomials, extensions of positive lin- 

ear functionals (Riesz-Markov theorem), boundary values of analytic functions, and the 
Herglotz-Nevanlinna-Riesz representation theorem all have their roots firmly in the study 
of the moment problem. 

This expository note attempts to present the main results of the theory with two 
main ideas in mind. It is known from early on (see below) that a moment problem is 
associated to a certain semi-infinite Jacobi matrix, A. The first idea is that the basic 
theorems can be viewed as results on the self-adjoint extensions of A. The second idea 
is that techniques from the theory of second-order difference and differential equations 
should be useful in the theory. 
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2 B. SIMON 

Of course, neither of these ideas is new. Both appear, for example, in Stone's treatment 
[39], in Dunford-Schwartz [5], and in Akhiezer's briUiant book on the subject [1]. But, 
normaUy, (and in [1] , in particular) , these ideas are presented as applications of the basic 
theory rather than as the central tenets. By placing them front and center, many parts 
of the theory become more transparent — in particular, our realization of the Nevanlinna 
matrix as a transfer matrix makes its properties (as an entire function of order 1 and 
minimal type) evident (see Theorem 4.8 below). 

Two basic moment problems will concern us here: 

Hamburger Moment Problem. Given a sequence 70, 71, 72, ... of reals, when is there 
a measure, dp, on (—00, 00) so that 



7^ = / x'^dpix) (1.1) 



and if such a p exists, is it unique? We let A^_ff (7) denote the set of solutions of (1.1). 

Stieltjes Moment Problem. Given a sequence 70, 71, 72, ... of reals, when is there a 
measure, dp on [0, 00) so that 

/•OO 

7„= / x^'dpix) (1.2) 

Jo 

and if such a p exists, is it unique? We let ^^5(7) denote the set of solutions of (1.2). 

We will not attempt a comprehensive historical attribution of classical results; see 
Akhiezer [1] and Shohat-Tamarkin [35] for that. 

We will always normalize the 7's so that 70 = 1. By replacing •^^ by 7n/7o, we can 
reduce general 7's to this normalized case. We will also demand that dp have infinite 
support, that is, that p not be a pure point measure supported on a finite set. This 
eliminates certain degenerate cases. 

One immediately defines two sesquilinear forms, H^ and Sn, on C^ for each A^ by 

Hn{0C,(3)= ^ an/3m7n+m (1-3) 

n=0,l,...,Af-l 
m=0,l,...,Af-l 

n=0,l,...,Af-l 
m=0,l,...,Af-l 

and corresponding matrices TIn and Sn so that H]^{a,l3) = {a^HNf^) and 5'iv («,/?) = 
(a, Sn(3) in the usual Euclidean inner product. Our inner products are linear in the 
second factor and anti-linear in the first. 
A standard piece of linear algebra says: 
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Lemma 1.1. An N x N Hermitean matrix A is strictly positive definite if and only if 
each suhmatrix At^'-^l = {aij)i<i^j<j has det(A[^"^l) > for J = 1, 2, . . . , A^. 

Proof. If A is strictly positive definite, so is each A^^''^], so their eigenvalues are all strictly 
positive and so their determinants are all strictly positive. 

For the converse, suppose that A^^'^"^! is positive definite. By the min-max principle, 
the (A^ — 1) eigenvalues of A^'^'^~^^ interlace the A^ eigenvalues of A^^'-^l = A. If all the 
eigenvalues of A^^'^"^! are positive, so are the A^ — 1 largest eigenvalues of A. If the A^*^ 
eigenvalue were non-positive, det(A) < 0. Thus, A^^'^~^^ strictly positive definite and 
det(A[^'^]) > imply A^^'^^ is strictly positive definite. An obvious induction completes 
the proof. D 

This immediately implies 

Proposition 1.2. {if7v}/^=i are strictly positive definite forms if and only ifdetiTi]^) > 
for A^ = 1, 2, . . . . Similarly, {Sn}'^^-^ are strictly positive definite forms if and only if 
det{SN) >0/or A^ = 1,2, .... 

Suppose that the 7„ obey (1.1). Then by an elementary calculation. 



N-l 2 



X 



E 

N-l 



N-l 2 



dp{x) = Hi^{a,a) (1-5) 



n=0 



dp{x) = SN{ct,a). (1-6) 



Taking into account that if / |P(x)p dp{x) = 0, p must be supported on the zeros of 
P, we have: 

Proposition 1.3. A necessary condition that (1.1) holds for some measure dp on 
(— oo, oo) with infinite support is that each sesquilinear form Hn is strictly positive defi- 
nite. A necessary condition that there be a dp supported on [0, cx)) is that each H^ and 
each Sn be strictly positive definite. 



Suppose now that each H]s[ is strictly positive definite. Let C[X] be the family of 

complex polynomials. Given P{X) = X]n=o '^n^^^ Q{^) = X]n=o f^n^^ (we suppose 
the upper limits in the sums are equal by using some zero a's or /3's if need be), define 

(P,g)=ifAr(a,/3). (1.7) 

This defines a positive definite inner product on C[X], and, in the usual way, we can 
complete C[X] to a Hilbert space H^'^' in which C[X] is dense. C[X] can be thought 
of as abstract polynomials or as infinite sequences (cto? «!,•••, ctAf-i, 0, . . . ) which are 

eventually via a ~ J2i=o ctj-^"'- 

We will start using some basic facts about symmetric and self-adjoint operators on 
a Hilbert space — specifically, the spectral theorem and the von Neumann theory of 
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self-adjoint extensions. For an exposition of these ideas, see Chapters VII, VIII, and X 
of Reed-Simon, volumes I and II [33, 34] and see our brief sketch at the start of Section 2. 
Define a densely defined operator A on Ti.^"^' with domain D{A) = C[X] by 

A[PiX)] = [XP{X)]. (1.8) 

In the sequence way of looking at things, A is the right shift, that is, A(ao, ai, . . . , oat, 0, . . . 
= (0, cto, cti, . . . , ctTv, 0, . . .). This means that 

{P,A[Q]) = SN{a,P) (1.9) 

and, in particular, 

(1,A-1)=7„. (1.10) 

Since Ti^ and Sn are real-symmetric matrices, A is a symmetric operator, that is, 
{A[P]jQ) = {P,A[Q]). Moreover, if we define a complex conjugation C on C[X] by 

^i^n=o ^riX'^) = Z]n=o ^^n-^", then CA = AC. It follows by a theorem of von Neu- 
mann (see Corollary 2.4 in Section 2) that A has self-adjoint extensions. 

If each Sn is positive definite, then (P, A[P]) > for all P, and it follows that A has 
a non- negative self-adjoint extension Ap, the Friedrichs extension (discussed further in 
Section 3). We thus see: 

Proposition 1.4. // all if at are positive definite, then A has self-adjoint extensions. If 
all Sn are positive definite, then A has non-negative self-adjoint extensions. 

Let A be a self-adjoint extension of A. By the spectral theorem, there is a spectral 
measure djl for A with vector [1] G T-C^'^\ that is, so that for any bounded function of A, 

{l,f{A)l) = Jf{x)dfi{x). (1.11) 

Since 1 G D{A^) C D{A^), (1.11) extends to polynomially bounded functions and we 
have by (1.10) that. 

In = / x^ djl{x). 



We see therefore that a self-adjoint extension of A yields a solution of the Hamburger 
moment problem. Moreover, a non-negative self-adjoint extension has supp((i/i) C [0, oo) 
and so yields a solution of the Stieltjes moment problem. Combining this with Proposi- 
tions 1.2, 1.3, and 1.4, we have the first major result in the theory of moments. 

Theorem 1 (Existence). A necessary and sufficient condition for there to exist a 
measure dp with infinite support obeying (1.1) is that det(7iAr) > for N = 1, 2, ... . 
A necessary and sufficient condition that also dp he supported on [0, oo) is that both 
det(H7v) > and det(5Ar) > Q for N = 1,2, ... . 

Historically, existence was a major theme because the now standard tools on existence 
of measures were invented in the context of moment problems. We have settled it quickly. 



THE CLASSICAL MOMENT PROBLEM 5 

and the bulk of this paper is devoted to uniqueness, especiaUy the study of cases of non- 
uniqueness. Non-uniqueness only occurs in somewhat pathological situations, but the 
theory is so elegant and beautiful that it has captivated analysts for a century. 

Henceforth, we will call a set of moments {7n}5^o with detiHN) > for all A^ a set 
of Hamburger moments. If both det(7i7v) > and det(iSAr) > for all N, we call them 
a set of Stieltjes moments. 

We will call a solution of the moment problem which comes from a self-adjoint ex- 
tension of A a von Neumann solution. The name is in honor of the use below of the 
von Neumann theory of self-adjoint extensions of densely- defined symmetric operators. 
This name, like our use of Friedrichs solution and Krein solution later, is not standard, 
but it is natural from the point of view of self-adjoint operators. As far as I know, neither 
von Neumann nor Friedrichs worked on the moment problem per se. While Krein did, 
his work on the moment problem was not in the context of the Krein extension we use 
to construct what we will call the Krein solution. What we call von Neumann solutions, 
Akhiezer calls A^-extremal and Shohat-Tamarkin call extremal. This last name is unfor- 
tunate since we will see there exist many solutions which are extreme points in the sense 
of convex set theory, but which are not von Neumann solutions (and so, not extremal in 
the Shohat-Tamarkin sense). 

Given the connection with self-adjoint extensions, the following result is reasonable 
(and true!): 

Theorem 2 (Uniqueness). A necessary and sufficient condition that the measure dp 
in (1.1) be unique is that the operator A of (1.8) is essentially self-adjoint {i.e., has 
a unique self-adjoint extension). A necessary and sufficient condition that there be a 
unique measure dp in (1.1) supported in [0, oo) is that A have a unique non-negative 
self-adjoint extension. 

This result is surprisingly subtle. First of all, it is not obvious (but true, as we will 
see) that distinct self-adjoint extensions have distinct spectral measures dfl, so there is 
something to be proven before multiple self-adjoint extensions imply multiple solutions 
of the moment problem. The other direction is even less clear cut, for not only is it not 
obvious, it is false that every solution of the moment problem is a von Neumann solution 
(Reed-Simon [34] has an incorrect proof of uniqueness that implicitly assumes every 
solution comes from a self-adjoint extension). As we will see, once there are multiple 
solutions, there are many, many more solutions than those that come from self-adjoint 
extensions in the von Neumann sense of looking for extensions in D(A). But, as we will 
see in Section 6, there is a sense in which solutions are associated to self-adjoint operators 
in a larger space. 

We also note we will see cases where the Stieltjes problem has a unique solution but 
the associated Hamburger problem does not. 

The Hamburger part of Theorem 2 will be proven in Section 2 (Theorems 2.10 and 
2.12); the Stieltjes part will be proven in Sections 2 and 3 (Theorems 2.12 and 3.2). If 
there is a unique solution to (1.1), the moment problem is called determinate; if there 
are multiple solutions, it is called indeterminate. It is ironic that the English language 
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literature uses these awkward terms, rather than determined and undetermined. Stieltjes 
was Dutch, but his fundamental paper was in French, and the names have stuck. Much 
of the interesting theory involves analyzing the indeterminate case, so we may as well 
give some examples that illuminate no n- uniqueness. 

Example 1.1. Let / be a non-zero C°° function on M supported on [0, 1]. Let g{x) = 
f{x), with / the Fourier transform of /. Then 

x'^gix) dx = v^(-z)" -4 (0) = 0. 

Let dpi{x) = (Reg)-^.{x) dx, the positive part of Reg, and let dp2 = {Reg)-{x) dx. By 
the above, 



x'^dpi{x) = / x^ dp2{x) 

D J — OO 

for all n. Since dpi and dp2 have disjoint essential supports, they are unequal and we 
have non-unique solutions of the moment problem. (We will see eventually that neither 
is a von Neumann solution.) The moments from pi,p2 may not be normalized. But we 
clearly have non-uniqueness after normalization. 

This non-uniqueness is associated to non-analyticity in a Fourier transform and sug- 
gests that if one can guarantee analyticity, one has uniqueness. Indeed, 

Proposition 1.5. Suppose that {'^n}'^=o is a set of Hamburger moments and that for 
some C,R> 0, 

hn\ < CEJ'nl (1.12a) 

Then the Hamburger moment problem is determinate. 
If {■yn}'^=o is a set of Stieltjes moments and 

l7n| < CR''{2n)\ (1.12b) 

then the Stieltjes moment problem is determinate. 

Proof. Let dp obey (1.1). Then x'^'^ G L^(R, rfp), and by the monotone convergence 
theorem, 



f X \ f°^ ^ f X \ 1 

cosh(^-Jrfp(x) = ^lim^y_^X:(^-J j^Mx 

<c lim y -] 



OO 



n=0 



Thus, e"^ G L^(R,dp{x)) for \a\ < jj^. It follows that Fp{z) = j e^^^dp{x) has an 
analytic continuation to {z \ \lm.z\ < -^}. If // is a second solution of (1.1), F^{z) is 
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also analytic there. But — P ^^^ (0) = 7n, so the Taylor series of Fp and F^ at z = 
agree, so Fp = Fp in the entire strip by analyticity. 

This implies that n = phy a, variety of means. For example, in the topology of bounded 
uniformly local convergence, (i.e., fn ^' f means sup||/n||oo < and fnix) -^ f{x) 
uniformly for x in any [—k, k]), linear combinations of the {e*^^ | y G M} are dense in all 
bounded continuous functions. Or alternatively, Gp{z) = J ^_J = i J_^ e~'^^F^(y) dy 
for Imz > 0, and fx can be recovered as a boundary value of Gp. 

The {2n)l Stieltjes result follows from the n! result and Proposition 1.6 below. D 

We will generalize Proposition 1.5 later (see Corollary 4.5). 
Proposition 1.6. Let {7n}5^o ^^ ^ ^^^ of Stieltjes moments. Let 

^2m = 7m, m = 0, 1, . . . 

r2m+i = 0, m = 0, 1,.... 

//{Fjj^^Q is a determinate Hamburger problem, then {7n}5^o ^-^ '^ determinate Stieljes 
problem. 

Proof. Let dp solve the Stieltjes problem. Let 

dp{x) = ^[X[o,oo){x) dp{x'^) + X{-oo,o]{x) dpix"^)]. 

Then the moments of p are F. Thus uniqueness for the F problem on (—00,00) implies 
uniqueness for the 7 problem on [0, 00). D 

Remark. We will see later (Theorem 2.13) that the converse of the last assertion in 
Proposition 1.6 is true. This is a more subtle result. 

Example 1.2. This example is due to Stieltjes [38]. It is interesting for us here because 
it is totally explicit and because it provides an example of non-uniqueness for the Stieltjes 
moment problem. Note first that 

) 

w'=w-'"^sin(27rlnw)rfw = 0. 

This follows by the change of variables v = — ^ + In-u, the periodicity of sin( • ), and 
the fact that sin is an odd function. Thus for any 6 G [—1, 1], 

-= / w'=w-i""[l + 6'sin(27rlnw)]rfw = e^(^+i)' 
VTT Jo 



(by the same change of variables) so 7^ = 6 4^^+^) is an indeterminate set of Stieltjes 
moments. Notice for 9 G (—1, 1), if dpeiu) is the measure with these moments, then 
sin(27rln(w))/l-t-6'sin(27rln(w)) is in L'^{dpg) and orthogonal to all polynomials, so dpg{u) 
is a measure with all moments finite but where the polynomials are not L^ dense. As 
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we shall see, this is typical of solutions of indeterminate moment problems which are not 
von Neumann solutions. 

Looking at the rapid growth of the moments in Example 1.2, you might hope that 
just as a condition of not too great growth (like Proposition 1.5) implies determinacy, 
there might be a condition of rapid growth that implies indeterminacy. But that is false! 
There are moments of essentially arbitrary rates of growth which lead to determinate 
problems (see the remark after Corollary 4.21 and Theorem 6.2). 

Example 1.3. There is a criterion of Krein [18] for indeterminacy that again shows the 
indeterminacy of Example 1.2, also of the example of the moments of exp(— |x|") dx (on 
(— oo, oo)) with ct < 1 and also an example of Hamburger [8], the moments of 



X[o,oo)(a^)exp(--^ ) dx. 

V In X + TT^ / 



Proposition 1.7 (Krein [18]). Suppose that dp{x) = F{x) dx where < F{x) < 1 and 
either 



(i) supp(F) = (— oo, oo) and 



^"(^(^^^ dx<oo (1.13a) 



oo i + iC 



or 

(ii) supp(F) = [0, oo) and 



HF{x)) dx ^ ^_ ^^_^3^^ 



/o (1 + a;) v^ 

Suppose also that for all n: 

/•oo 

\x\'^F{x)dx<oo. (1.13c) 

Then the moment problem {Hamburger in case (i), Stieltjes in case (ii)) with moments 

f x'^F{x) dx 
fF{x)dx 

is indeterminate. 

Remarks. 1. Hamburger's example is close to borderline for (1.13b) to hold. 

2. Since J^^ x^"" ex-p{-\x\") dx = 2a-^T{^^) ~ (^)! and (1.13a) holds for F{x) = 
exp(— |x|") if a < 1, we see that there are examples of Hamburger indeterminate moment 
problems with growth just slightly faster than the n! growth, which would imply by 
Proposition 1.5 that the problem is determinate. Similarly, since Jq a;"exp(— |x|") dx — 
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q,-1Y(^ii±1^ ~ (^)! and (1.13b) holds for F{x) = exp(— |a:|") if a < i, we see there are 
examples of Stieltjes indeterminate problems with growth just slightly faster than the 
(2n)! of Proposition 1.5. 

3. Since F{x) — ^e"'^' has moments 72n = (2n)! covered by Proposition 1.5, it is 
a determinate problem. The integral in (1.13a) is only barely divergent in this case. 
Similarly, F{x) = X[o,oo)(3^)e~^ is a Stieltjes determinate moment problem by Proposi- 
tion 1.5 and the integral in (1.13b) is barely divergent. This says that Krein's conditions 
are close to optimal. 

4. Krein actually proved a stronger result by essentially the same method. F need not 
be bounded (by a limiting argument from the bounded case) and the measure defining 
7tj, can have an arbitrary singular part. Moreover, Krein proves (1.13a) is necessary and 
sufficient for {e*"^}o<a<oo to not be dense in L^(M, F{x) dx). 

5. Krein's construction, as we shall see, involves finding a bounded analytic function 
G in C_(_ = {z E C \ Imz > 0} so that \G{x + iO)\ < F{x). From this point of view, 
the fact that F(x) cannot decay faster than e~'^' is connected to the Phragmen-Lindelof 
principle. 

6. The analog of Krein's result for the circle instead of the half plane is due to Szego. 

Proof. Suppose we can find G with ReG ^ with |G'(x)| < F{x) so that 

Then both r p^L ^^ and r pZ) ^x solve the 7„ moment problem, showing indeter- 
minacy. 

In case (i), define for z G C_(_ (see (1.19) and the definition of Herglotz function below): 

which is a convergent integral by (1.13a). Then, Im Q{z) > and by the theory of 
boundary values at such analytic functions [13], limejo Q{x + is) = Q{x + iO) exists for 
a.e. a; G M and Im Q{x + iO) — — lnF{x) . 

Let G{z) = ex.p{iQ{z)) for z G C_|_ and G{x) = ex.p{iQ{x + iO)). Then the properties 
of Q imply |G(-2)| < 1, lim£j^oG^(3^ + ^0) = G{x) for a.e. x in M, and |(j(a;)| = F{x). A 
standard contour integral argument then shows that for any e > and n > 0, 

x^'il - iexy-^Gix) dx = 0. 

Since (1.13c) holds and |G'(x)| < -F(x), we can take £ J, and obtain f x'^G{x) dx = 0. 
In case (ii), define 

Q{z) = / 5 Trdx. 

IT Jq X^ - Z^ 



10 B. SIMON 



Again, the integral converges by (1.13b) and 1biQ{z) > 0. Moreover, Q{—z) = —Q{z). 
We have for a.e. a; G M, lmQ{x + iO) = — InF(x^). Thus, if H{z) = exp{iQ{z)) for 
z e C+ and H{x) = exp{iQ{x + iO)) for x G M and G{x) = lniH{^/x) for x G [0, oo), 
then, as above. 



/CXD 
x^^'+^Hix) dx = 0, 
-CXD 



with Cnn > 






(1.14a) 


m = 0,1,2,... 


,n- 


- 1 


(1.14b) 
(1.14c) 



so since H{-x) = H{x), J^ x^^'+^Gix^) dx = or J^y''G{y)dy = 0. Since \H{x)\ = 
F(x2), we have that \G{x)\ < F{x). D 

Additional examples of indeterminate moment problems (specified by their orthogonal 
polynomials) can be found in [9, 10]. 

One point to note for now is that in all these examples, we have non-uniqueness with 
measures dp of the form dp{x) = G{x)dx (i.e., absolutely continuous). This will have 
significance after we discuss Theorem 5 below. 

To discuss the theory further, we must look more closely at the operator A given by 
(1.8). Consider the set {1,X, X^, ...} in H^^\ By the strict positivity of H^, these 
elements in 'H^'^' are linearly independent and they span Ti.^'^' by construction. Thus by 
a Gram-Schmidt procedure, we can obtain an orthogonal basis, {P„(X)}5^q, for H^'^\ 
By construction, 

Pn{X) = Cnn^^ + lower Order, 
(P P ) =0 

{Pn,Pn) = l. 

These are, of course, the well-known orthogonal polynomials for dp determined by the 
moments {'jn}^=o- Note that often the normalization condition (1.14c) is replaced by 
Cnn = Ij yielding a distinct set of "orthogonal" polynomials. There are explicit formulae 
for the Pn{X) in terms of determinants and the 7's. We discuss them in Appendix A. 

By construction, {Pj}^^q is an orthonormal basis for the polynomials of degree n. The 
realization of elements of H^"'^ as Xl^o ^^nPniX) with Xl^o l-^?^P < °^ gives a different 
realization of H^^' as a set of sequences A = (Aq, . . . ) with the usual ^^({0, 1,2,...}) 
inner product. C[X] corresponds to these A's with A^ = for n sufficiently large. But 
we need to bear in mind this change of realization as sequences. 

Note that span[l, . . . ,X"^] = span[Po, • • • , -Pn(-^)]- In particular, XPn{X) has an 
expansion in Pq, -Pi, • • • , Pn+i- But {XP^, Pj) = {Pn, XPj) = if j < n — 1 since then 
XPj is of degree at most n — 1, and (1.14b) holds. Thus for suitable sequences, {anjJ^O' 
{6^}- 0, and {c^}- (with P_i(X) = 0), 

XPn{X) = CnPn+l{X) + b^PniX) + a^-lPn-liX) (1-15) 

for n = 0, 1, 2, ... . Notice that by (1.14a), c^ > and 

Cn = \Pn+l,XPn) = {Pn,XPn+l) = ttn- 
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(1.15) thus becomes 

XP^iX) = a„P„+i(X) + b^P^iX) + an-iPn-i{X). 



(1.15') 



Since {Pn{X)} are an orthonornial basis for W'^\ this says that in this basis, A is given 
by a tridiagonal matrix, and D(A) is the set of sequences of finite support. 

Thus, given a set {7n}5^o o^ Hamburger moments, we can find 6o,6i,... real and 
ao, ai, . . . positive so that the moment problem is associated to self-adjoint extensions of 
the Jacobi matrix, 

/bo ao . . . \ 
tto bi tti 
A= ai 62 a2 ... . (1.16) 

a2 63 ••• 

V / 

There are explicit formulae for the 6^'s and a^'s in terms of the determinants of the 
■jnS. We will discuss them in Appendix A. 

Conversely, given a matrix A of the form (1.16), we can find a set of moments for 
which it is the associated Jacobi matrix. Indeed, if do is the vector (1, 0, . . . , 0, . . . ) in £2, 
then 

7^ = (5o,A"5o). 

Thus the theory of Hamburger moments is the theory of semi-infinite Jacobi matrices. 

So far we have considered the Pn{X) as abstract polynomials, one for each n. It is 
useful to turn this around. First, replace the abstract X by an explicit complex number. 
For each such z, consider the semi-infinite sequence tv{z) = {Po{z), Pi{z), P2{z), . . .). 
(1.15') now becomes (with P-i(z) interpreted as 0): 



anPn+l{z) + {bn " z)Pn{z) + ttn-lPn-liz) =0, ?! > 



(1.17) 



SO Pn{z) obeys a second-order difference equation. If, for example, 7r(a;o) G £^, then xo 
is formally an eigenvalue of the Jacobi matrix A. (Of course, tt is never a finite sequence 
because it obeys a second-order difference equation and tio{z) = 1. Thus, it ^ D{A). It 
may or may not happen that n G D{A), so the formal relation may or may not correspond 
to an actual eigenvalue. It will always be true, though, as we shall see, that if tt{xo) G P, 
then {A* - xo)7r(xo) = 0.) 

For convenience, set a_i = 1. For any z G C, the solutions of the equation 



a^Un+l + {bn - z)Un + tt^-lU 



n-1 



0, 



n> 



(1.18) 



are two-dimensional, determined by the initial data {u-i,uo) 
Pn{z) corresponds to taking 



u-i = 0, 



Uo = 1. 
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There is a second solution, Qn{z), taking initial conditions 

U-i = —1, Uq = 0. 

We will also define ^(z) = {Qo{z), Qi{z), . . .). It can be seen by induction that for n > 1, 
Qn{X) is a polynomial of degree n — 1. We will see later (Proposition 5.16) that the Q's 
are orthogonal polynomials for another moment problem. 

As we will see, this normalization is such that if some combination r]{z) = tn{z)+^{z) G 
£^, then (A* —z)r] = Sq and {Sq, rj) = t. Here dn is the Kronecker vector with 1 in position 
n and zeros elsewhere. We will have a lot more to say about Qn{z) in Section 4. 

There is a fundamental result relating the solution vectors tt, ^ to determinacy of the 
moment problem. In Section 4 (see Proposition 4.4 and Theorem 4.7), we will prove 

Theorem 3. Fix a set of moments and associated Jacobi matrix. Then the following 
are equivalent: 

(i) The Hamburger moment problem is indeterminate. 

(ii) For some zq E C with Iuizq ^ 0, 7r(zo) ^ ^2- 
(iii) For some zq E'C with Yvuzq ^ Q, ^{zq) G ^2- 
(iv) For some xq G M, both tx^xq) and ^(xq) lie in £2- 

(v) For some xq G M, both 7r(xo) and f^(a^o) ^^e in £2- 
(vi) For some xq G M, both ^{xq) and gf(a^o) H^ w £2- 
(vii) For all z E C, both 'jt{z) and ^(z) lie in £2- 

Remarks. 1. Appendix A has explicit determinantal formulae for X]n=o 1-^^(0)1^ ^^^ 
12n=o IQn(0)P providing "explicit" criteria for determinacy in terms of Umits of deter- 
minants. 

2. Theorem 3 can be thought of as a discrete analog of Weyl's limit point/limit circle 
theory for self-adjoint extensions of differential operators; see [34]. 

This implies that if all solutions of (1.18) lie in £^ for one z G C, then all solutions lie 
in £2 for all z E C. 

A high point of the theory, discussed in Section 4, is an explicit description of all 
solutions of the moment problem in the indeterminate case in terms of a certain 2x2 
matrix valued entire analytic function. It will suffice to only vaguely describe the full 
result in this introduction. 

Definition. A Herglotz function is a function $(-2) defined in C-|- = {z E C \ Imz > 0} 
and analytic there with 1bi^{z) > there. 

These are also sometimes called Nevanlinna functions. It is a fundamental result (see, 
e.g., [2]) that given such a $, there exists c > 0, d real, and a measure dn on M with 
/ i+x^ < 00 so that either c j^ or d^ j^ or both, and 



$(^) =cz + d + 



X 



X — z 1 + x^ 



dn{x). (1.19) 
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Theorem 4. The solutions of the Hamburger moment problem in the indeterminate case 
are naturally parametrized by Herglotz functions together with the functions ^{z) = t G 
M U {00} . This later set of solutions are the von Neumann solutions. 

If the Stieltjes problem is also indeterminate, there is to > so that the solutions of 
the Stieltjes problem are naturally parametrized by Herglotz functions, ^{z), which obey 
^{x+iO) G {to, 00) for X E (— oo,0) together with [to, oo)U {00}. The later set of solutions 
are the von Neumann solutions. 

We will prove Theorem 4 in Section 4 (see Theorems 4.14 and 4.18). From the ex- 
plicit Nevanlinna form of the solutions, one can prove (and we will in Section 4; see 
Theorems 4.11 and 4.17) 

Theorem 5. In the indeterminate case, the von Neumann solutions are all pure point 
measures. Moreover, for any t E M., there exists exactly one von Neumann solution, 
pL^^\ with l^'^^\{t}) > 0. Moreover, for any other solutions, p, of the moment problem, 
p({t})<^W({t}). 

Remark. The parametrization jj,^^' of this theorem is inconvenient since t hh^ ^(*) is many 
to one (in fact, infinity to one since each p^^' has infinite support). We will instead use 
a parametrization t ^-^ Ht given by 

dfxtix) 



X 



which we will see is a one-one map of IR U {00} to the von Neumann solutions. 

Examples 1.1 and 1.2 revisited. As noted, the explicit non-unique measures we 
listed were absolutely continuous measures. But there are, of necessity, many pure point 
measures that lead to the same moments. Indeed, the measures pt associated to the 
von Neumann solutions are each pure point. As we will see, there are many other pure 
point measures with the given moments. If z/ is a Cantor measure, then it can be proven 
that J Ht du{t) will be a singular continuous measure with the given moments. In the 
indeterminate case, the class of measures solving (1.1) is always extremely rich. 
Given a set of moments {7n}5^o ^^^ c eM., one can define a new set of moments 

7n(c) = X^('')c^7n-,. (1.20) 

For the Hamburger problem, the solutions of the {7n}^o ^^^ each {7^(c)}5^o problem 
are in one-one correspondence. If fi solves the {7n}^o Problem, then dp{x — c) = dpc{x) 
solves the {7n(c)}5^o problem, and vice- versa. But translation does not preserve the 
condition supp(/9) C [0, 00), so that as c decreases, the set of solutions of the Stieltjes 
problem shrinks. We will show that for any indeterminate Stieltjes moment problem, 
there is always a cq, so that for c > cq, {'^n{c)}'^=o is an indeterminate Stieltjes problem. 
For c < Co, {7n(c)}5^o ^^^ ^^ solution dp with supp((ip) C [0, 00), and {7n(co)} is a 
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determinate Stieltjes problem (but, of course, an indeterminate Hamburger problem). 
This means there are lots of examples of moments which are determinate Stieltjes but 
indeterminate Hamburger. The existence of cq is intimately connected to Theorem 5. 
Among all the von Neumann solutions, there is a distinguished one, dfip (the Friedrichs 
solution), with /o = inf(supp(/Ut)) maximal. One just takes cq = — /o- That /o is a pure 
point of dfiF is important in the analysis. 

Another topic in the moment problem concerns the theory of Fade approximants, or 
what is esssentially the same thing — continued fractions. Typical is Theorem 6 below, 
which we will prove in Section 5. 

Consider a sequence of numbers {Kn}'^=o, 



K^ = (-1)^ / x^'dpix) 
Jo 

for some p. We are interested in "summing" the formal series (called a series of Stieltjes) 

oo 

$^«n^", (1.21) 

n=0 

which is formally 

'^^. (1.22) 

1 + xz 

If the series (1.21) converges, then (1.22) is analytic in a circle, \z\ < R, which implies 
that p is supported in [0, -^]. Thus, if dp does not have compact support, then the series 
(1.21) will not converge for any z ^ 0. 

Definition. The [N,M] Fade approximant to the series (1.21) is the unique rational 
function of the form 

where A is a polynomial of degree A^ and S is a polynomial of degree M, and (as z ^ 0) 

N+M 
n=0 

Note that A/B has (A^ + 1) + (M + 1) - 1 = A^ + M + 1 free parameters and {k^}^+^ 
is A^ + M + 1 free numbers. There is an explicit solution for A, B in terms of determi- 
nants; see Baker and Graves-Morris [4]. We will say more about the definition of Fade 
approximants in Section 5. 
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Theorem 6. Let ^ k^^"" be a series of Stieltjes. Then for any z G C\(— oo,0), 
liuiN^oo f^^~^''^Kz) = f-{z) exists and liniN^oo f^^'^\z) = f+{z) exists. The con- 
vergence is uniform on compact subsets o/C\(— oo, 0]. Moreover, /_|_ = /_ if and only if 
the associated Stieltjes problem is determinate. If the problem is indeterminate and p is 
any solution of the associated Stieltjes problem, then 

for any x G [0, oo). In any event, for x > 0, f^^'^\x) is monotone decreasing in N and 
j^[N-i,N]^^-^ Z5 monotone increasing. 

In terms of the language of Section 3, /_ is associated to the Friedrichs solution and 
/-f to the Krein solution. For an interesting application of Theorem 6, see [22, 36]. In 
Section 5, we will also discuss series of Hamburger and prove the new result that the 
fi^'^\(^z) Fade approximants always converge in that case. 

As a by-product of the proof of Theorem 6, we will find 

Theorem 7. Let {7n}5^o ^^ ^ ^^^ ^/ Stieltjes moments. Then the Stieltjes problem is 
indeterminate if and only if 



n=0 

and 



5^|P,(0)|2<cx) (1.26) 

Qn{0) 



sup 

n 



Pn{0) 



< OO. (1.27) 



Theorem 8. Let {7n}5^o ^^ ^ ^^^ ^f Stieltjes moments. Then the Stieltjes problem is 
determinate while the Hamburger problem is indeterminate if and only if 

oo 
J2\Qnm^<00 (1.28) 

n=0 

and 

lim 

n— »oc 



QniO) 



PniO) 



OO. (1.29) 



We will see that (1.26)/(1.27) are equivalent to a criterion of Stieltjes. 

Theorem 6 is proven in Section 5 as parts of Theorem 5.2, Froposition 5.6, Froposi- 
tion 5.8, and Froposition 5.11. Theorem 7 is Theorem 5.21 and Theorem 8 is Froposi- 
tion 5.22. 

In Section 6, which synthesizes and extends results in Akhiezer [1], we will study when 
the closure of the polynomials has finite codimension, and among our results proven as 
part of Theorem 6.4 is 
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Theorem 9. Let p lie in A4^ {'^) and let Hq be the closure of the set of polynomials 
in Hp = LF' (R, dp) . Then Hq has finite dimension if and only if the Herglotz function, 
$, in the Nevanlinna parametrization of Theorem 4 is a real rational function {i.e., a 
ratio of real poynomials) . Equivalently, if and only if the measure p of (1.19) has finite 
support. 

In fact, we will see that dini(7io") is the degree of the rational function $. 
One consequence of this will be the following, proven in Appendix B (see Theorems B.l 
and B.4). 

Theorem 10. M.^ {'^) is a compact convex set {in the weak topology as a subset of the 
dual o/C(RU {oo})). Its extreme points are dense in the set. 

Remark. Do not confuse this density with the fact that the Krein-Millman theorem says 
that the convex combinations of the extreme points are dense. Here we have density 
without convex combinations, to be compared with the fact that in many other cases 
(e.g., probability measures on a compact Hausdorff space), the extreme points are closed. 

Here is a sketch of the contents of the rest of this paper. In Section 2, we review von 
Neumann's theory of self-adjoint extensions and use it to prove the uniqueness result 
(Theorem 2) in the Hamburger case. In Section 3, we review the Birman-Krein-Vishik 
theory of semi-bounded self-adjoint extensions and use it to prove the uniqueness result 
(Theorem 3) in the Stieltjes case. In Section 4, we turn to a detailed analysis of the 
polynomials P and Q in terms of transfer matrices for the difference equation associated 
to the Jacobi matrix and prove Theorems 3, 4, and 5. In Section 5, we discuss Pade 
approximations and prove Theorems 6, 7, and 8. In Section 6, we discuss solutions, p, 
where the closure of the polynomials has finite codimension in L^(]R, (ip). Appendix A 
is devoted to a variety of explicit formulae in terms of determinants of moments, and 
Appendix B to the structure of ^^^(7) as a compact convex set. Appendix C summarizes 
notation and some constructions. 



§2. The Hamburger Moment Problem as a Self-Adjointness Problem 

Let us begin this section with a brief review of the von Neumann theory of self- 
adjoint extensions. For further details, see [2, 5, 33, 34]. We start out with a densely 
defined operator A on a Hilbert space 7Y, that is, D{A) is a dense subset of TC and 
A : D{A) ^ Ti a, linear map. We will often consider its graph T{A) d Ti xTi given by 
T{A) = {((/?, Aip) I (f G D{A)}. Given operators A, B, we write A G B and say B is an 
extension of A if and only if T{A) C r(i?). 

One defines a new operator A* as follows: r] G D{A*) if and only if there is a i/' G 7i 
so that for all (/? G D{A), ('(/',</') = {v^^v)- We set A*r] = i/j. In essence. A*, called the 
adjoint of A, is the maximal object obeying 

(A*77,y:) = (r7,Avp) (2.1) 

for all 77 G D{A*), ip G D{A). 
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An operator is called closed if and only if T{A) is a closed subset of H x H and it is 
called closable if and only if r{A) is the graph of an operator, in which case we define A, 
the closure of A, by T{A) = r{A). Thus, D{A) = {i^ e H \ 3(/?^ G D{A), so that Lpn -^ 
ip and Aifn is Cauchy} and one sets Aip = lim Aip^. 

Adjoints are easily seen to be related to these notions: A* is always closed; indeed, 

T{A*) = {{Aip, -ip)enxn\ipe D{A)}^. (2.2) 

A* is densely defined if and only if A is closable, in which case (by (2.2)), A = (A*)*. 

An operator is called symmetric if A C A* ( equivalent ly, if ((/?, Ai/j) = {A^p, i/j) for 
all (/?,'(/'€ D{A)), self-adjoint if A = A*, and essentially self-adjoint if A is self-adjoint 
(equivalently, if A* is self-adjoint). Notice that symmetric operators are always closable 
since D{A*) D D{A) is dense. Note also that if A is symmetric, then A* is symmetric if 
and only if A is essentially self-adjoint. 

Von Neumann's theory solves the following fundamental question: Given a symmetric 
operator A, when does it have self-adjoint extensions, are they unique, and how can they 
be described? If S is a self-adjoint extension of A, then B is closed, so A G B. Thus, 
looking at self-adjoint extensions of A is the same as looking at self-adjoint extensions of 
A. So for now, we will suppose A is a closed symmetric operator. 

Define /C± = ker{A* =f i), that is, }C+ = {^p e H \ A*ip = ip}. They are called the 
deficiency subspaces. Using \\{A ± i)v?|p = H^y'lP + ll'/'P, it is easy to see that if A is 
closed, then Ran(A ± i) are closed, and since ker(A* =1= z) = Ran(A ± z)^, we have 

Ran(A±z) = /Ci. (2.3) 

Let d± = dim(/C-|-), the deficiency indices of A. 

Place the graph norm on D{A*), that is, ||</7||^. = {p,p) + {A*(pj A*ip). This norm 
comes from an inner product, (99, i{j)a* = {'^^ i^) + {A*'^^ A*!/;). 

Proposition 2.1. Let A be a closed symmetric operator. Then 

L>(A*) = D(A)©/C+©/C_, (2.4) 

where © means orthogonal direct sum in the ( ■ , ■)a* inner product. 

Proof. If ip E /C+ and ifj G /C_, then {(p,'i/j)a* = {v,'^) + {i^.^ii^) = so /C_|_±^./C_. 
If p e D{A) and i) E /C±, then {A*Lp,A*xp) = (Av?, ±i^) = {ip,±iA*^Jj) = -(</?, ^), so 
L>(A)±A'/C+©/C_. 

Let?7 G D{A*). By (2.3), Ra.n{A + i) + IC+ = 7^ so we can find v? G D{A) and V e /C+, 
so [A* + i)r] = {A + i)ip + 2i'il). But then {A* + i)[r] - p - if)] = Q, so r] - p - if) e K,-, 
that is, ?7G L'(A) + /C+ + /C_. D 
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Corollary 2.2. (i) Let A be a closed symmetric operator. Then A is self-adjoint if and 
only if d+ = d- = 0. 

(ii) Let A be a symmetric operator. Then A is essentially self-adjoint if and only if 
d+ = d- = 0. 

If A C B, then B* C A*, so if B is symmetric, then Ac B C B* C A*. Thus, to look 
for symmetric extensions of A, we need only look for operators B with A G B G A* , 
that is, for restriction of A* to D{B)^s with D{A) G D{B). By Proposition 2.1, every 
such D{B) has the form D{A) + S with S G IC+ + IC-. On D{A*) x D{A*), define the 
sesquilinear form (sesquilinear means linear in the second factor and anti-linear in the 
first) 

Proposition 2.3. Let A be a closed symmetric operator. Then 

(i) The operators B with A G B G A* are in one-one correspondence with subspaces 

SofIC+ + /C_ under D{B) = D{A) + S. 
(ii) B is symmetric if and only if Q \ D{B) x D{B) = 0. 
(iii) B is symmetric if and only if Q \ S x S = 0. 

(iv) B is closed if and only if S is closed in /C+ © /C_ in D{A*) norm. 
(v) v? e D{A*) lies in D{B*) if and only ifQ{^,ij) = for all ip e D{B). 
(vi) Let J : /C+ © /C_ ^ /C+ © /C_ by J(v?, ^) = (v?, -^) ■ If D{B) = D{A) + S, then 

D{B*) = D{A) + J[S]^, where ^ is in JC+ © /C_ in the D{A*) norm. 
(vii) /C+(S) =}C+r]S^, /C_(S) = }C-nS^ {with ^ in the ( , )a* inner product). 

Proof. We have already seen (i) holds and (ii) is obvious, (iii) holds since if (/? G D{A) 
and t/j G D{A*), then Q{(f,t(j) = by definition of A* and symmetry of A. Thus, if 
(pi,(p2 G -D(A) and 'i/'i, '02 e 5", then Q{(pi + ipi,(p2 + ^2) = Qi^Pi, '^2)- (iv) is immediate 
if one notes that T{B) = T{A) © {((/?, A*(p) \ ip e S} with Q) inHxH norm and that 
the Ti X Ti norm on {{ip,A*ip) | </? G 5"} is just the D[A*) norm, (v) follows from the 
definition of adjoint. 

To prove (vi), let ry = (?7i,?72), H^ = {'^ij'^2) G A^+ © A^-- Then, direct calculations 
show that 

Qiv, V) = 2z[(??i, v?i)7^ - {V2, V2)h] (2.5) 

(??, ^)a' = 2[(?7i, (pi)n + (?72, ^2)h]- (2.6) 

Thus, Q(?7, (/?) = if and only if rj-LJip. (v) thus implies (vi). 

To prove (vii), note that by (vi), /C±(S) = /C± n D{B*) = /C± n J[S]^ = }C± n S^ 
since </? G /C± lies in J [5*]^ if and only if Jip lies in S-^ if and only if ip lies in S^ by 

J(/7 = ±(p if (p G /C±. D 

With these preliminaries out of the way, we can prove von Neumann's classification 
theorem: 
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Theorem 2.4. Let A he a closed symmetric operator. The closed symmetric extensions 
B of A are in one-one correspondence with partial isometrics of IC+ into /C_, that is, 
maps U : /C+ -^ /C_ for which there are closed subspaces U+ C /C+ and U- C /C- so 
that U is a unitary map of V(+ -^ U- and U = on U^ . If B corresponds to U , then 
D{B) = D{A) + {^ + U^ I ^ e U+}. Moreover, IC±{B) = IC± n U^ {H orthogonal 
complement). In particular, A has self-adjoint extensions if and only if d^ = d- and 
then the self-adjoint extensions are in one-one correspondence with unitary maps U from 
/C_(_ to /C_ . 

Proof. Let S be a symmetric extension of A. Let </? G D{B) fl (/C_(_ © /C-), say </? = {^pi + 
(fi2)- By (2.5), Q{ip,(fi) = 2z(||(/?i|p — ||(/?2|P), so B symmetric implies that ||(/7i|| = \\'^2\\j 
so D{B)n{IC^Q)}C-) = {some ((/?i, (/?2)}5 where ||(/?i|| = ||(/'2||- Since il>(i?) is a subspace, 
one cannot have ((/7i,(/?2) and {Lpi,ip2) with (/?2 7^ '^2 (for then (0, (/?2 — '^2) ^ D{B) so 
\\'^2 ~ '/'2II = 0)- Thus, D(B) n (/C_|_ © /C_) is the graph of a partial isometry. On the 
other hand, if t/ is a partial isometry, Q{lp + U(p, V + ^'0) = 2«(((/', V') ^ (^V'? ^V')) = 
so each U does yield a symmetric extension. 

That IC±{B) = }C± nU^ follows from (vii) of the last proposition if we note that if 
ip G /C+, then </? ±a* {('0 + Uifj) \ ifj G U+} if and only if </? ±a* {1/^ \ i^ & 1^+} if and only 
ii^±nl^+ (by (2.6)). 

Thus -B is self-adjoint if and only if t/ is a unitary from /C_|_ to /C_, which completes 
the proof. D 

Recall that a map T is anti-linear if T {a(p+b'i/j) = dT{Lp)+bT{'i/j) for a, 6 G C, (/?, V' G '^; 
that T is anti-unitary if it is anti-linear, a bijection, and norm-preserving; and that a 
complex conjugation is an anti- unitary map whose sequence is 1. 

Corollary 2.5. Let A be a symmetric operator. Suppose there exists a complex conju- 
gation C -.n^n so that C : D{A) -^ D{A) and CAip = Adp for all (p G D{A). Then 
A has self-adjoint extensions. If d+ — 1, every self-adjoint extension B is real, that is, 
obeys C : D{B) -^ D{B) and CBip = Bdp for all (p G D{B). If d+ > 2, there are 
non-real self-adjoint extensions. 

Proof. C is an anti-unitary from /C_|_ to /C_ so d+ = d-. 

We claim that if B is self-adjoint and is associated to U : /C-|- -^ /C_, then B is real if 
and only if CUCU = 1 for it is easy to see that C : D{A*) -^ D{A*) and CA*ip = A*Cip, 
so B is real if and only if C maps D{B) to itself. Thus, C{(p -\- U(p) must be of the form 
ifj + Uifj. Since ifj must be CUip, this happens if and only if Cip — UCUip, that is, 
CUCU if = ip. This proves our claim. 

If rf+ = 1 and ip G /C+, CUip = e^V for some 6. Then {CUCU)ip = Ce'^Uip = 
e~^^CUip = ip showing that every self-adjoint extension is real. 

If (i_i_ > 2, pick (pji/j E /C+ with </? ± -0 and let U^p = Cifj, Uif^ — idp. Then 
CUCU<p = CUCCip = CUi/j = CiCi/j = -i^, so CUCU ^ 1. Thus, the B associated to 
this U will not be real. D 

Next we analyze a special situation that will always apply to indeterminate Hamburger 
moment problems. 



20 B. SIMON 

Theorem 2.6. Suppose that A is a closed symmetric operator so that there exists a 
complex conjugation under which A is real. Suppose that (i_|_ = 1 and that ker(A) = {0}, 
dimker(A*) = 1. Pick v? G ker(A*), C^ = ip, and r] E D{A*), not m D{A) + ker(A*). 
Then {(fi,A*r]) ^ and i/j = {r] — [{rj, A*rj)/{(p, A*rj)]Lp}/{(p, A*r]) are such that in </?,•!/' 
basis, {■ , A* ■) has the form 

<■-*•■>= (o I)- <'■'' 

The self-adjoint extensions, Bf, can be labelled by a real number or oo where 

D{Bt) =D{A) + {a{ty^ + tlj)\aeC} teR 
= D{A) + {acp \aeC} t = oo. 

The operators Bt are independent of which real V' in D{A*)\D{A) is chosen so that (2.7) 
holds. 

Proof. If ((/?, A*ri) = 0, then the {■ , A* ■) matrix with basis </?, ry would have the form 
( " ) with a G M, in which case Q = on the span of </? and rj, which is incompatible with 

the fact that d.^. = 1. Thus, {ip,A*r]) ^ and (2.7) holds by an elementary calculation. 

Once (2.7) holds, it is easy to see that the subspaces 5" of D{A*) with D{A) C 5", 
dira{S/D{A)) = 1, and Q t -5 x 5 = are precisely the D{Bt). 

If -(/; is a second if) for which (2.7) holds, then Q{il) — t/j^p) = for all p G D(A*), so 
•i/; — "0 G D{A) and the D(i?t)'s are the same for -0 and t/j. D 

Remarks. 1. It is easy to see that if ker(A) = {0}, then dimker(A*) < d+, so if 
(i-i- = 1 and ker(A) = {0}, dimker(A*) is either or 1. The example of A = —-^j-^ 
on L^(0, oo) with D{A) = Cq°{0,oo) shows that it can happen that A is real, d+ = 1, 
but dimker(A*) = 0. 

2. One example where this formalism is natural is if A > o: > (for a G M), in which 
case dimker(A*) = d^. One takes tj = Ap </?, where Ap is the Friedrichs extension. 
We will discuss this further in Section 3. A second example, as we will see, is the 
indeterminate moment problem, in which case one can take ip = n (0) , i/j = ^{0) . 

We can now turn to the analysis of the symmetric operator A on Ti*^^-* described in 
Section 1. In terms of the explicit basis Pn{X) in W'^', we know A has the Jacobi matrix 
form (1.16). We will first explicitly describe A* and the form Q. Given any sequence 
s = (so. Si, . . . ), we define a new sequence J-'{s) (think of JF for "formal adjoint") by 

^, , f boso + aosi if n = 

J^[S)n = S 

[ an-lSn-l+bnSn + anSn+1 if n > 1. 

Given any two sequences, Sn and tn, we define their Wronskian, 

W{S, t){n) = an{Sn+ltn - Sntn+l)- (2.8) 
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The Wronskian has the foUowing standard property if 

anSn+l + CnSn + ttn-lSn-l = (2.9) 

and 

ttnin+l + dntn + ttn-ltn-l = 0- (2.10) 

Then multiplying (2.9) by t„ and subtracting (2.10) multiplied by s^, we see that 

W{S, t){n) - W{S, t){n - 1) = (rf, - Cr,)Sntn. (2.11) 

In particular, if d = c, then W is constant. 

Theorem 2.7. Let A be the Jacobi matrix (1.16) with D(A) = {s \ Sn = for n 
sufficiently large}. Then 

D{A*) = {se^2\ :F{s) e £2} (2.12) 

with 

A*s = J^{s) (2.13) 

for s e D{A*). Moreover, if s,t e D{A*), then 

lim W{s,t){n) = {A*s,t) - {s,A*t). (2.14) 

n — >oo 

Proof. Since the matrix (1.16) is a real-symmetric matrix, it is easy to see that for any 
t G D{A), and any sequences 

(s,At) = (^(s),t). (2.15) 

Since t and At are finite sequences, the sum in ( ■ , ■ ) makes sense even if s or J-'{s) are 
not in £2- Since D{A) is dense in £2, (2.15) says that s E £2 lies in D{A*) precisely if 
J-^{s) G £2 and then (2.13) holds. That proves the first part of the theorem. 
For the second, we perform a calculation identical to that leading to (2.11): 

N N 

J] [T^t^-s;;{J^t^)] = W{s,t){0) + J2 W{s,t){n) - W{s,t){n - 1) 

n=0 n=l 

= W{s,t){N). (2.16) 

If s, t G D{A*), then the left side of (2.16) converges to {A*s, t) - (s, AH) as A^ ^ 00, so 
(2.14) holds. D 

Lemma 2.8. If ^ & D{A*), [A* — z)(p = for some z and (^q = 0, then ^p = 0. 

Proof. Essentially, this is so because {A* — z)s = is a, second-order difference equation 
with s-i = 0, so solutions are determined by sq- Explicitly, suppose A*s — zs and 
So = 0. Then 

Sn+l = a~^[{-bn + z)Sn - ttn-lSn-l], n > 1 

= aQ^[{-bo + z)so], n = 0. 

By a simple induction, s = 0. D 
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Corollary 2.9. The operator A associated to a Hamburger moment problem always has 
either deficiency indices (1, 1) or deficiency indices (0,0). 

Proof. We have already seen that dj^ = d- so it suffices to show that ker{A* — i) has 
dimension at most 1. By Lemma 2.8, the map from ker(A* — i) to C that takes s ^-^ sq 
is one-one, so ker{A* — i) is of dimension or 1. D 

We are now ready to prove the more subtle half of the Hamburger part of Theorem 2. 

Theorem 2.10 (one quarter of Theorem 2). Let A be essentially self-adjoint. Then 
the Hamburger moment problem has a unique solution. 

Proof. Pick z with Imz > and a measure p obeying (1.1). Since A is essentially self- 
adjoint, {A — i)[D{A)] is dense in £^, and by an identical argument, {A — z)[D{A)] is 
dense in i'^. Thus, there exists a sequence Rn (X) of polynomials so that 

||(X-^)i?W(X)-l||^0 

in H^'^\ and thus 

''[{x-z)Ri^\x)-lfdp{x) ^0 



since dp{x) realizes H^~'^ on polynomials. Now ^^^ is bounded for a; G M since Im^ > 0. 
Thus 

2 

dp{x) —^ 0. 



I Ri:^^)-^ 



It follows that 



GJx)= /"M^= lini I R\^\x)dpix) 

J X — Z n^oo J 



is independent of p. Since Gp{x) determines p (because p{{a}) = liuigio elm Gp{a + is) 
and p((a, h)) + p([a, h]) = - limg^o / IiiiGp(?/ + is) dy), all p's must be the same. D 

Theorem 2.11. Let A be a Jacobi matrix with D{A) the sequences of finite support. 
Suppose A is not essentially self-adjoint and B, F are distinct self-adjoint extensions of 
A. Then 

'so, -^ So] ^ (so. -J— 6o] . (2.17) 



B — I J \ L — I 

Remarks. 1. The proof shows for any z G C\]R, (5o, [B — z)~^5q) ^ {So{F — z)~^do). It 
also works for ^ G M as long as z is in the resolvent set for both B and F. 

2. In Section 4, we will have a lot more to say about the possible values of (5o, -^zt^Sq) 
as B runs through the self-adjoint extensions of A. 

Proof. We first claim that do ^ Ran(A — i). For suppose on the contrary that there exists 
r] G D{A) and So ^ {A - i)r] and that {A* + i)ip = 0. Then {So,(p) = {{A - i)?7, v) = 
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(?7, {A*+i)^) = 0. By Lemma 2.8, (/? = 0. Thus, if 5o G Ran(A-z), then ker(A*+z) = {0}, 
so A is essentiaUy self-adjoint. By hypothesis, this is false, so 6o ^ Ran(A — i). 

Thus, {B - i)-^6o and (F - i^^So are in D(A*)\D(A). Since dim{D{B)/D{A)) = 
dini{D(F)/D{A)) = 1, if these vectors were equal, D{B) would equal -D(-F), so B = A* \ 
D{B) would equal F = A* \ D{F). Thus, {B - i)-^5o ^ (F - i)-^5o. 

Let r]={B- i)-^5o - (F - i)-^5o. Then 

(A* - i)rj = {A* -i){B- i)-^do - {A* -i){F- i)-^do = 5o - 5o = 0, 

so ?7 7^ implies ryo 7^ by Lemma 2.8. Thus, {Sq,!]) = (5o, {B — i)~^do) — (5o, {F — 
i)-Mo)^0. D 

As a corollary of Theorem 2.11, we have 

Theorem 2.12 (two quarters of Theorem 2). A Hamburger moment problem for 
which A is not essentially self-adjoint is indeterminate. A Stieltjes moment problem for 
which A has multiple non-negative self-adjoint extensions is indeterminate. 

Proof. Theorem 2.11 implies that distinct self-adjoint extensions lead to distinct spectral 

measures since 

dn^{x) 



{6oAB-i)-'6o) = j 



X 



where ^^ is the solution to (1.1) associated to B. Positive self-adjoint extensions yield 
solutions of the Stieltjes moment problem. D 

With this machinery available, we can prove the converse of Proposition 1.6: 

Theorem 2.13. Let {7^)5^0 ^^ ^ ^^^ ^f i^oments for which the Stieltjes problem has 
solutions. Let 

r2m = 7m, m = 0, 1,... (2.18a) 

r2^+i = 0, m = 0,l,.... (2.18b) 

Then {FjI^^q is a determinate Hamburger problem if and only ^/{7n}5^o ^-^ ^ determi- 
nate Stieltjes problem. 

Proof. The proof of Proposition 1.6 shows there is a one-one correspondence between 
solutions of the Stieltjes problem for 7 and those solutions dp of the Hamburger problem 
for F with dp{—x) = dp{x). We will call such solutions symmetric. Thus it suffices 
to show that if the Hamburger problem for F is indeterminate, then there are multiple 
symmetric solutions. 

Let U map C{X) to itself by U[P{X)] = P{-X). By T2m+i = 0, t/ is unitary, and 
so extends to a unitary map from W^' to itself. Clearly, U maps D{A) = C{X) to itself 
and UAU-^ = -A. 

Thus, U maps D{A*) to itself and UA*U~^ = —A*. This means that U maps /C+ to 
/C_. 
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Let C be the complex conjugation under which A is real. Then UC = CU is an 
anti-unitary map of /C-(- to itself. So if (/? is a unit vector in /C_(_, then Udp = e ip for 

some real 6. Thus, UC{e'^^/'^Lp) = e~'^^/'^UCip = e*^/^(/? so there exists •i/' 7^ in /C_|_ with 

ucij = V'. 

In particular, since UC = CU, U{i/j ± Cif;) — ±(-0 ± Cifj) and therefore, if B = A* \ 
D{B) and F = A* \ D{F), where D{B) ^ D (A) + [ij + Cij] and D{F) = D{A) + [i;-Ci;], 
then U leaves both D{B) and D{F) invariant, and so B and F are distinct extensions 
with UBU~^ = —B and UFU~^ = —F. Their spectral measures thus yield distinct 
symmetric extensions. D 

Remark. Since U{ijj + e'^^Ctp) = e'^^{ijj + e~^^Cip), we see that if 6* 7^ 0, tt, then U does not 
leave D{B0) invariant. Thus among von Neumann solutions of the F problem, exactly 
two are symmetric. 

As a consequence of our identification of A* , we can prove one small part of Theorem 3. 
Recall that 7r(xo) is the sequence {-Pn(3^o)}5^o ^^^ Ci^o) is the sequence {<5n(a^o)}5^i- 

Theorem 2.14. Let {'^n}'^=o be a set of Hamburger moments. If either of the following 
holds, then the Hamburger moment problem is indeterminate. 

(i) For some xq E M, 7r(xo) and f^(a^o) H^ in £2. 
(ii) For some xq E M, ^{xq) and -^{xq) lie in £2- 

Remark. Since Pniz) is a polynomial, ^ makes sense as a sequence. 
Proof. In terms of the formal adjoint JF, we have 



J^{n{xo)) =xo7r(xo) 



d 



TV 



^{^{xo)) = So + a;o^(xo) 
^f o- (3^0) j = ^(a^o) + Xq -— . 



(2.19) 
(2.20) 
(2.21) 
(2.22) 



Thus, if (i) holds, we conclude tt, ^ are in D{A*) and 



{A* -xo)tv{xo)=0, {A*-xo)-^ 



n{xo). 



x=xo 



If A* were self-adjoint, then 

lk(xo)f = 



/ 4 * N ^^ 

TV, [A -xo)-^ 



x=xo 



/ 4* ^ ^^ 

(A -a:o)7r, — 



X = Xo 



which is impossible since Po{xo) = 1. Thus, A* is not self-adjoint and the problem is 
indeterminate. 
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If (ii) holds, t H are in D{A*) and 



{A* - xo)^{xo) = 5o, {A*-xoj 






If A* were self-adjoint, then 

\\axo)f = { axo)AA*-xo) 



{A* - xo)^(xo 
dQo{x 



dx 
dx 



x=xo 



x=xo 



dx 



-0 

x=xo 



since Qo{xo) = 0. This is impossible since Qi{xo) ^ and so again, the problem is 
indeterminate. D 

While we have used the theory of essential self-adjointness to study the moment prob- 
lem, following Nussbaum [30], one can turn the analysis around: 

Definition. Let i? be a symmetric operator on a Hilbert space 7Y and let </? G C^^B) = 
r\nD{B^). (f is called a vector of uniqueness for B if and only if the Hamburger moment 
problem for 

(yp, i?» 

^" = II 112 {Z.Z6) 

\m\ 

is determinate. 

Theorem 2.15 (Nussbaum [30]). Let B he a {densely defined) symmetric operator on 
a Hilbert space Ti. Suppose that D{B) contains a dense subset of vectors of uniqueness. 
Then B is essentially self-adjoint. 

Proof. Let U be the vectors of uniqueness for B. Let ip eU and let A^ be the restriction 
of B to the closure of the space spanned by {B'^ip}'^^^. A^ is unitarily equivalent to 
the Jacobi matrix for the moment problem (2.23). Thus, A<^ is essentially self-adjoint 
on Hip, the closure of the span of {B^(p}'^^Q. Therefore, {A^ + i)[D{Aip)] = Hip and. 



in particular, ip G {A^ + i)[D{A^)] C {B + i)[D{B)]. It follows that U C {B + i)[D{B)] 
and thus, Ran(i? + i) = H. Similarly, Ran(i? — i) = H, so B is essentially self-adjoint 
by Corollary 2.2. D 

The motivation for Nussbaum was understanding a result of Nelson [28] . 

Definition. Let S be a symmetric operator. </? G C°° is called 

analytic if Jl-B""!/?!! < CR^nl 
semi-analytic if ||i?"^(/?|| < CR^{2n)l 
quasi- analytic if Xl^i ||-B^"'</7||"^/^"' < oo. 
Stieltjes if ^^=1 WB^'pf/'^'' < oo. 
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Corollary 2.16. If B is a symmetric operator, then any analytic or quasi- analytic vector 
is a vector of uniqueness. If B is also hounded from below, then any semi-analytic or 
Stieltjes vector is a vector of uniqueness. In particular, 

(i) (Nelson [28]) If D{B) contains a dense set of analytic vectors, then B is essen- 
tially self-adjoint. 

(ii) (Nussbaum [31], Masson-McClary [23]) If D{B) contains a dense set of semi- 
analytic vectors and is bounded from below, then B is essentially self-adjoint. 

(iii) (Nussbaum [30]) If D{B) contains a dense set of quasi- analytic vectors, then B 
is essentially self-adjoint. 

(iv) (Nussbaum [31], Masson-McClary [23]) If D{B) contains a dense set of Stieltjes 
vectors and is bounded from below, then B is essentially self-adjoint. 

Proof. The first two assertions follow Proposition 1.5, Corollary 3.4, and Corollary 4.5. 
Theorem 2.14 completes the proof. D 



§3. The Stieltjes Moment Problem as a Self-Adjointness Problem 

One goal in this section is to prove the remaining part of Theorem 2, namely, that if the 
operator A has a unique non-negative self-adjoint extension, then the Stieltjes problem 
is determinate. In the indeterminate case, we will also introduce two distinguished non- 
negative von Neumann solutions: the Friedrichs and Krein solutions. 

The name Friedrichs and Krein for the associated self-adjoint extensions is standard. 
We will naturally carry over the names to the measures dfip and d/iK that solve the 
Stieltjes problem. Those names are not standard. Krein did important work on the 
moment problem, but not in the context of what we will call the Krein solution because 
of his work on self-adjoint extensions. 

We begin with a very brief summary of the Birman-Krein-Vishik theory of extensions 
of strictly positive operators (see [3] for a complete exposition). Suppose A is a closed 
symmetric operator and for some a > 0, (</?, A(f) > a||(/'|P for all </? G D{A). One closes 
the quadratic form analogously to closing the operator, that is, </? G Q{Af) if and only 
if there is (pn G D{A) so that </?„—> </? and {^pn — V'm, ^(Vn — V'm)) -^ as n — > oo. One 
then sets (</?, Ai?(/?) = lira^LpnjAipn)- The quadratic form (</?, Ai?(/?) on Q{Af) is closed, 
that is, Q{Af) is complete in the norm \/{(p, Ap(p). 

Not every quadratic form abstractly defined is closable (e.g., on L^{M.,dx), the form 
with form domain Co^(]R) given by (/? i-^ |(/7(0)p is not closable), but it is not hard to 
see that if a quadratic form comes from a non-negative symmetric operator, it is closable 
(see Theorem X.23 in [34] or [12], Chapter VI). 

While a closed symmetric operator need not be self-adjoint, it is a deep result (see 
[12, 34]) that a closed positive quadratic form is always the quadratic form of a self- 
adjoint operator. It follows that Ap is the quadratic form of a self-adjoint operator, the 
Friedrichs extension of A. It follows that inf spec(Ai;') = a > 0. 

Define N = ker{A*) which is closed since A* is a closed operator. Then D{A*) = 
D{A) + A^ + Ap [N], where -I- means independent sum. In particular, A^ ^ {0} if A is 
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not self-adjoint (which we suppose in the rest of this paragraph). The Krein extension is 
the one with D{Ak) = D{A) + N. Then Q{Ak) = QiAp) + N. Since A^ C D{Ak): we 
have that inf spec{AK) = so Ak 7^ Ap- The set of non-negative self-adjoint extensions 
of A is precisely the set of self-adjoint operators B with Ak < B < Ap (where Q < C < D 
if and only if Q{D) C Q{C) and (</?, C*/?) < {ifi^Dip) for all (f G Q{D); equivalently, if 
and only if (L> -|- 1)~^ < (C + 1)~^ as bounded operators). These S's can be completely 
described in terms of closed but not necessarily densely defined quadratic forms, C, on 
A^. One can think of C = 00 on Q{C)^. Ap corresponds to C = 00 (i.e., Q{C) = {0}) 
and Ak to C = 0. In general, Q{Ac) = Q{Af) +Q{C) and if v? G Q{Af) and r] G Q{C), 
then {ip + ?7, Aci(p + v)) = (</?, ^f^) + (??, Cr]). 

Two aspects of Jacobi matrices (defined on sequences of finite support) make the 
theory special. First, dim(A^) < 1. Second, by Theorem 5 (which we will prove in 
the next section without using any results from the present section), if A is not essen- 
tially self-adjoint, Ap has a discrete spectrum with the property that if t G spec(Ai?), 
then fXF{{t}) > /u({t}) for any other solution, //, of the Hamburger moment problem. 
Moreover, if t G spec^Ap), t ^ spec(i?) for any other self-adjoint extension, B. 

Proposition 3.1. Let A be a Jacobi matrix defined on the sequences of finite support, 
and suppose that A has deficiency indices (1,1) and that A is non-negative {i.e., that 
the moment problem associated to A has at least some solutions of the Stieltjes moment 
problem). Then A has a unique non-negative self-adjoint extension if and only ifO is an 
eigenvalue in the spectrum of Ap. 

Proof. Let a = inf spec{Ap). By hypothesis, a > 0, and by Theorem 5, ct is a discrete, 
simple eigenvalue of Ap. We want to show that if ct > 0, then A has additional non- 
negative self-adjoint extensions; and contrary-wise, if ct = 0, then it has no additional 
non-negative self-adjoint extensions. 

Ifa>0,A>Q;>0, then A has a Krein extension distinct from the Ap, so there are 
multiple non-negative self-adjoint extensions. 

Suppose a = and that B is some non-negative self-adjoint extension. Since Ap is 
the largest such extension, < S < Ap, so 

(Af + 1)-'<(S + 1)-^<1. (3.1) 

Since a = 0, \\{Ap + l)-^\\ = 1, so (3.1) implies that \\{B + l)-^\\ = 1. Since B has 
discrete spectrum by Theorem 5, this means is an eigenvalue of B. But, by Theorem 5 
again, Ap is the unique self-adjoint spectrum with a zero eigenvalue. Thus, B = Ap, 
that is, there is a unique non-negative self-adjoint extension. D 

Theorem 3.2 (last quarter of Theorem 2). Suppose that A is a Jacobi matrix on 
the vectors of finite support and that A is non-negative. If A has a unique non-negative 
self-adjoint extension, then the associated Stieltjes moment problem is determinate. 

Proof. Clearly, if A is essentially self-adjoint, then the Hamburger problem is determinate 
and so, a fortiori, the Stieltjes problem is determinate. Thus, we need only consider the 
case where A has multiple self-adjoint extensions, but only one that is non-negative. 
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Then A has deficiency indices (1,1). By Proposition 3.1, a = inf spec(Air) = 0. 
Moreover, by Theorem 5, is a pure point of fip- Let (p = P{o}(^-f)'^o. By Lemma 2.8, 
(/? 7^ 0. Then there exists (pn ^ D{A) so i/?^ ^ </? with Hv'nll = llv^ll and {(fn, Apipn) -^ 
(AiT'V?, v?) = 0. If ;Uf({0}) = r, then ||(/7|p = (</?, 5o) = r so ||(/7n|| = Vr and ((/?„, 5o) ^ r. 

Suppose p is some solution of the Stieltjes moment problem. Since (fn ^!- (fin L'^^dup), 
(fn is Cauchy in L^((ip), and soipn ^^ f inL'^{dp) and ||/|| = a/t. Since j x(fin{xY dp — > 0, 
we conclude that fxpdp = 0. Since p is a measure supported on [0, oo), we conclude 
f{x) = for a; 7^ 0, and thus p has a pure point at zero also. Since {^Pm^o) -^ t, we 
see that f f dp = lira f ipri{x) dp{x) = lim((/7^,5o) = r. But J f^ dp — ||/|p = r and 
f f dp = T. Thus /(O) = 1 and so p({0}) = r. But Theorem 5 asserts that any solution 
of the Hamburger problem distinct from hf has /ii?({0}) < r. We conclude that p = hf, 
that is, that the solution of the Stieltjes problem is unique. D 

We will call the solution of the Stieltjes problem associated to the Friedrichs extension 
of A the Friedrichs solution. If the Stieltjes problem is indeterminate, then A > a > 
by Proposition 3.1 and Theorem 3.2. Moreover, since the Hamburger problem is 
indeterminate, A^ = ker(A) = span[7r(0)] has dimension 1, so the Krein extension exists 
and is distinct from the Friedrichs extension. We will call its spectral measure the Krein 
solution of the Stieltjes problem. 

As discussed in the introduction, every indeterminate Stieltjes problem has associated 
to it an example of a determinate Stieltjes problem where the associated Hamburger 
problem is indeterminate. For by Proposition 3.1 and Theorem 3.2, the indeterminate 
Stieltjes problem is associated to an A where the bottom of the spectrum of Ap is some 
/o > 0. The moment problem {7m(-/o)} with 7m(-/o) = ZlJLo {"j)lji-fo)'^~^ wiU 
be determinate Stieltjes (since the associated Ap is the old Ap — /o) but indeterminate 
Hamburger. 

We summarize with 

Theorem 3.3. Suppose that {7n}5^o ^-^ ^ ^^^ of Stieltjes moments. Then, if •^ is Ham- 
burger indeterminate, there is a unique cq < so that 7(co) is Stieltjes determinate. 
Moreover, 7(c) is Stieltjes indeterminate if c > cq and 7(c) are not Stieltjes moments 
if c < Co- In particular, if 'y{c) is Stieltjes determinate for two values of c, then 7 is 
Hamburger determinate. 

Corollary 3.4. // {7n}5^o ^•^ '^ ^^^ ^f Stieltjes moments and obeys 

|7n| < CR^{2n)\ (3.2) 

then {7n}5^o ^•^ Hamburger determinate. 
Proof. By (3.2), 

l7n(c)| < CR''{2n)lf2(f)c' =C[R{c + l)r{2ny. 
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SO by Proposition 1.5, {7n}5^o ^^ Stieltjes determinate for all c > 0. By Theorem 3.3, it 
is Hamburger determinate. 



§4. Transfer Matrices and the Indeterminate Moment Problem 

In this section, we will analyze the indeterminate Hamburger problem in great de- 
tail, using the second-order difference equation associated to the Jacobi matrix A. In 
particular, we will prove Theorems 3, 4, and 5. 

Throughout, we will fix a set of Hamburger moments {'Jm}m=o- ^^ ^^^^ t)e convenient 
to define a linear functional E{-) on polynomials P{X) by E{X^) = 7^. Where we 
have polynomials P of several variables X,Y, . . . , we will indicate by Ex ( ■ ) the obvious 
map to polynomials of y, . . . with 

i?x(x"^i""^...) = 7^y^^... 

In (1.14), we defined the set of polynomials Pn{X) by the three conditions, deg(P^(X)) = 

n, E{Pn{X)Pm{X)) = Snm, and Pn{X) = CnnX"- H with Cnn > 0. We noted that 

with {a^}5^05 {^n}5^05 ^^^ '^-1 — 1' ^^^ elements of the Jacobi matrix associated to 7, 
the sequence Un = Pn{z); n = 0, 1, 2, . . . obeys 

ttnUn+l + {bn - z)Un + ttn-lUn-l = (4.1) 

with 

w_i = 0, uo = 1. (4.2) 

We also introduced the sequence of polynomials Qn{X) of degree n — 1 by requiring 
that Un = Qn{z) solve (4.1) with the initial conditions 

u-i = -1, uq = 0. (4.3) 

Notice that the WronskianIi^(g.(z),P.(2))(-l) = a-i[Qo{z)P-i{z) -Q-i{z)Po{z)] = 
1. Thus, by (2.11) we have: 

Proposition 4.1. 

ak-i[Qk{z)Pk-i{z) -Qk-i{z)Pk{z)] = 1 

The following is a useful formula for the Q's in terms of the P's. Note that since 
'^x-Y) ~ X]?=o X^Y'^~^~^ for any polynomial P, x-y ^^ ^ polynomial in X and 
Y. 



Theorem 4.2. For n > 0, 



gYi ''"'^j:y""^' i=Q„(n 
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Proof. ^ " (x_Y) is a polynomial in X and Y of degree n — 1, so we can define 
polynomials 



X -Y 
Subtract the equations, 

a^P^+i{X) +b^P^{X) + a^-iPn-iiX) = XP^iX) 

ar,Pr,+ l{Y)+b^Pn{Y)+a^-lP^-l{Y)=YPr,{Y) 

and see that for n = 0, 1, 2, . . . with R_i{Y) = 0: 

'XP^{X)-YP^{Y) 



anRn+i{Y) + bnRn{Y) + an-lRn-l{Y) = Ex 



X-Y 



^ fYPJX)-YPJY) ^ ,,,, 
= Ex( x-Y + ^- ( ^ ) 

= YR^{Y)+5no. 

Replace the abstract variable Y with a complex variable z and let Rn{z) = Rn{z) if 
n > but R-i{z) = -1. Then 

anRn+l{z) + bnRn{z) + an-lRn-l{z) = zRn{z) 

(since the 5^0 term becomes —an-iRn-i{z) for n = 0). Moreover, 

R.,{z) = -1, ^o(^) = Ex (y3y) = ° 

so by uniqueness of solutions of (4.1) with prescribed initial conditions, Rn{z) = Qn{z), 
that is, for n > 0, Rn = Qm as claimed. D 

One consequence of Theorem 4.2 is that Pn{X) and QniX) have the same leading 
coefficient, that is, if Pn{X) = CnnX^ -\ , then Qn{X) = CnnX'^~^ + ■ ■ ■ . 

The following calculation is part of the standard theory [1] and plays a critical role in 
what follows: 

Theorem 4.3. Let p solve the moment problem and assume z G C\]R. Set ( = Gp{z) = 
/fg. Then 

{{X - z)-\P^{x))L2idp) = Quiz) + CPniz). (4.4) 

In particular, 

f^\Q^(z) + CPn{z)\'<P^ (4.5) 

^-^ imz 

n=0 
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with equality if p is a von Neumann solution. 

Remark. We will eventually see (Proposition 4.15 and its proof) that the "if in the last 
sentence can be replaced by "if and only if." 



Proof. 



^'^^'^^ dp{x) = CPniz) + f ^"^"^^ ^"^^^ dp{x) = CPniz) + Qn{z) 



X — Z / X — z 



by Theorem 4.2. This is just (4.4). (4.5) is just Parseval's inequality for the orthonormal 
set {Pn{x)}^^Q in L'^{dp) if we note that 



dp{x) 1 , , . ^ 

1 H = — - / Mx) 

\x — zr z — z 



1 1 



X — Z X — z 



ImC 
Ivnz 



If p comes from a self-adjoint extension, by construction of IH}'^', {-Pn(a^)}5^o i^ an 
orthonormal basis for LF'{dp), so equality holds in (4.5). D 

Remark. In l-[^'^\ thought of as limits of polynomials, the vector with components 
Qn{z) + (Pn{z) is represented by 

J2[Qn{z)+CPniz)]Pn{X). 
n 

This product of P's may seem strange until one realizes that Ylin=i-Pn{y)Pn{x) is a 
reproducing kernel in L'^{dp) for polynomials of degree A^. This links the construction 
to ideas in Landau [19]. 

Given Theorem 2, the following proves the equivalence of parts (i), (ii), (iii), and (vii) 
for z G C\R of Theorem 3. 

Proposition 4.4. Suppose that zq G C\M. Then the following are equivalent: 
(i) The Jacobi matrix A is not essentially self-adjoint. 
(ii) 7v{zo) = {Pn{zo)}?^=o isine. 
(iii) ^(^o) = {Qn(^o)}^=o ^sme. 
(iv) Both TV (zq) and ^(zq) are in l'^ . 

Moreover, when any of these conditions holds, there is a closed disk of positive radius, 
V^zq) , in the same half-plane as zq so that for any solution p of the moment problem, 
C = Gp^zo) G V{zq). The values of ( when p is a von Neumann solution lie on dV{zo) 
and fill this entire boundary. Every point in T>{zo) is the value of Gp{zo) for some p 
solving the moment problem. 

In addition, if these conditions hold, 7r(zo) 077.^^(^0) H^ in D{A*)\D{A). 

Remarks. 1. One can show (using (2.11)) that the center of the disk is 



J. Qn{zo) Pn-l{zo) - Qn-l{zo) Pn{zo] 



n — >oo 



Pn{zG)Pn-l{zo) — Pn-l(^o) -Pn(^o) 
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and the radius is 



^ ^ (4.6) 



but these exphcit formulae will play no role in our discussions. 

2. We will see later (Theorem 4.14) that if p is not a von Neumann solution, then 

Proof. Since A is real, A fails to be essentially self-adjoint if and only if there is a non-zero 
solution of {A* — zq)u = 0. By Theorem 2.6 and the unique solubility of second-order 
difference operators given U-i = and -uq, every such solution has Un = uoPn{zo) 
so (i) is equivalent to (ii). Let ( = Gp{zo) for some von Neumann solution. Then 
C{zq) + ('n'{zo) G £^, so (ii) is equivalent to (iii) or (iv). 
If (i)-(iv) hold, then (4.5) has the form 

a\Cf + bC + bC + c<0, 

where a = ||7r(2o)|P, c = ||^(-2o)|p7 and b = 2{^{zo),7t{zo)) — 2imz • "^^^ ^^^ where this 
inequality holds is always a disk, "^(-20) 7 although a priori, it could be empty depending 
on the values of a, 6, c. However, by Theorem 2.11, we know that there are multiple C's 
obeying (4.5) so the disk must have strictly positive radius. By Theorem 4.3, C's for 
von Neumann solutions obey equality in (4.5) and so lie in dV{zo). 

The self-adjoint extensions are parametrized by a circle (unitaries from /C+ to /C_) 
in such a way that the map from the parametrization to values of {5o, {Bf — zo)~^6o) is 
continuous. By Theorem 2.11, this map is one-one. Every one-one continuous map of 
the circle to itself is surjective, so all of dV{zo) occurs. 

Given a point C in ^(^o)'"*, find Co and Ci in dV{zo) and 9 E (0, 1) so C = ^Co + (l-^)Ci- 
If ^Oj A*i are the von Neumann solutions that have G^^(zo) = G (^ = 0, 1), then p = 
9^0 -|- (1 — 9)fXi is a measure solving the moment problem with Gp{zo) = C- 

If (iv) holds, then 7r(zo), ^{zq) E D{A*). By (2.14) and Proposition 1.4, {A*7i{zo), C{zo)) 
— {tv{zo), A*^{zo)) = —1, so neither w{zo) nor ^{zq) lies in D{A). D 

Remark. The Cayley transform way of looking at self-adjoint extensions says that for 

any zq E C\M, 

{Bt - zo){Bt - zo)-' = V{zo) + e^^(^0'*)X(^o), 



where V{zo) is the t-independent map {A—zo){A—zo)~^ from Ran(A — zq) to Ran(A — zq) 
and X{zo) is any isometry from ker{A* — zq) to ker(A* — zq) extended to 7i as a partial 
isometry. 6* is a parametrization depending on zo^t (and the choice of X{zo)), but the 
theory guarantees that as t varies through all possible values, e runs through the full 
circle in an injective manner. Now (Bt — zo)(Bt — ^o)"^ = 1 + 2z(Im2o)(-St — ^o)"^, so 

i5oABt-zo)-'6o) = [2i{lmzo)]-'[-l + {6o.V{zo)6o)+e''^'^-'\5o,X{zo)6o)] 
is seen directly to be a circle. Since one can take ^(^o) to be 

fX Vt^ \X \ — Prn{Zo)Pn{Zo) 
[0n,^[Z0 0^ — =35 — --^ 

Ej=ol^j(^o)r 

and Po{z) = 1, we even see that the radius of the circle is given by (4.6). 
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Corollary 4.5. // a Hamburger problem is indeterminate, then Xl^o I'^nl""'^ < oo. In 
particular, ifYl^=o \C'n\~'^ = oo, {e.g., an = 1); then the Hamburger problem is determi- 
nate. If 

CXD 

E-l/2n 
l2n = OO: 

n=l 

then the Hamburger problem is determinate. If 



E7n^^'" = «^ 



n=l 



for a set of Stieltjes moments, that problem is both Stieltjes and Hamburger determinate. 

Remarks. 1. The last pair of assertions is called Carleman's criterion. It generalizes 
Proposition 1.5 and proves uniqueness in some cases where growth doesn't imply analyt- 
icity of the Fourier transform. 

2. See Corollary 5.24 for the Stieltjes analog of the first assertion. 

Proof. By Proposition 4.1 and the Schwarz inequality, 

N /N+1 xl/2/iV+l X 1/2 

Y,\an\-'<2iY,\Pr.{z)n {Y.\Qn{z)n . 

n=0 ^n=0 ^ ^ n=0 ^ 

Thus, divergence of X^^o I'^nl"^ implies that either vr or ^ (or both) fail to be £2, and 
so determinacy 

Consider Carleman's criterion in the Hamburger case. By induction in n and (4.1) 
starting with Pq{x) = 1, we see that 

Pn{x) = (tti . . . an)~^x'^ + lower order, 

so ((ai . . .an)~^x^,Pn{x)) = 1 and thus, by the Schwartz inequality, 

1 < 72n(ai---ar^~^ 
hence 



"nj 



72^ 



^ S \ai. . .On) 



Therefore our result follows from the divergence criteria proven at the start of the proof 
and the inequality 

n n 

J2iai...a,)-'/^<2eJ2a-' 
(which is actually known to hold with e in place of 2e) . 
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To prove this, note first that 1+x < e^ so {1+n ^)^ < e so using induction, n^ < e"n! 
so by the geometric-arithmetic mean inequahty, 

{a^... a,)-V^- = {a^'2a^' . . .jaj')'/^ {j\)-'/^ 



J 

J 1 — 1 



Thus, 



3 = 1 fc=l S=fc J = l 



Since 



oo 



k ^l , f°^ dy 



y—<- + k ^<2. 

By Proposition 1.6, Carleman's criterion in the Hamburger case means that if a set of 
Stieltjes moments obeys X]^i(7n)~^^^"^ = oo, then it is Stieltjes determinate. Thus by 
Theorem 3.3, it suffices to show that 



J2^1n{c))-'/'->D{c)J2ln'^ 



L/2r 
In 

n=l n=l 



for c > and D{c) < cx) to conclude Hamburger deter minacy. 

Note first that ^^^^^^ > -^^ by the Schwarz inequality. Thus, if a = — , we have that 

In- 3 < InOl^ ■ 



Therefore, 



and so 



7n(c) <ln^y.\ c^a^ = (1 + ca)'^-fr, 

j=0 ^^^ 



n=l n=l 

To complete the proof of Theorem 3, we need to show that if every solution of [{A* — 
zo)u]n = (n > 1) is £^ for a fixed z = zq, the same is true for all 2 G C and that ^, 
^ are in £^. We will do this by the standard method of variation of parameters. So we 
write a general solution of [{A* — z)u\n = (n > 1) in terms of Pn(-2o) and Qn{zo) using 

""+'). a„f"«<f)+/3„(*;«<f). (4,7) 

Un J V Pn[Zo) J V Qn[Zo) J 
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Since W{P.{zq)j Q.{zq)) = —1, the two vectors on the right of (4.7) are hnearly indepen- 
dent, and so they span C^ and {an, (3n) exist and are unique. Indeed, by the Wronskian 
relation: 



an = W{Q.{zo),u){n) (4. 

(3n = -W{Rizo),u){n). (4.8b) 

A straightforward calculation using the fact that P,Q obey (4.1) with z = zq, (2.11), 
and (4.8) shows that (4.1) at n is equivalent to 

2^=[l + {z-zo)S{n,zo)](^2~'y (4.9) 



where 



q(^ y \ — I -Qn{Zo)Pn{Zo) -Qn{Zo)Qn{zo) \ (A^^\') 

''^'"'^O^ ^ Pn{zo)Pn{zo) Pn{zo)Qn{zo) '' ^ ^ 



For example, using (2.11) and (4.7) for n ^ n — 1: 

ttn = W{Q.{zo),u){n) 

= W{Q.{zo),u){n- 1) + {zo - z)Qn{zo)un 

= ttn-l - {Z- Zo)Qn{zo){an-lPn{zo) + (3n-lQn{zo)} ■ 

Notice that 

S^ = det{S) = Tr{S) = 0. (4.11) 

The following is obvious and explains why £^ solutions are so natural: 
If nizo) = {Pn{zo)}^=o and ^(^o) = {Qn(^o)}^=o are both in l^, then 

oo 

Y,\\S{n,zo)\\<^. (4.12) 

n=0 

Lemma 4.6. Let A^ be a sequence of matrices with Yl'^=o \\^n\\ < oo. Let Dj^{z) = 
(1 + zAj^){l + zA]^-i) . . . (1 + zAq). Then Doo{z) = UuiDn^z) exists for each z E C 
and defines an entire function of z obeying 

||^oo(^)|| < Ceexp{e\z\) 

for each e > 0. 

Proof. Notice first that 

11(1 + B^)...{1 + Bo)\\ < n(l + 115.11) < exp(5^ ||S,|| ) (4.13) 
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and that 



1\ / iV \ 

11(1 + S^) . . . (1 + So) - 111 < n(l + ll^.ll) - 1 < exp(5^ lis, II j - 1. 
From (4.14), we see that 

r / ^^ \ 1 

\DNiz)\\ 



(4.14) 



D 



N+j 



:z)-dn(z)\\< 



< 



N+j 

N+l 

oo 

expi \z\ J2 \\^j\ 

N+l 



N 

expi I^I^IIAjii p 



from which it follows that Dn{z) is Cauchy uniformly for z in balls of C Thus, Dc 
exists and is entire in z. By (4.13), 

TV ^ oo X 

poo(^)ii<n(i+i^iii^j-ii)^^p(i^i E ii^^-iij' 



TV 



SO given e, choose A^ so J27Ln+i ll^ill — f ^'^'^ ^^^ ^^^ ^^'^^ ^^^^ ^^^ polymonial ni=o(-'- + 
I^IIIAjll) can be bounded by Cg exp(ie|2|). D 

Given Proposition 4.4 and Theorem 2.14, the following completes the proof of Theo- 
rem 3: 

Theorem 4.7. Let A be a Jacobi matrix and consider solutions of (4.1) for z E C If 
for some zq G C, all solutions are in £'^ , then that is true for each z E C Moreover, tv{z) 
and ^{z) are analytic i2-valued functions so, in particular, ^, gf G £2 for all z. The 
disk 'D{zo) varies continuously as z runs through C_|_. 

Proof. Define 

T{n, -1; z, Zq) = {1 + {z- zo)S{n, ^o)) • • • (1 + (^ - ^o)5'(0, ^o))- 
By (4.12) and Lemma 4.6, 



sup||T(n, -1,;2, 2o)|| < expf \z - zo\Y\\S{j, z)\\ ) 

V ,=0 / 



< 00. 



Thus for any initial (^ ^), if ti has the form (4.7), then sup^ || (^") || < 
sup„ ||T(?i, -1;^, 2o)||||(^:^')|| = C. Thus, u^ = anPn{zo) + (3nQn{zo) has \Un\'^ < 
C'^[Pn{zo)'^+Qnizo)'^] whichisiu^i. Since n{z) is associated to T{n, —1; z, zq) (q) and ^(2) 
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is associated to T{n, —1; z, zq) (-^), the claimed analyticity holds because sup„ \\T{n, —1; z, zq) 
is bounded as z varies in bounded sets. 

Continuity ofV{zo) follows from the formula for dV{zo), viz. w G dV{zo) if and only 
if 

\\7r{zoW\w\^ + 2Re\(2{C{zo),7v{zo)) - wt^—]w\ + ||^(^o)f = 0. D 



By Theorem 3 and Lemma 4.6, if the Hamburger problem is indeterminate, 

T(oo, — 1; z, Zo) = lim T(?i, — 1; 2, ^o) (4-15) 

n — >cxD 

exists. We define four functions A{z), B{z), C{z), D{z) by 

r(.=,-l;.,.„ = 0)=(-^^(~) -^). (4.16) 

The Nevanlinna matrix is defined by 

Theorem 4.8. 

(i) Each of the functions A, B, C , D is an entire function obeying 



for each e > 0. 
(ii) For z small, 

and 

with a > 0. 
(iii) AD-BC=1. 



\f{z)\ < Csexp{e\z\) 

B{z) = -l + 0{z) (4.18a) 

D{z)=az + 0{z^) (4.18b) 



Proof, (i) follows from Lemma 4.6. (ii) follows if we note that 

oo 

T(oo, -l;z,zo = 0) = l + zY^ S{n, zq = 0) + 0{z^). 

n=0 

(4.18a) is immediate and (4.18b) follows since (4.10) says that a = J2^=o ^■"■{0)'^ > 0. 
(iii) holds since by (4.11), detT = 1. D 

For our purposes, the formula for A, B, C, D as matrix elements of an infinite product 
is sufficient, but there is an infinite sum formula connected to standard equations for 
perturbed solutions (see, e.g., [32, 11]) that lets us make contact with the more usual 
definitions [1]. 
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Theorem 4.9. If the Hamburger problem is indeterminate, then, 

oo 
A{z)=zY,Qn{0)Qr,{z) 
n=0 

CX) 

B{Z) = -1 + zJ2Qn{0)Pn{z) 

n=0 

oo 

C{z)=l + zJ2Pn{0)Qn{z) 
n=0 

oo 

D{z)=zJ2Pn{0)Pn{z). 
n=0 

Proof. 

T{n, -1; z, ^0 = 0) = (1 + zS{n, ^o = 0)) . . . (1 + ^^(0, zq = 0)) 

^l + J2 ^^(^^ -^0 = 0) n (1 + ^Sik, zo = 0)) 

j=0 k=0 

SO it suffices to sfiow tliat 

S{j, zo = 0)(1 + zS{j - 1, zo = 0)) . . . (1 + zS{0, zo = 0)) 

-Q,iO)P,{z) -Q^{0)Q^(z) 
PAO)P,{z) Pj{0)Q,{z) 



By definition, 



(1) (2) 



wfiere 



(1 + zS{j - 1, ^0 = 0)) . . . (1 + zS{0, ^0 = 0)) = ( J-/ {-^ 

Pj-i Pj-i 



(1) ( PM V .(1) f QM \^( PA^) 



(4.19) 



and 

^(2) ^^.(0) V .(2) /g,(0) \^(QA^) 

^-1 \p,-m) ^'-' \Qj-m) \Qj-M 

Tlius by (4.10), 

LHS of (4 19) = (^ -^^•(0)[4-i^:^(0) + (^f-iQM\ -Q,(0)[af_\P;(0) + /3f_W,(0)] 

\pjm4\pM+^f\Q3m Pjmc^f\pM + p^"\Qr'' 

= RHSof (4.19). D 
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7r(0) = {Pn(0)} and ^(0) = {Qn{0)} span D{A*)/D{A) by Proposition 4.4. Notice 
that A*n{0) = while A*^{0) = 5o, so (7r(0), A*^{0)) = 1 and (2.7) holds. Thus, we can 
parametrize the self-adjoint extensions Bt of A by 

D{Bt) = D{A) + {a{tn{0) + ^(0)) | a G C} (4.20a) 

if t < oo and 

D{B^)=D{A) + {an{0)}. (4.20b) 

This is equivalent to defining t by 

t = {6o,Br'6o). (4.20c) 

Theorem 4.10. For each t G M U {oo} and z G C\R, 

(5o,(5,-.) M = -^(^)^^5(^)- (4.21) 

Proof. Let us sketch the proof before giving details: 

T(-l, oo; z, zo = Q)= T(oo, -1; z, zq = 0)"^ (4.22) 

Then 

^]|) =T(-l,oo;^,zo = 0)Q (4.23) 

is such that the u^ associated to (a^^) obeys (4.1) and is asymptotically t7r(0) +^(0) and 

so in D{Bt). {Bt - z)-Ho wiU thus be -^^^ and (5o, {Bt - z)-Ho) wiU be -^^y 
But 'Ut(O) = a_i and «*(— 1) = — /9-i so {5q, Bt — z)~^5q) will be ^^, which is given by 
(4.21). 

Here are the details. If </? G D{A*), then by (2.14), </? G -D(-Bt) if and only if 

lim W((/:,t7r(0) + ^(0))(n) = 0. (4.24) 

n — >oo 

Suppose u solves (4.1) for a given z. Then u G D{A*) and u has the form (4.7), where 
by (4.15), Hm^^oo (J^;;) = (J^) exists. Clearly, 

W{u, t7r(0) + ^(0))(?i) = -a^ + /?^t 

so (4.24) holds if and only if 

Thus, if T(-l,oo;z,zo = 0) is given by (4.22), then g-^) given by (4.23) is initial 
data for a solution u of (4.1) that has u G D{Bt). A solution u of (4.1) has 

{A* — z)u — — w_i5o 
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and thus, if u is associated to the data in (4.23), 

(Bt - z)u = -u-i5o 



and 



But 



so 



{do,iBt-z)-'6o) 



Uq 



U-i 






(So,(Bt-z)-'do) = ^. 

P-i 

Since T(oo, —1; z, zq = 0) has the form (4.16) and has determinant 1, its inverse has the 
form 

' C{z) A{z) 
-D{z) -B{z) 

and so a_i = C{z)t + A{z), /3_i = -D{z)t - B{z). Thus (4.21) is proven. D 

Remark. Our convention for the parameter t differs from that in Akhiezer [1], which is 
the traditional one! If s — —t~^, then 

Ct + A As-C 



Dt + B Bs-D 

The parameter s is what he caUs t. Our ^{z) later and his ^{z), which I will call '^{z), 
are related by ^{z) = —'^{z)~^. Since this transformation takes Herglotz functions to 
themselves, we both deal with Herglotz functions. See [9] for an interesting alternate 
reparametrization. 

We turn next to half of Theorem 5. We will prove that each Bt has point spectrum 
with any two disjoint. The condition that f^t{{i}) > pi{i}) for any other solution of the 
moment problem will wait until after we have proven Theorem 4. 

Theorem 4.11 (half of Theorem 5). Suppose the Hamburger problem is indetermi- 
nate. Each Bt has pure point spectrum only. The eigenvalues of the different Bt 's are 
distinct and every x eR is an eigenvalue of some Bt- If x is an eigenvalue of Bt, then 

Moreover, if Xn{Bt) are the eigenvalues of Bt {ordered in some convenient way), then 
for any p > I, 

J2\^n{Bt)\-P < ^. (4.26) 
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Proof. Let Ht be the spectral measure for Bf. By Theorem 4.10, 

d^t{x) _ C{z)t + A{z) 
X- z ~ ~ D{z)t + B{z) 

is a meromorphic function of z, since A, B, C, D are entire. That imphes that ^t is pure 
point only and these points are given by solutions of 

D{z)t + B{z) = 0. 

Since A, B, C, D are real when z = x is real, for any x, x E spec(i?t) for t = —jyt§\- Since 
AD — BC = 1, B and D have no simultaneous zeros, so the eigenvalues are distinct for 
different t's. (4.26) is a standard result on the zeros of an entire function like D{z)t+B{z), 
which obeys (i) of Theorem 4.8. 

To prove (4.25), note that if x is an eigenvalue of Si, then the corresponding normalized 
eigenvector ip obeys A*ip = xip. It follows that (fn = foPn{x) and then by normalization, 
that ifl = 1/ Er=o \Pn{x)\^. Thus, iJt{{x}) = ifl is given by (4.25). D 

Now define the map F{z) : C U {oo} ^ C U {oo} by 

D[z)w + B[z) 

F{z) as a fractional linear transformation is one-one and onto, and takes M U {oo} onto 
a circle or straight line. By Proposition 4.4 and (4.21), F(2)[]RU {oo}] is precisely the 
circle dV{z) when Im^ 7^ 0. 

Proposition 4.12. Iflmz > 0, F{z) maps the upper half-plane C+ = {z E C \ Imz > 
0} onto the interior of the diskV{zy^^ . 

Proof. For each z G C\M, F{z) must map C_|_ onto either ©(2)'"* or C\D{z) since M is 
mapped to dV{z) and a fractional linear transformation maps CU{oo} to itself bijectively. 
Thus it suffices to show F(2)[z] G 2^(2)™* for all z with Imz > 0. But F(^)[z] moves 
analytically, and so continuously. It can only move from V{zy^^ to C\D[z) by lying 
on dV{z), which is impossible if Im^ > since then dV{z) — F{z)\R U {00}]. Thus it 
suffices to show F(2)[z] G V{zy^^ for z = ix with a:: > and small. 

Now, F{z)\i] G 1^(2)'"* if and only if F{z) takes the lower half-plane to CU{oo}\V{z), 
which holds if and only if F{z)~^[oo] is in the lower half-plane. Thus it suffices to show 
for X small and positive that F{ix)~^[oo] G C_. 

F{ix)~^[oo] is that w with D{ix)w + B{ix) = 0, that is, w = —jyf^- But by (4.18), 
B{ix) = — 1 -|- 0{x) and D{ix) = iax + 0{x'^) with a > 0, so w = —ia~^x~^ + 0(1) lies 
in C_, which suffices to prove the proposition. D 

We now follow the classical argument of Nevanlinna [29] to obtain Theorem 4. We 
need to relate solutions of (1.1) to asymptotics of the Stieltjes transform of p following 
Hamburger [8] and Nevanlinna [29]. 
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Proposition 4.13. Let M^{-f) = {p \ p obeys (1.1)}. Let Gp{z) = J 
any N as y ^> oo, 



dp{x) 



N+1 



N 



G,{ty) + J2(- 



.^)n+l^-n-l^^ 



n=0 







Then for 



(4.28) 



uniformly for p inM-^i^). Conversely, if G is a function on C+ withImG{z) > there 
and so that (4.28) holds for each N, then G{z) = Gp{z) for some p in M.^{'^). 

Proof. By the geometric series witli remainder, ttie left side of (4.28) is 






so, using \x — iy\ ^ < y '^, 



\RNip)\<y-' / kr+Mp( 



[X 



X — ly y 



iN+iy~ 



(4.29) 



if N is odd 



<hhN + 7N+2]y ^ if A^ is even. 



This shows that (4.28) holds and the convergence is uniform for p G ^^^(7). 

For the converse, first use the Herglotz representative theorem which says that if G 
maps C-i- to C+, then for some measure dp, some c > and some real d: 



where, a priori, p only obeys 



1 



X 



X — z 1 + x^ 



(4.30) 



dp{x) 

1 + x^ 



< 00. 



By (4.28), 



y[G{iy)] -^ iio- 

If (4.30) holds, then y~^G{iy) — > zc, so c = 0. Since c = 0, (4.30) says 

,2 



(4.31) 



y lmG{iy) 



2 I 2 ^^(^) 
a;^ + y 



so (4.31) and the monotone convergence theorem implies that j dp{x) — 70. Once this 
is true, (4.30) implies that as y — * 00, 



Re G{iy) ^ d — j dp{x) — 



X 



+ X' 



(4.32) 



so (4.31) implies the right side of (4.32) is zero, and thus (4.30) becomes 

dp{x) 



G{z) 



X — z 



(4.33) 
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We will now prove inductively that p obeys (1.1), that is, p G 7Vf^(7). Suppose that 
we know that (1.1) holds forn = 0, 1, . . . , 2M - 2. (4.28) for N = 2M then implies that 

dp[x) + Z?/72M-1 -" -12M- 



X — ty 
Taking real and imaginary parts, we see that 

72M-1 = lim / ^ — 7^dp{x) (4.34) 

y^oo J X"^ + y^ 

72M = hm / ^- — :rdp{x). (4.35) 

y^oo J a;^ -\- y^ 

(4.35) and monotone convergence implies that f x'^^ dp{x) < oo, so that dominated 
convergence and (4.34), (4.35) imply that (1.1) holds for n = 0, 1, 2, . . . , 2M. D 

Let JF be the set of analytic maps from C+ to C+ U {oo}. By the open mapping 
theorem for analytic functions, $ G JF either maps C+ to C-|- or $ has a constant value 

inMU{cx)}. 

Theorem 4.14 (Nevanlinna [29]). Let 'j be a set of indeterminate Hamburger mo- 
ments. Then there is a one-one correspondence between solutions, p, of the moment 
problem and functions $ G JF given by 

^ fdp{x)_ C{zMz)+A{z) 
'- J x-z D{z)^z) + B{z) - ^ ' 

In particular, if p is not a von Neumann solution, Gp{z) G D{zy^^ for all z G C-|-. 

Remarks. 1. By (4.21), the von Neumann solutions pt correspond precisely to the case 
$(z) =t gMU{oo}. 

2. We will call the Herglotz function associated to some p G A4^ {'^) the Nevanlinna 
function of p and denote it by $p. Conversely, given $ (and 7), we denote by p$ the 
associated solution of the moment problem. 

3. In Section 6, we will discuss the p$'s associated to the case where the measure p 
in the Herglotz representation (1.19) is a measure with finite support. 

Proof. Let $ G JF and let G{z) denote the right side of (4.36). In terms of the map F{z) 
of (4.27), G{z) = F{z){^{z)), so by Proposition 4.12, gIz) G V{z). By the uniformity in 
Proposition 4.13, the fact that the Ait G M^ {-f) and that {G^^{z) \ t G MUJoo}} = dV{z), 
this implies that G{z) obeys (4.28). Thus by Proposition 4.13, G{z) = Gp{z) for some 
peM^{-f). 

Conversely, let p G Ai^ {•y) and consider Gp{z). By Theorem 4.3, Proposition 4.4 
and Proposition 4.12, $(2) = F{z)-^{Gp{z)) G C+ U {00}, and so $ G J^. Thus, 
Gp{z) = F{z){^{z)), that is, (4.36) holds. 

If p is not a von Neumann solution, then lin.^p{z) > 0, so Gp{z) = F{z){^p{z)) is in 
the interior of 1^(2). D 
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Proposition 4.15. Let p G A4^{'y). If p is a von Neumann solution, then the poly- 
nomials are dense in L'^(K,dp) and for each zq E C+, there is no other p in A4^ {'j) 
with Gfj,{zo) = Gp{zo). If p is not a von Neumann solution, then the polynomials are 
not dense in L'^{M.,dp) and for each zq G C+, there are distinct p's in A4^ {'^) with 
G^{zo) = Gp{zo). 

Proof. If p is a von Neumann solution, ttie polynomials are dense by construction and 
the uniqueness result holds since G^(-2o) = Gp{zo) if and only if $^(2:0) = ^p(^o) and 
^p{zo) is then real. 

Conversely, suppose p is not a von Neumann solution. Then, Gp{zo) G D(z)™*, so by 
the proof of Theorem 4.3, {x — zq)~^ is not in the L^ closure of the polynomials. Thus, 
the polynomials are not dense. If $p is not a constant, set "^(z) = —\^p{zo)\'^^p{z)~^ + 
2Re^p{zo). If^p{z) is a constant, q (of necessity with Img > 0), set "^(z) = cz + d with 
c = Y^^ (> 0) and d = Req — cRezo. In either case, \1/ 7^ $, \1/ is Herglotz, and "^{zq) = 
$(zo). Thus, if p = p*, we have Gp{zo) = F{zo){^{zo)) = F{zo){^{zo)) = Gp{zo). □ 

To complete the proof of Theorem 5, we need the following fact about Herglotz func- 
tions: 

Theorem 4.16. Let ^{z) be a Herglotz function so 

1 X 



*(.)=c. + d+/d,(x) 



X — z 1 + x^ 



with d real, c>0, and either c > or dp j^ 0. Suppose for some xq, ^{xq + ie) ^ t G 
Then either I '■ ^° . ^^' — - 1 ^ 00 or else f , {^ < 00 and 



^{xo+ie)—t C dp{x) 

slo ie J {x — Xq) 



lim "^-^"~_ ^~ =c+ / ,^:^^\":, . (4.37) 



Remark. Do not confuse the p in the Herglotz representation for $ and the measure p of 
(4.36). In particular, p is allowed to have finite support even though we are suppposing 
that p does not. 



Proof. Note first that 



On the other hand, 



liii^{xo + ie) , f dp{x 

{x — xo)^ + e 



c+l ,.. ""';;^^,. . (4.38) 



9$ , . , f dpix) 2e(x — Xq) 

^— {xo + ie)= / ., ^ ^ ,^ -^^ 

dz ^ 'J [{x-xoY+e' 



Im ^ (xo + ^e)= I Z} 'ZZIZ ' (^.39) 



If J t^_^ \2 = 00, then (4.38) implies Im [_L£o±l£i — ] ^ 00 , so if the limit is finite, then 
/ {x-x V ^ ^^' (4-39) and the dominated convergence theorem implies that Im ^{xq + 
ie) ^ so that [R^'^("o+^^)-*] = i JJini||(xo + iy) dy -^ 0. This and (4.38) implies 
(4.37). D 
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Theorem 4.17 (end of Theorem 5). Let 7 be an indeterminate moment problem. 
Let p E J^^ {'-f) correspond to a ^ which is not a constant in M U {00} {so that p is not 
a von Neumann solution). Suppose a = p{{xo}) > for some point xq G M. Then there 
is a von Neumann solution pt so that pt{{xo}) > a. 

Proof. We will suppose that D{xo) 7^ 0. If D{xo) = 0, then B{xo) 7^ and we can use 

C{z)^z)+A{z) _ A{z){-^{z))-^-Ciz) 
~ D{z)^z) + B{z) ~ ~ B{z){-^z))-^ - D{z) 

and ■~^{z)~^ in place of ^{z) to give an identical argument. Define t — ^^^ 

1 



Since AD - BC = 1, 



' D{xo) 



C{xo)t + A{xo) 



D{xo) 



(4.40) 



is non-zero and has the same sign as D. 

C{xo + ie)^{xo + is) + A{xq + is) 



{—ie)G p{xQ + ie) = (ie) 



a. 



t. Then 



D{xo + ie)^{xo + ie) + B{xo + ie 
This is only possible if Gp{xo + ie) ^ 00, which requires that $(xo + ie 

C{xo + ie)^{xo + ie) + A{xo + ie) -^ C{xo)t + A{xo) 

while 

D{xo + ie)^xo + ie)+B{xo + ie) , , ^xq + ie) - 1 
: > D {xo)t + B (xo) + Dixo) lim ■. 

te eio le 



Thus using (4.40), 



1 



a 



D{xo)' lim,^o ''^"tf^"' + D{xo)[D'ixo)t + B'{ 
On the other hand, taking $(x) =t to describe pt, 

1 



X] 



lJ't{{xo]) 



D{xo)[D'{x^)t + B'{xo)] 



must be positive. It follows from Lemma 4.16 that J 



d/i(x) 
{x-xaY 

a-^ > Pt{{xo})~^ . □ 



< 00 and that 



Finally, we turn to the Nevanlinna parametrization in the indeterminate Stieltjes case. 
As we have seen, among the Bt there are two distingtuished ones, Bp and Bk- The Krein 
extension has ker(A*) C D{Bk) so 7r(0) G D{Bk), which means t — 00 in the notation 
of (4.20). The Friedrichs extension is Bp — Bt^ where to = (5o, Ap 5o). 
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Theorem 4.18 (Stieltjes case of Theorem 4). // the Stieltjes problem for 7 is in- 
determinate, its solutions obey (4.36) for precisely those $ that either are a constant in 
[to,oo)U{oo} or that obey ^{z) is analytic in C\[0, 00) with^{x) G [to, 00) ifxE (— oo,0). 
Hereto = (5o,Ap^5o). 

Remark. Such $'s are precisely those of the form 

$(.)=c^+r^^^, (4.41) 

Jo X — z 

where d >to and J^ ~^+i < ^^• 

Proof. Clearly, the solutions of the Stieltjes problem are precisely the solutions p of the 
Hamburger problem with linisioGp{—y + is) real for each y G (0, 00). Since the Bts 
obey t ~ {5q, B^^Sq) and Bk < Bt < Bp means to < (^o, B^^6o) < 00, we know that for 
y G (-00, 0), F{y) maps [to, 00) into (0, 00) (F( ■ ) is given by (4.27)). It follows that if $ 
obeys $(a;) G [to, cxo) for x G (— cxo, 0), then Gp{y + is) has a limit in (0, cxo) on (—00, 0), 
that is, $ defines a solution of the Stieltjes moment problem. The converse follows from 
this argument and the next theorem. D 

Theorem 4.19. Let dfxp o,nd dpK be the Friedrichs and Krein solutions of an indeter- 
minate Stieltjes moment problem and let dp be another solution. Then for all y G (0, 00), 

dprix) ^ /•"" dp{x) ^ /•"" dpKJx) ,^^2) 



x + y Jo x + y Jo x-\-y 

This will be proven below as part of Theorem 5.2. 
Corollary 4.20. If p E M^i'j) and dp 7^ dpp, then for y G [0, 00 



dpF{x) ^ /""" dp{x) 



x + y Jo x + y 

{Note that the inequality is strict.) 

Proof. For each y G (0, 00), F{—y){ • ) is a strictly monotone map of [to, 00) to [G^p (—?/), 
Gpj^{—y)) by the proof of Theorem 4.18 and the fact that any fractional linear map 
of R to M is monotone in any region where it is finite. By a limiting argument, the 
same is true for y = 0. Thus, (4.43) follows from the relation of Nevanlinna functions 
$p(y) > ^ij.p{y) = to- This is immediate from (4.41). D 

As a corollary of this, we see 

Corollary 4.21. Let 7 6e a set of indeterminate Stieltjes moments and let dpp be its 
Friedrichs solution. Let r = f x~^dnF{x) and let 

70 = 1; 7n = r"Sn-i, n = l,2, 
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Then {7j}°^o ^-^ ^ determinate Hamburger moment problem. 

Remark. This lets us produce Stieltjes determinate moment problems with essentially ar- 
bitrary rates of growth for 7^. By using Theorem 2.13, we obtain determinate Hamburger 
problems with arbitrary fast growth also. 

Proof. r~^x~^dfXF = dXo{x) is a Stieltjes measure whose moments are 7. Let dX{x) be 
another such Stieltjes measure and let dp{x) = rxdX{x). Then the moments of dp are 
7 and j x~^dp{x) = r. So, by Corollary 4.20, p = pp- Thus, A = Aq. Thus, Aq is the 
unique solution of the Stieltjes problem. But ^ supp((iAo). Thus by Proposition 3.1, 
the Hamburger problem is also determinate. D 



§5. Pade Approximants, Continued Fractions, and Finite Matrix Approxima- 
tions 

In this section, we will primarily consider aspects of the Stieltjes moment problem. 
We will discuss finite matrix approximations which converge to the Friedrichs and Krein 
extensions of A and, as a bonus, obtain Theorem 4.19 and Theorems 5, 6, 7, and 8. We 
will see that the convergence of the Pade approximants ( equivalent ly, continued fractions 
of Stieltjes) can be reinterpreted as strong resolvent convergence of certain finite matrix 
approximations to Ap and Ak- 

Before we start, we note that our continued fraction expressions are related to but 
distinct from those of Stieltjes. Our basic object (the approximations to Ap) are of the 
form: 

2 (5.1a) 



-z + bo 



«6 



at 

-z + bi 



-2 + 62 H 

with two terms (an affine polynomial) at each stage while Stieltjes' are of the form 



1 

ciw H 



(5.1b) 



1 

C2 + — 

C3W + 



with a single term in each step and alternate constant and linear terms. This gives 
him twice as many continued fractions, so his alternate between the two sets of rational 
approximations we get. We will say more about this later. 

Until Theorem 5.29 below, we will suppose {7n}5^o ^^ ^ ^^^ ^^ Stieltjes moments. 

Given the Jacobi matrix (1.16), we will consider two N x N approximations. The first 
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is obvious: 



A 



[N] 





B. 


SIMON 




fbo 


ao 




"] 


tto 


bi 










bN-2 


aN-2 


\ 




aN-2 


bN-i) 



(5.2a) 



The second differs by the value of the A^A^ coefficient: 



A 



[N] 
K 



(bo ao 



V 



\ 



bN-2 ttN-2 . 

aN-2 bN-1—C^N-l/ 



(5.2b) 



where aN-i is chosen so that Aj^ has a zero eigenvalue. That such an aN-i exists is 
the content of 

Lemma 5.1. There is a unique aN-i > so that (5.2b) has a zero eigenvalue. A"^ > 
and aN-i obeys 

{bN-l - aAr_i)PAr_i(0) + aAr_2PAr_2 (0) = 



and 



Moreover, 



Pn{0) 

OiN-l = —CiN-l 75 TTTT • 



{bN - <^n){<^n-i) - a^-1 = 0. 



(5.3) 
(5.4) 
(5.5) 



Proof. Any solution of A'^ 'u = must obey the eigenfunction equation (4.1) and so be 
a multiple of -Pj(O), that is, 

Uj=Pj{0), j = 0,l,...,iV-l. 



This obeys the condition at the lower right of the matrix if and only if (5.3) holds. 

If Pn-i{0) were 0, then Ap ' would have a zero eigenvalue, which is impossible 
since Ap is strictly positive definite (since the form Sn-i is strictly positive definite). 
Thus, P/v-i(0) 7^ and so (5.3) has a unique solution. Since Ap is positive definite, 
this solution must have aN-i > 0. 

The eigenfunction equation for P/v(0) then yields (5.4). (5.3) for A^ ^ A^ + 1 and (5.4) 
yield two formulas for p ^ ^ .L . Setting these to each other yields (5.5). D 

The first main result of this section, towards which we are heading, is 
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Theorem 5.2 (includes Theorem 4.19). Let 

f+(z)^{5o,{A^;^^-z)-'5o) (5.6) 

and 

f-{z)^{6o,{A^F^-zr'So). (5.7) 

Then for x G (0, oo), 

fNi-^) > fN+ii-^) converges to (5o, {Ak + x)~'^6o) = /+(-x) 

while 

/^(-a:) < f^+ii-x) converges to {Sq, {Ap + x)"^(5o) = f~{-x). 

The convergence holds for any z G C\[0, oo) and is uniform on compacts. Moreover, for 
any solution p of the moment problem and x > 0, 

f-{-x)< /^</+(-x). (5.8) 

J x + y 

Remarks. 1. Notice that (5.8) is (4.42). 

2. Let A^^\/3) be (5.2b) with a^-i replaced by — /3. Then A^^^ is quite a natural ap- 
proximation: It is that A'-^' {(3) which is minimal among the non-negative matrices — rea- 
sonable to capture Ak, the minimal non- negative self-adjoint extension. From this point 
of view, Ap seems less natural. One would instead want to consider lim^^oo A'-^'{(3). 
In fact, one can show this limit is "essentially" Ap in the sense that 

lim (5o, (A[^](/3) - z)-'do) = (5o, (4^-'^ - z)-'do). 

Thus, the A^^\/3) interpolate between A^' and A^p~^\ 

Proposition 5.3. Let A'^ be the operator on i"^ which is A'^ ©0. Then A'^ converges 
to Ak in strong resolvent sense. 

Remark. For self-adjoint operators {An}'^^i and Aoo, we say An converges to A^o in 
strong resolvent sense if and only if for all z G C\M, {A^ — z)~^ip -^ {A^o — z)~^f for all 
vectors </?. It is a basic theorem (see Reed-Simon [33]) that it suffices to prove this for a 
single z G C\M. The same proof shows that if all A^ and A^o are non-negative, one then 
also has convergence for z G (— oo, 0). Moreover, one has convergence uniformly for z in 
compact subsets of C\[0, oo). 

Proof. Suppose first that the Hamburger problem is indeterminate. Thus, 7r(0) = 
(Po(0),...) g£2- Let PW(0) be the vector (Po(0), . . . , Pn-i(O), 0, 0, . . .). Then A^^pW 
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= SO (A^' - 2)-ipW = -z-ipt'^l for any z E C\R. Since ||(A^^ - z)-^\\ < \lmz\-^ 
for any such z, and ||p["l — 7r(0)|| -^ 0, we conclude that 

(A^l - ^)-V(O) ^ -z-V(O) = (Ak - ^)-V(O) (5.9) 

since Ai^7r(0) = 0. 

Suppose If = {A — z)r] for r] G -D(A), that is, a finite sequence. Then [Ak — z)r] = 
if = {A^^' — z)rj for n large, and thus lim„^oo(^^ — z)~^ (p = rj = {Ak — z)~^(p. Thus, 
{A^ — z)~^ip ^ {A — z)~^ip for if E Ran(A — z). In the determinate case, that is all (/j's; 
while in the indeterminate case, we claim that Ran(A — z) + [7r(0)] is all of £^ so (5.9) 
completes the proof. 

That leaves the proof of our claim that Ran(A — z) + [7r(0)] is £'^. We first note that 
7r(0) ^ D{A) by Proposition 4.4. {Ak -z)-^{ Ran(A - z) + [n{0)]) = D{A) +_[7t{0)] since 
{Ak - z)-^)tt{0) = -z-^Tv{0). But D{Ak)/D{A) has diml and 7r(0) ^ D{A) so {Ak - 
z)~^ {Ran{A - z) + [n{0)]) = D{Ak) and thus, since Im^ 7^ 0, Ran(A - z) + [n{0)] = £2, 
as claimed. D 

Proposition 5.4. Let A^j^ ^ he the {N + 1) x (A^ + 1) matrix which is A^j^ ' © 0, that is, 
it has zeros in its last row and column. Then 

ilf]<ASf+^l. (5.10) 

In particular, for x > 0, 

(5o, (Aif +^] + x)-'do) < (5o, (iSf ^ + x)-'do) = (5o, {A^^^ + x)-'do). 



Remark. Propositions 5.3 and 5.4 prove the part of Theorem 5.2 involving monotonicity 



of f^ and its convergence. 



Proof. B = Aj^ — Aj^ is an (A^ + 1) x (A^ + 1) matrix which has all zeros except for 
a 2 X 2 block in its lower right corner. This block is 

OAT-I «Ar-l 

aTv— 1 bN — ctjv 

which is symmetric with a positive trace and determinant zero (by (5.5)). Thus, the 
block is a positive rank one operator so B is positive, proving (5.10). D 

The monotonicity of the Ap s and their convergence is a special case of monotone 
convergence theorems for forms. Consider non-negative closed quadratic forms whose 
domain may not be dense. Each such form qs with D{qB) is associated to a non- 
negative self-adjoint operator B on D{qB). We define {B — x)~^ to be on D{qB)^ ■ 
This essentially sets S = 00 on D{qB)^ consistent with the intention that D{qB) is the 
set of (p for which qB{'p) < 00. 
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Given two forms qs and qc, we say qs < qc if and only if D{qB) ^ D{qc) (intuitively 
qc{'~p) < oo implies qsif) < oo) and for all Lp G qc, qci.^) > Qb{'^)- It is a fundamental 
result [12, 37] that if qB < qc, then for all x > and all </?, 

{if, {C + x)- V) < {^, {B + a:)- V)- (5.11) 

One monotone convergence theorem for forms [12, 37] says the following: Let {qAr^}'^=i 
be a sequence of quadratic forms with qA^ > ?a„+i- Then {A^ + x)~^ converges to 
the resolvent, {Aoo + x)~^, of an operator A^o- It is the operator associated to the 
form qA^ defined to be the closure of qA^ where D{qA^) = UD(qA„) with g^^((/?) — 
lim^i^oo qA„ i'^) provided that the form cja^ is closable. 

We can apply this to the Ap . Let qp be defined to be the form with D{qp ') = 
{{ifo,...,ipN-i,0,...,0,...)} and q^p\ip) = {ip,Aip). Then D{q^^^) C D{q^^+^^) and 
Q?\^) = qf^'Hv^) if V^ e D{qf^). Thus, q^^^ > qf+'^ (in essence, i^fl = A^^ © 
oo on C^"*"^). The monotone convergence theorem applies. A^o is just the closure of 
ip 1-^ {(f,A(f) on finite sequences, that is, the Friedrichs extension Ap, so we have the 
convergence and monotonicity part of Theorem 5.2 for /^ (using (5.11)): 

Proposition 5.5. Let x > 0. Then {6o,{A^p^^^ +x)-^6o) > (5o, (^jf ^ +x)-^6o) and 
this converges to (5o, {Ap + x)~^dQ) as N ^' oo. 

To head towards the remaining part of Theorem 5.2, viz. (5.8), we need expressions 

for ((/?, {Ap — z)~^ip) and (</?, (A'^ — z)~^Lp), which are of independent interest. 

Proposition 5.6. 

(«„,(4'"-^)-*) = -^ (B.12) 

Proof. The only possible eigenfunctions of Ap are the vectors n'-^^z) = {Po(z), . . . , 
Pn-i{z)). By the eigenfunction equation (4.1), this is an eigenfunction of Ap if and 
only if Pn{z) = 0. Thus, the A^ poles of (5o, {Ap — z)~^5q) are precisely the A^-zeros 
of Pn{z). The zeros are distinct because the eigenvalues are simple (since n^^\z) is the 
only potential eigenfunction, since eigenfunctions must have -uq 7^ 0). 

Now let the Bp be the (A^ — 1) x (N — 1) matrix obtained by removing the top row 

and left column of Ap . By Cramer's rule, (5o, {Ap — z)~^So) = [^] — —, so the 

N — 1 zeros of (5o, {Ap — z)~^5o) are precisely the A^ — 1 eigenvalues of Bp . The only 

possible eigenfunctions of Bp are ^{z)^^~^'^ = {Qi{z), . . . , Qn-i{z)) (since Qq{z) = 0, 

Qi{z) = — ). Thus, the eigenvalues of Bp are precisely ^'s with Qn{z) = 0. It follows 

that 

{do, {A'p - z) bo) = , 

Pn{z) 

and we need only show that cLn = — 1. 
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Since {6Q,{Af^^^-z)-^6o) = (bo-z)-^ a.nd Qi{z) = ^, Pi{z) = ^^^, we see that 
di = —1. On the other hand, Proposition 4.1 imphes that 

Qk{z) Qk-i{z) _ 1 _ fl 



Pk{z) Pk-i{z) ak-iPk{z)Pk-i{z) \z 
Moreover, (5o, (A — z)~^)do) = —z~^ + 0(z~^), so dk = dk-i, that is, d^ — —1 for aU 

Remarks. 1. The proof shows that for suitable cat (indeed, an induction shows that 

Cat = (— l)^(ao . . .aN-i)~^), we have 

Pn{z) = CN det{A^p^ - z) (5.13) 

and 

Q^{z) = -CN det(sSf ^ - z). (5.14) 

2. Since Ap is strictly positive definite, (5.13) shows once more that P/v(0) 7^ in 

the Stieltjes case. Since Bp is also strictly positive definite (as a submatrix of Ap '), 
(5.14) shows that Qn{0) 7^ also in the Stieltjes case. 

(5.13)-(5.14) imply some facts about the zeros of P/v and Qn: 

Proposition 5.7. AU the zeros of each Pn{z) and each Qn{z) are real. Moreover, there 
is exactly one zero ofQ^^z) and one zero of Pn-i{z) between any pair of successive zeros 

ofPNiz). 

Proof. The first assertion follows from the fact that the eigenvalues of a real-symmetric 
matrix are all real. The second assertion follows from the fact that if X is an A^ x A^ 
real-symmetric matrix and yisaniV — IxA^ — 1 submatrix, then by the min-max 
principle, there is exactly one eigenvalue of Y between any pair of eigenvalues of X. D 

Remark. Since the Qat's are orthogonal polynomials for another moment problem (see 
Proposition 5.16), between any two successive zeros of Qn^z), there is a single zero of 

Qn-i{z). 
Now define 

Mn{z) = Pn{z) - Pn{0) ^^^^ (5.15) 

Nm{z) = Qn{z) - Pn{0) ^^^ (5.16) 

since Ap is strictly positive, is never an eigenvalue and -Pj(O) 7^ for all j, so M is 
well-defined. Notice Mn{0) = 0. 
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Proposition 5.8. 

(*.(-*l^"' --'-'*"' -Ifi <-^' 

Proof. Let cn be the constant in (5.13) and let Bj^ be A^^ with the top row and 
left-most column removed. We will prove that 

Mn{z) = cn det(A^' - z); Nn{z) = -cn det{B^^^ - z), (5.18) 

from which (5.17) holds by Cramer's rule. 
Clearly, 

det(Ajf 1 -z)= det(yllf J - z) - a^v-i det{A^^~^^ - z) 



so 



CAT det(AiJ'^ -z) = Pn{z) - I3nPn-i{z), 



where 

Pn — «jv-i ■ 

CN-1 

But det(A'^ ^) = by construction. So P/v(0) — /3nPn-i{0) = and thus, cn det(A'f^ ' — 
z) = Mn{z). 

In the same way, 

CN det(s{^] -z) = -Qn{z) + PnQn-i{z) = -Nn{z). D 

Because of our convergence theorem, we have 
Corollary 5.9. Let {Xf^}fLi be the zeros of Pn{z) and let uf^ = -^f^^. Then 

'I2i=i^i "(-^"-^i ) converges to dupi)^), the Friedrichs solution of the moment problem. 
A similar formula holds for the Krein solution, dfixi^), with P,Q replaced by M,N . 

The following concludes the proof of Theorem 5.2. It extends the calculation behind 

(4.4): 

Theorem 5.10. Let dp solve the Stieltjes moment problem. Then for x > and N >1, 

QnJ-x) ^ I dp{y) ^ Nn{-x) 



Pn{—x) J x + y Mn{—x) 



PN{y)-PN{z) 



Proof. ^^ _ ^^ ' is a polynomial in y of degree A^ — 1 so 



dp{y)PN{y) 



PN{y)-PN{z) 

y- z 



(5.20) 
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On the other hand, (4.4) says that 



dpyy)- — - = Pn[z) / - — --^Qn{z). 



y- z 



y- z 



(5.21) 



Thus for z = —X with a; > 0, 

y+x J y+x J x+y 



+ Qn{-x)Pn{-x). 



Dividing by Pn{—x)'^, we obtain the left-most inequahty in (5.19). 

Similarly, since MAr(O) = 0, — ^^^ is a polynomial of degree A^ — 1 and so [y~^MN{y) — 
z~^MN{z)]/y — 2 is a polynomial of degree N — 2 iny, which is orthogonal to PN{y) and 
P/v_i(y) and so to M^i^y). Thus 



/ 



dp{y) 



y- z 



MN{y) Mn{z) 



0. 



Since for A^ > 0, 



Qn{z) = I dp{y) 
for each z, we see that for A^ > 1, 



PN{y)-PN{z) 

y- z 



(5.22) 



(5.23) 



Nn{z) = J dp{y) 

Therefore, for x > 0, 

dp{y) MN{y? Mn{-x 



MN{y)-MN{z) 
y- z 



(5.24) 



0< 



y + x 



y 



dp{y) 



MN{y) 



by (5.22)) 



(— x) J ' "" ' y -\- X 
Mn{-x)^ f dp{y) Mn{-x)Nn{-x) 



i-x) 



x + y 



+ 



i-x) 



Dividing by 



MAf(-x)- 



we obtain the second half of (5.19). D 



Next, we turn to Theorem 6 and the connection of Pade approximants to these finite 
matrix approximations. Given a formal power series Xl^o'^^'^"' ^^ define the Pade 
approximants f^^'^'{z) as follows (see Baker and Graves-Morris [4] for background on 
Pade approximants): We seek a function f'-^'-^^z) so that 



f[N,M]^,^ 



Bi^^^]{z] 



(5.25) 
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where A'-^'-^^ (z) is a polynomial of degree at most A^, B^^'^' (z) is a polynomial of degree 
at most M, as ^ ^ 0, 

N+M 

i=o 

and 

S[^'^1(0) = 1. (5.27) 

There is at most one solution, /, for these equations since if A, B are another pair of 
such polynomials, by (5.26), AB — AB = 0{z-^~^^~^^) must be zero since AB — AB is a 
polynomial of degree A^ + M. Thus, 

A[^'^l(^)S[^'^l(z) - i[^'^](^)S[^'^](^) = 0. (5.28) 

This implies A/B = A/B, showing / is uniquely determined as an analytic function. If 

deg A = N, degB = M, A, B relatively prime, (5.29) 

then (5.28) shows A — A and B — B, so A and B are uniquely determined. It can be 
shown ([4]) that if 

det{{KN-M+i+j-i)i<i,j<M) 7^ (5.30) 

then A,B exist and obey (5.29). 

There are degenerate cases where a Pade approximant exists, but A, B are not unique 
(e.g., for Xl^o ^n-2" = 1 + z, f^'^'^\z) = 1 + z can be written as A{z) = (1 + z){l + az), 
B{z) = (1 + az) for any a). In any event, if a solution of (5.25)-(5.27) exists, we say the 
[N,M] Pade approximant exists and denote it by f^^'-'^^z). 

In the context of Theorem 6, we are interested in the Pade approximants for the Taylor 
series of 

where p is a measure with finite moments. Without loss, normalize p so J dp = 1. Recall 
the context of Theorem 6. We have 

/•CO 

Jo 
and want to formally sum near z = 0: 



oo 

_z' 

n=0 
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If we define 

/•CX) 

7n = / x^ dp{x) 
Jo 

then near w = oo. 



oo 



G»= r^^--y^l^w—\ (5.32) 

Jo x-w ^^ 

Thus formaUy, 

Az) = \g,[-^-^. (5.33) 

We begin by noting: 
Proposition 5.11. As \w\ -^ oo, 

-^-E^.«»--'+o(«»-^") (5-35) 

^ ' j=o 

Proof. Let A be the Jacobi matrix for the moment problem {'jn}'^=o- Then 

'yn = {So,A^So). (5.36) 

On the other hand, we have an expansion converging near infinity for 

oo 

{do, (4^^ - w)-'6o) = - Y,(6o, {A^?^ydo)w-^-\ 

so (5.34) foUows from 

(5o, {Af^y5o) = (5o, A^o), J = 0, 1, . . . , 2iV - 1. (5.37) 

To prove (5.37), note that 

{A^^^y5o = A^5o (5.38) 

if J = 0, . . . , A^ — 1 and for B symmetric, 

{do,B''^o) = \\B'So\\ 

and that 

(5o, B^^+^6o) = {B^So, BiB^o))- (5.39) 
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To obtain (5.35), note that if Ap is replaced by A}^ , (5.38) still holds for j = 
0, . . . , A^ — 1, but (5.39) fails for j = A^ — 1 and so the analog of (5.37) only holds for 
i = 0,...,2iV-2. D 

Remark. While (5.37) holds for j = 0, 1, 2, . . . , 2N — 1, it never holds for j = 2N. For if 
it did, we would have for any polynomial, R, of degree 2N that 

{5o.RiA^^%) = {5o.R{A)5o). 
But this is false for R{x) = x^P/v_i(x)^ since 

{5o,R{A)5o) = \\A5n-i\\^ = a%_^ + a%_2 + b%_i, 

while 

{5o.R{A^^%) = \\A^PSN-if = a%_2 + b%_,. 

Proof of Theorem 6. By (5.33) and (5.34) as ^ ^ 0, 

-^fy4= Y.-^-'+0{z-^). (5.40) 

Now z^~^Qn{—^) is a polynomial of degree A^ — 1 since QAr(O) 7^ and z^ Pn{ — ^) is 
a polynomial of degree N since PAr(O) 7^ 0. Moreover, lim^^o ^^ Pn{—j) 7^ since Pn 
has degree N. Thus (5.40) identifies /[^"I'^l and -z^-^Qn{-^)/z^ Pn{-1)- Sim- 
ilarly, noting that since Mn{0) = (but M^(0) 7^ since A^ has simple eigenval- 
ues), z^ Mn{—-) is a polynomial of degree A^ — 1, and we identify f\.^~^'^~^^(^z) and 



-z 



N- 



-^Nn{--)/z^Mn{-^). Theorem 5.2 thus implies Theorem 6. D 



Remarks. 1. Since Pn and Qn are relatively prime (by Proposition 5.7) and similarly 
for M^v, Nn, we are in the situation where the numerator and denominator are uniquely 
determined up to a constant. 

2. Suppose 7n < Ca"' so p is unique and supported on some interval [0, R]. Then if R 
is chosen as small as possible, Xl^o '^nZ'^ has radius of convergence R~^, but the Pade 
approximants converge in the entire region C\(— 00, — i?~^]! 

Before leaving our discussion of Pade approximants for series of Stieltjes, we will 
consider the sequences fl^+^'^'ii^z) for j fixed. We will show they each converge as 
N ^ 00 with limits that are all distinct (as j varies) if the Stieltjes moment problem 
is indeterminate and all the same if a certain set of Stieltjes moment problems are all 
determinate. 

Given a set of moments {7j}°^i and £ = 1,2,..., define for j = 0, 1, . . . 

^f = ^^. (5.41) 
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Proposition 5.12. Suppose that {7j}°^o ^•^ ^ •^^^ ^f Stieltjes {resp. Hamburger) mo- 
ments. Then {7- }^o ^-^ ^^^ ^^^ ^f fnoments for a Stieltjes {resp. Hamburger) problem 
for i. = 1,2,... {resp. i = 2,4,6,...). This problem is indeterminate if the original 
problem is. 

Proof. We will consider the Stieltjes case. Let p E M.^{'y). Then 7^ x^ dp{x) = dp^^\x) 
solves the moment problem for 7*^^^. If pi,p2 G -^"^(7)7 then pi\ p2 are both in 
M.^{'~^^^^). li p\ = P2 7 then pi — p2 is supported at zero, which means it is zero 
since J dpi{x) = 70 = J dp2{x). Thus, if the 7 problem is indeterminate, so is the 7*^^^ 
problem. D 

Example. As Corollary 4.21 shows, there can be determinate moment problems so that 
7*-^-* is indeterminate (in the language of that corollary, 7*^^-* = 7). Thus, the converse of 
the last statement in Proposition 5.12 fails. 

As an aside, we present a proof of a criterion for indeterminacy of Hamburger [8] , espe- 
cially interesting because of a way of rewriting it in terms of a single ratio of determinants 
(see Theorem A. 7). Note this aside deals with the Hamburger problem. 

Proposition 5.13. A necessary and sufficient condition for a set of Hamburger mo- 
ments to be indeterminate is that 

(i) Er=ol^.(0)P<ooaW 
(ii) Er=ol^^'^(0)P<oo 
where P\ {x) are the orthogonal polynomials for the 7*^^-' moment probl 



,(2)- 



Proof. If {7j}°^o ^s indeterminate by Proposition 5.12, so is {'Jj^'JJ^q, and then (i), (ii) 
follow by Theorem 3. The converse is more involved. The intuition is the following: If dp 
is an even measure, then P2J {x) = y/^ ^'^"^ , so Pm {0) = y/^ — ^~"(0) (both sides 
are zero if m is odd). So (ii) is equivalent to E^o I Sx \'^ < ^^ ^^'^ condition (v) of 
Theorem 3 says that (i), (ii) imply indeterminacy. Our goal is to prove in general that 
when (i) holds, P- (0) and q^^ (0) are essentially equivalent. 
Let Sr,{x) = pP{x), r]n = ^(0), and define 



ctn = xSn{x) dp{x) (5-42) 

(3n= [[xSn{x)][Pn+l{x)]dp{x), (5.43) 



where p is any solution of the 7 moment problem. Since Sn is an orthogonal polynomial 
for dp^ \ we have 

xSn{x)]x-' xdp{x) = 0, J = 0, 1, . . . , n — 1 
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and thus 

/ xSn{x)Pj{x)dp{x) = a^PjiO), j = 0, l,...,n 

so the orthogonal expansion for xSn{x) is 

n 
xSr^ix) = /3r,Pn+l{x) +anJ2PjiO)Pji^)- (5-44) 

Since xSn{x) vanishes at zero, 

n 

Since 



we have by (5.44) that 



^2=Pl + alJ2\PMf- (5-46) 

i=o 



Moreover, taking derivatives of (5.44) at a; = 0, 

n 

Sr,{0)=Pr,r]r,+i+anJ2Pj(^H- (5.47) 

By (5.46), (3n is bounded, so by (5.45), 

|«n|<Ci|Pn+i(0)|. (5.48) 

By hypothesis, YlTLo l-Pj(0)P < C)0, so Pn+i{0) -^ and thus by (5.46), (3n -^ -Jyi as 
n ^ oo. Since /3^ > for aU n (this follows from the fact that both xS^ix) and Pn^\{x) 
have positive leading coefficient multiplying a;"^"'"^), we see that /3~^ is bounded. Using 
(5.48), the Schwarz inequality on the sum in (5.47), and X^^o l-^j(0)P < oo again, we 
see that 



Iryn+il' < c{\S^m^ + |Pn+i(0)P Y. I^^l') 

<C(|5.(0)|2 + |P.+i(0)|2) l + 5^|r7,| 

^ j=0 
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It follows that 



n+l 



1 + Ei^^-i 



and so by induction that 

n+l 



<[l + C(|5.(0)|2 + |P,+i(0)|2)] 



1 + Ei^^-i 



n+l 



sup J]|r7,f < sup ll[l + C{\S,.^{0)\' + \P,m^)]<oo 



l<n<oo 



i=o 



l<n<oo • 1 



since YlTLo \Pj(^)\'^~^Y17Lq I'S'j(0)P < cx) by hypothesis. Thus, if (i) and (ii) hold, r]n G £2 
and thus, the problem is indeterminate by Theorem 3. D 

Theorem 5.14. Let {7n}5^o ^^ ^ ^^^ ^f Stieltjes moments. Fix i > 1 and let Pj^ (z), 
Q]^ (z), Mj^ (z), and N^'{z) be the orthogonal polynomials and other associated poly- 
nomials for the 7*^^^ moment problem. Let f^^'^\z) be the Fade approximants for the 
series of Stieltjes J^TLo'^j^'' where Kj = { — lY'^j. Then 



£-1 









and 



J=0 









(5.49) 



(5.50) 



In particular: 

(1) For each i, {—l)^f^^~^^~^'-^\x) is monotone increasing to a finite limit for all 
X G (0, 00). 

(2) For all z E C\(-oo,0], lim^v^oo fl^+^-^'^]{z) = fi{z) exists. 

(3) The 'j^^' moment problem is (Stieltjes) determinate if and only if fi{z) = /^_|_i(z). 

(4) // the 7*^^^ moment problem is determinate, then fo{z) = fi{z) = ■ ■ ■ = fe+i{z). 

(5) // the 7*^^^ moment problem is indeterminate, then as x -^ 00, 



W, 



f,+^{x) = (-l)^^,^^^({0})x^ + O(x^-^) 



(5.51) 



where /j)^ is the Krein solution of the 7*^^'' problem. In particular if i > 1, 
—fg+i{x) is not a Herglotz function {as it is if the problem is determinate). 

Remarks. 1. Thus, it can happen (e.g., if 7 is determinate but 7*^^^ is not (cf. Corol- 
lary 4.21) that /[^~i'^] and /t^'-^l have the same limit, but that /t-^+i'^l has a different 
limit. 
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2. If the {7j}^o obey a condition that imphes uniqueness and which is invariant 
under 7^ — > 'jj+i (e.g., the condition of (1.12b) of Proposition 1.5), then aU 7*^^^ are 
determinate; see Theorem 5.19 below. 

3. Even for i = 1 and consideration of /[-^~i'-^l and fl^'^\ we have something new 
— the remarkable fact that xdiiK{x)/'^i^ is d^p {x). Multiplication by x kills the point 
measure at a; = and produces a measure supported on some set [R, 00) with R > Q. 
We also see that monotonicity for /^ implies monotonicity of /^. The direction of 
monotonicity flips because of the minus sign ((—1)^ = —1) in (5.49). 

4. In particular, we have Pj^'{x) = cnx~^Mn^i{x), that is, for any i j^ m, 
J x~^Mm{x)M£{x) dp{x) = 0, something that can be checked directly. 

5. This theorem implies that all the /I-^'^l with N > M — 1 exist and their denomi- 
nators obey a three-term recursion relation. 

6. The connection in Remark 3 extended to xdpj^ (x) = cid^p (x) implies that if 
any 7*^^^ is indeterminate, then duxix) is a point measure. Thus if k^ = (—1)"^ J^ x"^ dp{x) 
where p is associated to a determinate Stieltjes problem and dp is not a point measure, 
then all the /^(2)'s, £ = 0, 1, 2, ... , are equal. 

Proof. It is easy to check that the right side of (5.49) and (5.50) are ratios of polynomials 
whose degrees are at most of the right size. (Because P^ (0) 7^ 7^ M^^ (0), it is 
easy that the denominators are always precisely of the right size.) Since Mj^ (z) and 
Pp^ (z) are of degree n, the rationalized denominators do not vanish at - = 0. By 
Proposition 5.11 for the "f^^' problem, they have the proper asymptotic series to be 
fW+i-^^m(^z) and /[^+^-i'^-i](2), respectively. That proves (5.49) and (5.50). 

Assertions (1), (2), (3) are then just Theorem 6 for the "f^^' moments. To prove asser- 
tion (4), note that if 7*-^-* is determined by Proposition 5.12, so are 7*-^~^-*, X^~'^\ • • • , 7*-^ , 
7, and so by (3), we have Z^+i = ft = ■ ■ ■ = fo- 

To prove assertion (5), note that if the 7*^^-* moment problem is indeterminate, dp\^ {{0}) 
> 0. Thus, 

lim /i?^ = ^«{0}>0. (5.52) 

' ooj 1 + xy 



X 



By (5.50), 






til r^ A.M) 



i)^7.x^- + (-1) v^ r ^f^ 

Jo l + xy 



so (5.52) implies (5.51). D 

To deal with the sequences fl^-'^+^'^l (^z) with £ < 0, we need to introduce yet another 
modified moment problem associated to a set of Hamburger (or Stieltjes) moments. To 
motivate what we are looking for, let p G ^^^(7) and let Gp{z) be the associated Stieltjes 
transform. Then —Gp{z)~^ is a Herglotz function which has an asymptotic series to all 
orders as ^ ^ zoo. Since Gp{z) ~ —^(1 + 71-2""*^ + 12Z~^ + 0{z~^)), 

-Gp{z)-^ ~ ^ - 71 - (72 - 7i)^"' + 0{z-^) 
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SO (by the proof of Proposition 4.13) the Herglotz representation of —Gp{z)~^ has the 
form 

~G,{zr' = . - 71 + (72 - 7?) / ^ , (5.53) 

J X z 

where dp is a normahzed measure. Its moments, which we will call 7- , only depend 
on the asymptotic series for Gp{z), and so only on the original moments. We could do 
everything abstractly using (5.53), but we will be able to explicitly describe the relation 
between the moment problems. 

The key is the following formula, a Ricatti-type equation well known to practitioners 
of the inverse problem [6] (which we will, in essence, prove below): 

—mo{z)~^=z — bo + aQmi{z)j (5.54) 

where mo{z) = {do,{A — z)~^do) is just Gp{z) for a spectral measure, and mi{z) = 
{5oj (AI^I — z)~^5o), where A'-^^ is obtained from A by removing the top row and left-most 
column. ( 5.54) is just (5.53) if we note that 60 = 7i and ag = 72 — 7i . Thus, we are led 
to define 7- as the moments associated to the Jacobi matrix A'-^' obtained by removing 
the top row and left-most column, that is. 
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Proposition 5.15. Let {7j}°^o ^e a set of Hamburger moments. Then the {7J Jj^o 
Hamburger problem is determinate if and only «/{7j}j^o ^•*- 

Proof. Let A^^' be A^^l with a row of zeros added at the top and column of zeros on the 
left. Since A^^^ = 0©^!^], A^^^ is essentially self-adjoint if and only if At^l is. A-A^^^ is a 
matrix with three non-zero elements and essential self-adjointness is provided by bounded 
perturbations. Thus, A is essentially self-adjoint if and only if A^^l is. By Theorem 2, we 
have the equivalence of determinacy for the Hamburger problem. D 

Remark. As we will see shortly, this result is not true in the Stieltjes case. 

Let Pp^ (x), Qjy (x), and fp^{z) = (077-7 be the polynomials and finite matrix 

approximation for the 7''^-' problem. Then 

Proposition 5.16. 

(i) Pi^\x)=aoQN+i{x) (5.56) 

(ii) q!J)(x) = -- Pn+i{x) + ^— ^ Qn+i{x) (5.57) 

(iii) - f]^+,{zr' -z-bo + alff{z) (5.58) 
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Remarks. 1. (i) implies that the Qa^'s are orthogonal polynomials for some measure. 

2. (5.58) is (5.54). 

Proof. For each fixed x, Pjy {x),Qj^ (x) obey the same difference equation as Pn+i{x), 
Qn+i{x) so we need (5.56), (5.57) at the initial points A^ = —1,0 where it is required 
that 

Go 

{Q_i{x) = —-^ since we don't make the convention a_{ = 1, but rather a_{ = gq 
consistent with On = an+i)- This is consistent with (5.56) if we note that 

Qo{x)=0, Qi{x) = -; Po(x) = l, Pi{x) - ^"^ ~ ^'^ 



This proves (i), (ii). 
They in turn imply 



Go ao 



2Qn\^) , PN+lix) 



X — bo — 



which, by (5.12), proves (iii). D 

Proposition 5.17. Let {7j}°^o ^^ ^ ^^^ of Stieltjes moments. If the Hamburger problem 

is indeterminate, then the 7*^*^^ Stieltjes problem is indeterminate {even if the 7 Stieltjes 
problem is determinate). 

Proof. As we will see below (Proposition 5.22 = Theorem 8), a set of Stieltjes mo- 
ments which is Hamburger indeterminate is Stieltjes determinate if and only if L = 
lim^^^[-%M] is infinite. By (5.58), 



Pn (0) ■ 



alL(-^'^=bo-L-\ 



Since L > 0, L^^' is never infinite. D 

Remarks. 1. Thus, if 7 is Hamburger indeterminate but Stieltjes determinate, 7 and 7*^'^'' 
have opposite Stieltjes determinacy. 

2. The spectral way to understand this result is to note that if 7 is Hamburger 
indeterminate, the Friedrichs extensions of A and A'-^' have interlacing eigenvalues. Given 
that A > 0, A[^] must be strictly positive. 

For £ < and integral, we let 

(0) 

7f^ = [7^°^];-'^ - ^'-' 



J L' Jj (0) ■ 



Then: 
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Theorem 5.18. Let {7n}5^o ^^ ^ ^^^ ^f Stieltjes moments. Fix i < and let Pf^ {z), 
Qj^ (z), Mj^ (z), and Nj^' {z) he the orthogonal polynomials and other associated polyno- 
mials for the "^'^^'^ moment problem. Let f^^'^\z) he the Fade aproximants for the series 
of Stieltjes Y^^qKjZ^ where Kj = (— l)-'7j. Then 



and 

r[N+£-l,N+l]^^^ 



-£-1 









VzF 



N+e- 



-i(-i)J 



(5.59) 



-£-1 



\£+iA0)-i+2 



1-71^- E i-iyifz'+'H-iY+W-iz- 

J=0 



N 



N+l- 



-i(-i) 



^<l.(-i) J 



(5.60) 
In particular, 

(1) For each i, {—l)^f^^~^^~^'-'^\x) is monotone increasing to a finite limit for all 
X G (0, oo). 

(2) For all z G C\(-oo,0], liniTv^oo fl^+^-^^^]{z) = fi{z) exists. 

(3) The 7*^^^ moment prohlem is Stieltjes determinate if and only if fg{z) = fe-i{z). 

(4) // the 7*^^^ moment prohlem is determinate, then fo{z) = f-i{z) = ■ ■ ■ = fi{z) = 
f£-l{z). 

(5) // the 7*^^^ moment prohlem is indeterminate, then as x ^ oo, 

/,_i(x) = (-l)^+i(7i°])-V-VSJ^({0})-i+O(x^-3), 



where HjJ is the Krein solution of the 7^^^ prohlem. In particular, —fi-ii 
not a Herglotz function {as it is if the problem is determinate). 

Proof. By (5.58), we have a relation between formal power series: 



X] IS 



i=o 



from which one obtains 



2\2 



= 1+712;- (72 -7i)^ 



'-j=o 



-ly^fz^ 



(5.61) 



flM,N+2]^^^ = (l+^i^-(^2-7lV/^°^f^'''k^))-'- 

Thus, (5.59) is just (5.49) and (5.60) is (5.50). The consequence (l)-(5) follows as in the 
proof of Theorem 5.14. D 



We summarize and extend in the following: 
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Theorem 5.19. Let YlT^o '^j-^"' ^^ '^ series of Stieltjes. Then for each £ = 0, ±1, ±2, . . . , 

fi{z) = liuiN^oo f^^^^''^\z) exists uniformly for z in compact subsets o/C\(— oo, 0]. 
Moreover, 

(1) // any two fg 's are different, then all of them are meromorphic. 

(2) // the \k,j\ = 7j moment problem is determinate and the measure solving the 
moment problem is not a discrete point measure, then all fi 's are equal. 

(3) // \kj\ < C^{2j)l, then all fi's are equal. 

Proof. The first assertion is Theorem 5.14 for i > 1 and Theorem 5.18 for £ < 0. If 
some fi ^ fi+i, then some 'j^^' are indeterminate, so the corresponding dpp^s are pure 
point and /^'s meromorphic. This proves (1). Under the hypothesis of (2), /o is not 
meromorphic, so by (1), all /'s are equal. To prove (3), an induction using (5.61) proves 

\lj I ^ C'"'(2j)!, so by Proposition 1.5, all the "f^^' moment problems are determinate. D 

This completes our discusion of Pade approximants for series of Stieltjes. We return 
to consequences of Theorem 5.2 for the study of the Stieltjes moment problem. 

As we have seen for any x > 0, — p ;~% is monotone increasing and — now— r- 

Qiv(o) 

Pjv(O) 
M^(0) 

Proposition 5.20. 



is monotone decreasing. By taking limits (recall Mn{0) = 0), we see that — pfjoj is 
monotone increasing and iJ^ (J) is decreasing. 



(i) lim -— = V ^ dupiy) 

where fxp (resp. hk) is the Friedrichs {resp. Krein) solution. 
Proof. By Theorem 5.2 and Proposition 5.6, for x > 0, 

Pn[-x) J 
so taking x to zero and A*" to infinity, we see that 

lim -^421 < fy-u^Ay)- 



On the other hand, since — p (_^) = (5o, {Ap + x)~^6o) is monotone increasing as x 
decreases, we have for each A^ and x > 0, 

Qn{-x) ^ ^.^ Qn{0) 



Pn{-x) ~ N^oo Pn{0) 
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SO taking A^ to infinity for fixed a; > and using Theorem 5.2, 

Qn{0) 



(x + y) ^diiriy) < lim 



Pn{0) 
Taking x to zero, we see that 

Qn{0) 



J N^oc 



Pn{0) 



so (i) is proven. 

The proof of (ii) is similar. By Theorem 5.2 and Proposition 5.8, 

xNis[{—x) ^ f X 



> [ ^—dfiK{y)>fiK{{0}). 
J y + x 



Mn{-x) 
so taking x to zero and N to infinity. 

On the other hand, since ~ ^m^{-x) ~ i^Qi^i^K + x)~^5q) is monotone decreasing as 
X decreases, we have for each A^ and x > 0, 

xNn{-x) ^ ^.^ Nn{Q) 



Mn{-X) - N^OO M'jy{0) 

so taking A^ to infinity for fixed a: > and using Theorem 5.2, 

d^Kiy) > lim ^ 



Taking x to zero, we see that 

^k({0}) > ^lim ^^^^ , 

so (ii) is proven. D 

This leads us to define (note that we use j^^L? , not jj^ /q-^ so M, L E (0, oo) U {oo}): 

L= hm -^4^ (5.62) 

M= lim ^^^, (5.63) 
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SO Proposition 5.20 says that 

By (4.25), ^^({0}) = l/Er=o I^^WP, so 

oo 

M = J2\Pn{0)\^. (5.64) 

n=0 

Note. (4.25) was only proven in the indeterminate case, but the argument apphes in 
the determinate case also. If x is an eigenvalue of a solution fx of the moment problem 
associated to a self-adjoint extension, then fx{{x}) = 1/Xl^o iPnix)]"^. 

Theorem 5.21 (= Theorem 7). Let {7n}5^o ^^ ^^^ moments of a Stieltjes problem,. 
Then the problem, is indeterminate if and only if 

L < oo and M < oo. 

Equivalently, the problem is determinate if and only if 

L = oo or M = oo. 

Proof. If L < oo and M < oo, then j y~^dnF{y) < oo, while hk{{0}) > 0, so clearly, 
fJ-F 7^ fJ-K and the problem is indeterminate. Conversely, if the problem is indeterminate, 
then by Proposition 3.1, a = inf spec(Ai;') > 0, so / y~^ dfipiv) < ct~^ < oo and L < oo. 
Moreover, since the Stieltjes problem is indeterminate, so is the Hamburger problem; and 
thus by (5.64), M< oo. D 

Theorem 5.22 (= Theorem 8). Let {'jn}'^=Q be a set of Stieltjes moments. Then 
the Stieltjes problem is determinate while the Hamburger problem is indeterminate if and 
only if 



CO 

n=0 



_|Qn(0)r<oo (5.65) 

n=0 

and L = oo. 



Proof. For a set of Stieltjes moments — pLl is positive and monotone increasing, so 

|Qn(0)| > '^'|P,\^]|°^' - Since Qi(0) ^ 0, we see that (5.65) implies also that M < oo. 
Thus, (5.65) is equivalent to the Hamburger problem being indeterminate. Given that 
M < oo. Theorem 5.21 says that determinacy of the Stieltjes problem is equivalent to 

L = oo. D 
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As our next topic, we will further examine the conditions M < oo and L < oo to see 
they are conditions of Stieltjes and Krein in a different form. Define 

and 

^^.^_^^, .>2 (5.67a) 



iV,(0) iV,-i(0) 

mm 

iVi(O) 



rm-^. (5.67b) 



By the monotonicity properties of this section, in > 0, rrin > 0. By definition and 

Qo(o) = 0, 

Qn{0) ^ 



Pn{0) 



n=l n=l 



-n 



?Br 



Proposition 5.23. For N >l, 



mN = \PN-m? (5.68) 

£^ = -[a^_iPAr(0)PAr_i(0)]-\ (5.69) 



\m 



Remark. Since A'p > 0, Pn{z) has no zeros on (— oo,0] and thus, Pn{z) has the same 
sign near — oo and 0. But Pn{z) = c^z^ + lower order with cn > 0, so {—1)^Pn{0) > 
(cf. (5.13)). Thus, Pis[{0)Pn-i{0) < and the minus sign in (5.69) is just what is needed 
to ensure that In > 0. 

Proof. A^j^' has {Po(0), . . . , -P/v-i(0)} as its eigenfunction with eigenvalue zero. Thus, 
lim,^oa:(5o, {A^J^^+x)-^6o) = |Po(0)|7 Ef="o' l^j(0)p. Since Po(0) = 1, we see that 

^^g|P,(0)P. (5.70) 

which implies (5.68). (5.69) follows immediately from Proposition 4.1. D 



^^^^^" Pn(0) ' 


/ 6 -^ V7 


^^ ^" Pn(0) ' 


n> 


^^(a:) = -a„_iM„(-x)P,_i(0), 


n> 1 


if,(a;)=a,_iiV,(-a;)P,_i(0), 


n> 1. 
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Corollary 5.24. If A given by (1.16) is the Jacobi matrix associated to an indeterminate 
Stieltjes moment problem, then 

oo 
n=0 

In particular, ifYl^=o an = oo for the Jacobi matrix associated to a Stieltjes problem, 
then the problem is determinate. 

Proof. If the problem is indeterminate, by Proposition 5.23 and Theorem 5.21, 
an ^^^ I Pn(0)Pn-i (0)1"^/^ and |Pn(0)Pn-i(0)|^/^ both he in £^ so their product lies in 

Now we will define, following Stieltjes and Krein, functions 

(5.72) 

(5.73) 
(5.74) 

We claim that 
Proposition 5.25. 

(i) t/n(0) = 1 

(ii) Gn(0) = 

(iii) i/n(0) = 1 

(iv) Un{x)H^{x) - Vn{x)Gn{x) = 1 forn > 1 
(v) f+{-x) - f-{-x) = ^„(,)G„(,) forn > 1 

Proof, (i) is immediate from the definition (5.71), (ii) from Mn(0) = 0, and (iii) follows 
from Proposition 4.1 and the definition (5.16) of N^ (and explains why we multiply M 
and N by an-i-Pn-i(O)). To prove (iv), we note that for some constants an,/9n,7n, 

Un{x) = anPn{-x), Vn{x) = ~anQn{-x) 

Gn{x) = (3riPni-x) + 'JnPn-li-x), Hn{x) = -(3nQn{-x) - 'JnQn-li-x), 

SO 

Un{x)Hn{x) - Vn{x)Gnix) = -an'yn[Pn{-x)Qn-l{-x) - Qn{-x)Pn-l{-x)] 

is constant as x is varied for fixed n. But at x = 0, by (i)-(iii) this combination is 1. (v) 
follows from (iv) and the definitions. D 

The following says that for x > 0, U,G can be associated with the equation of motion 
of a string of beads (see [1]). This is the starting point of deep work of Krein [17]. 
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Theorem 5.26. 

(i) Unix) -Un-i{x) = i^Gnix), n>l 
(ii) Gn+i{x) - Gn{x) = nin+ixUnix), n>l 
(iii) Gi{x) = niix Uo{x) = 1 

Proof, (i) By definition of Un, 

Unix) - Un-i{x) = ^p^~^^ = lnGn{x) 

by (5.69) and (5.73). 
(ii) We have for n > 1, 

anPn+l{-x) + bnPn{-x) + an-lPn-l{-x) = -xPn{-x). (5.75) 

(5.75) for X = implies that 

bn = -anPn+iiO)Pn{0)-^ - an-iPn-i{0)Pn{0)-\ (5.76) 

which we can substitute into (5.75) to obtain 

anMn+l{-x) - an-lPn-liO)Pn{0)-^Mn{-x) = -xPn{x). 

Muhiplying by — ^^,(0) and using the definitions (5.71)/(5.73), we see 

Gn+l{x) - Gnix) = xPniOfUnix) 

SO (5.68) implies the result. 

(iii) Uo{x) is a constant so Uo{0) = 1 implies Uo{x) = 1. Gi is a linear polynomial with 

Gi(0) = 0, so Gi{x) = G[{0)x. For any n, ^^(0) = lim.^o ^^7^ = E^^i^j- □ 

We can use these equations to get further insight into Theorem 5.2. First, we obtain 
explicit bounds on the rate of convergence oi f^{—x) — f^{—x) to zero in the determinate 
case, where either ^ rrij or ^ £j or both diverge. 

Theorem 5.27. For x > 0, 

f^{-x) - f^{-x) < rnix^(J2^j)[J2^A 



Proof. It follows from Theorem 5.26 and induction that Un{x) and Gn{x) are non- 
negative for X > and then that they are monotone in n. In particular, Un{x) > 1, 
Gnix) > niix from which, by (i), (ii). 



Un{x)>mi(j2^j) 
Gn{x) > i^mjjx. 

^7 = 1 ^ 



J = 

(v) of Proposition 5.25 then completes the proof. D 

And we obtain directly that in the indeterminate case, /^ are meromorphic functions. 
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Theorem 5.28. (i) For any z E C, we have that for n > 1, 

n n 

\Un{z)\<l[{l + ej)l[{l + m,\z\) 

n— 1 n 

iG.(^)i<n(i+^.)n(i+^^-i^i)- 

j=i i=i 

(ii) If M < oo and L < oo, then for all z E C, Un{z), Vn{z), Gn{z), Hn{z) converge to 
functions Uoo{z),Voo{z),Goo{z), Hoo{z), which are entire functions obeying 

\f{z)\<C,eMe{z)) 
for each e > 0. 

Remark. In terms of the Nevanlinna functions A, B, C, D, one can see (using the fact 
that the Friedrichs solution is associated to Bt with t = L) that Gooi—z) = —D{z), 
Hooi-z) = C{z), Uooi-z) = -B{z) - LD{z), and Vooi-z) = A{z) + LC{z), where 

Proof, (i) foUows by an elementary induction from Theorem 5.26. Similarly, one follows 
the proof of that theorem to show that Vn, Hn obey 

Vn{x) -Vn-l{x) = £nHn{x), n>l 

Hn+l{x) - Hn{x) = mn+lXVnix), n > 1 

Hi{x) = l 1/0(2^) = 
and obtains inductively that 

n n 

n— 1 n 

\Hr,{z)\<\[{l + l,)\[{l + mM)- 
i=i j=2 

Thus, if L, M < cxD, we first see that U, V, H, G are bounded and then, since Ylij>n ^j ~^ 
0, J2j>n''^j ~^ 0' ^^^^ each sequence is Cauchy as n — *> 00. The Cg bound is easy for 
products n^i(l + "^jkD with '^'^ nij < 00. (iii), (iv) are then immediate from the 
definitions. D 
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To make the link to Stieltjes' continued fractions, we need to note the relation between 
(In, bn, in, and nin- Immediately from Proposition 5.23, we have 

an = (C+iV™n+i"^n+2) , n > (5.77) 

and by (5.76), 

bn = m-l,iC+i-l,), n>l (5.78) 

6o=m^^£^^ (5.79) 

(Parenthetically, we note that these equations can be used inductively to define rrij , ij 
given ttj, bj. We have mi = 1, so (5.79) gives £i. Given £i, . . . ,£j and rrii, . . . , nij, we 
can use aj-i and (5.77) to find TOj_|_i and then bj and (5.78) to get ij+i-) 

Stieltjes' continued fractions are of the form (5.1b). Let 

do = C^\ dn = {CnCn+l)~^ (5.80) 

for n > 1 so (5.1b) becomes 

dn 



di 

w -\ 



1+ ^^ 



w + . 



Now use the identity 

, Pi ^n /3l/32 

^ ^ ^~ =W+ Pi- n , f^ ^ 

(^2 P2 + f{w) 

to see that (5.1b) has the form (5.1a) if w = —z, do = 1, and 

bo = di (5.81a) 

bn = d2n+l + d2n (5.81b) 

fln = d2n+ld2n+2- (5.81c) 

Thus, (5.80) and (5.81) are consistent with (5.77)-(5.79) if and only if 

C2J-1 = ruj, C2j = ij. (5.82) 



Thus, we have seen that (5.1a) is (5.1b) if w = —z and the c's are given by (5.82). 

■'3- 



Stieltjes' criterion that depends on whether or not Xl^i <^j < oo is equivalent to L < oo 



and M < oo. 
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The connection (5.82) is not surprising. In (5.1b) if ty = 0, the continued fraction 
formaUy reduces to C2 + C4 + ■ ■ ■ consistent with Xl^i ^3 — ~ I^a^^oo p(o) ■ ^^ the 
other hand, if we multiply by w, the continued fraction formally becomes 

1 



1 

ci + 



C2W H 

C3 H 

which, formally at w = 0, is (ci +cs + - ■ ■ )~^, consistent with J27Li ^j = liiiiAf^oo liiiixj.o 

L MNi-x) i- 

We conclude this section by discussing Pade approximants for the Hamburger analog 
of series of Stieltjes. A series of Hamburger is a formal power series Xlf^o '^j^'' with 



Kj = {—ly / X-' dp{x) 

J — 00 

for some measure p. We will begin with the "principal" Pade approximants p^~^'^\z): 
Theorem 5.29. Let YlTLo ^i^"' ^^ ^ series of Hamburger. Then: 

(i) The /[^~^'^l(z) Pade approximants always exist and are given by 



N-l 



QnH: 



/[^-'^](^) - (5o, (1 + .<J)-^5o) = - ^"7,;^ . (5.83) 

(ii) // the associated Hamburger moment problem is determinate, then for any z with 
Ivnz 7^ 0, 

hm /[^-i'^](^) = -f^!4^!4T^ (5-84) 

exists and equals J jrj^ for the unique solution, p, of the moment problem. 

(iii) The sequence f^^~^'^\z) is pre-compact in the family of functions analytic in 
C_|_ (in the topology of uniform convergence on compact sets). 

(iv) Any limit of f'-^~^'^^ (z) is of the form J j^§^ where p is a von Neumann solution 
of the Hamburger problem. 

Proof, (i) The proof that (5.12) holds (Proposition 5.6) is the same as in the Stieltjes 
case, and from there we get (5.34) as in the Stieltjes case. Since lim^^^o -2^-P/v(-) 7^ 0, 
we see f^^~^'-^\z) exists and is given by (5.83). 
(iii) Since 

|^;[^-i,A^](^)| = |(5o,(4^] +^-i)-i5o)| < \Im{z-')\-\ 



74 



B. SIMON 



we see that 



|/[^-i'^l(z)|<|^||Im(^)|-i 



is bounded on compacts. By the Weierstrass-Vitah theorem, such functions are pre- 
compact. 

(ii), (iv) By the proof of Proposition 4.13, if dp is any measure whose first M + 2 

moments are 70, 7i, . . . , 7m+2 5 then 



,M+1 



M 



G,{ty) + J2i- 



_^^n+ly-n-l^^ 



n=0 



< 



ay 



(5.85) 



[Nh 



with a only depending on {7} (either a = 7M+1 or a = 2 [7m+7m+2])- Since {Sq, {Ap ySo) 
= 7j for j < 2N — 1, we see that for M fixed, (5.85) holds uniformly for A^ large if 
Gn{z) = — -2~^/f'^~^'^l(— ^), so any limit point of the G's obeys (5.85) and is Herglotz. 
Thus by Proposition 4.13, it is of the form Gp{z) with p G Ai^ {•y). Thus the limit for 

y[7V-i,Ar]^^-j as A^ — s> oo, call it /(-z), has the form 



fiz) 



dp{x) 

1 + xz 



with p G 7^1^(7). In the determinate case, p must be the unique solution, so all limit 
points are equal. Thus compactness implies convergence, and we have proven (ii). 
In the indeterminate case, we note that 6j = Pj{Ap )So for j = 0, . . . , N — 1 so 



(-41?" 



N-l 



-'So = J2 (PMf%, (4^^ - z)-'So)Sj. 



dm 



As in the proof of Theorem 4.3, we conclude that (with Gn{z) = (5o, {Ap — z) ^5o)), 

N-l 



j=0 



Im Gn{z) 
Im^; 



Thus, since Q,P E £2 if Gn{z) -^ G{z) through a subsequence, the limit ( = G{z) 
obeys equality in (4.5). By Theorem 4.14, any such solution of the moment problem is a 
von Neumann solution. D 

We want to show that in many cases, lim/[-^~^'-^J(z) will not exist. Consider a set of 
Stieltjes moments, {'jn}^=o ^^^ the associated even Hamburger problem with moments 



^2n — In, ^2n+l — 0. 

Given p G M^i'j), let p G A^-'^(r) be given by 

dp{x) = ^{X[o,oo){x) dp{x'^) +X{-oo,o]{x) dp{x'^)) 



(5.86) 



(5.87) 



Gp{z) : 


= zGp{z^) 


P2n{x) ■■ 


= Pn{x') 


P2n-l{x) ■■ 


_ dr,Mn{x^) 
X 


Q2n{x) - 


= xQn{x'^) 


Q2n-l{x) ■■ 


= dnNr^{x^). 
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SO that Theorem 2.12 says that this sets up a one-one correspondence between even 
integrals in A4^ (T) and A4^{'y). Let Pn.Qn be the orthogonal polynomials for F and 
P„, Qm Mn, Nn for 7. Define dn by d'^ f Mn{x)'^x~^ dp{x) = 1. Then we claim that 

(5.88) 
(5.89) 

(5.90) 

(5.91) 
(5.92) 

To prove (5.89) and (5.90), we note that the right sides have the correct degree, and 
since P„(a;^) is even and — nlELd jg odd, they are dp orthogonal to each other. Easy 
calculations prove the dp orthogonality of two Pn(a^^)'s and two — "^^ ' 's. The formulae 
for Q and Gp follow by direct calculations using Theorem 4.2 for Q. 

Theorem 5.30. Let 7,r be related by (5.86) and let dp^^\dp^^^ be the images of 
dpFjdpK {for 7) under the map p t-^ p given by (5.87). Let f^^'-'^^z) be the Fade 
approximants for J2'^=o(~^)^^nZ'^ ■ Then 

hm /[2A^-i'2N](^)= f^P^^M (5.93) 

li^ ^~[2iV,2iV+l](^)^ [ dp^'^Hx) ^ 
N^co J 1+XZ 

In particular, if the T problem is indeterminate, lim/['^~-'^'^l(z) does not exist. 

Proof. This follows directly from Theorem 6 for the 7 problem and (5.88)-(5.92). D 

Of course, in the special case, the non-convergence of /[^~^'^1 (2) is replaced by some- 
thing almost as nice since the even and odd subsequences separately converge. But 
there is no reason to expect this to persist if F is not even and indeed, simple numeri- 
cal examples ([14]) illustrate non-convergence is the rule. One can understand why this 
happens. In the indeterminate case, there is a one-parameter family of extensions, each 
with discrete spectrum. There is, in general, nothing to tack down Ap so it wanders. In 
the even case, two self-adjoint extensions are even, one with an eigenvalue at and one 
without, and the Ap for odd A^ and even A^ track these. This suggests that A^^ , which 
is tacked down to have a zero eigenvalue, should converge in strong resolvent sense and 
the corresponding /l^'^l's should converge. This is true and is not a previously known 
result: 
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Theorem 5.31. Let {7n}5^o ^^ ^ ^^^ of Hamburger moments and let X]^o(~-'-)"^T"-'^"^ 
be the associated series of Hamburger. Then the /^^'^■'(-z) Fade approximants exist if 
and only if 

Pn{0)j^O. (5.95) 

Moreover, (5.95) holds if and only if there exists un so that the matrix A^j^ given by 
(5.2b) has a zero eigenvalue, and then 

/[^'^l(^) = (5o,(l + ^At?^+'^)-i5o). (5.96) 

If N ^> oo through the sequence of all N 's for which (5.95) holds, then 

lim f[N,N]^^^^ f Mxl ^5_g^^ 

iV^oo ^ ' J 1 + XZ ^ ' 

JV:Pjv(0)/0 

where dp is the unique solution of the moment problem in case it is determinate and the 
unique von Neumann solution with /u({0}) > if the moment problem is indeterminate. 

Remarks. 1. If Pn{0) = 0, then, by the second-order equation Pn{0) obeys, we have 
that P/v+i(0) 7^ 7^ Pn-i{0), so (5.95) can fail at most half the time. Typically, of 
course, it always holds (e.g., under a small generic perturbation of the Jacobi matrix, 
one can prove that it will hold even if initially it failed). 

2. If the 7's are a set of even moments, then (5.95) holds exactly for A^ even and 

.[2M,2M]/ X ^ ^ V2M+U--J 
^ ^ ^ Z^^^^P2M+l{-\) ' 

SO we get convergence of those measures we called dp above. 

3. The non-existence of /I^'-^l is using the Baker definition that we have made. There 
is a classical definition which would define f^^'^\z) even if P;v(0) = 0, (but (5.26) and 
(5.27) would both fail!). Under this definition, f^^^^\z) would be f[N-i,N-i]^^-^ ^^en 
Pn{^) = and (5.97) would hold using all A^ rather than just those A^ with Pn{^) 7^ 0. 

Proof. Define 

Mj,{z) = Pn-i{0)Pn{z) - Pn{0)Pn-i{z) (5.98) 

Nn{z) = PN-iiO)QN{z) - Pn{0)Pn-i{z). (5.99) 

Then by the Wronskian calculations (see Proposition 4.1); also see (5.70)), 

iV^(O) = ^ (5.100) 

ttN-l 

MUO) = J2 ^^"^ ^ (5.101) 
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and, of course, 

Miv(O) = 0. (5.102) 

Define 

An{z)=z^Nn+i{-^J (5.103) 

BN{z)=z''+'MN+^{-l). (5.104) 

By (5.100), deg^AT = A^ and by (5.101/5.102), degSjv = A^. Moreover, by tlie definition 
(5.98), 

SAr(O) = ^^ P7v(0) = (5.105) 

and 

Bn{0) = ^ P'^iO) ^ (5.106) 

since Pn{0) = ^ Pn-iW t^ 0. 
Moreover, we claim that 

2N 

Am{z) - Bn{z)J27jz' = 0(^2^+i). (5.107) 

This can be seen either by algebraic manipulation as in the remark following the proof of 
Theorem A. 5 or by noting that we have proven it in (5.39) for Stieltjes measures. So if 
p G A^^(7) and po is a Stieltjes measure and pt = {1 —t)po + tp, then Bn{—z), An{—z) 
are real analytic functions oft and (5.107) is true for all small t in (0, 1) since Ap (t) > 
for small t (and fixed A^) and that is all the proof of (5.106) depends on. 

If Pm{0) ^ 0, B^iO) ^ 0, so S[^'^!(z) = 1^ and A[^'^l(z) = ^ yield denomi- 
nator and numerator of f'-^'^^t); so /[-^'-^l exists. 

Moreover, the proof of Lemma 5.1 shows that an a^r exists giving Aj^ a zero 
eigenvalue if and only if Pn{0) 7^ 0. In that case, the proof of Proposition 5.8 is valid, 
and we have that (5.17) holds. Noting that ^f- = ^^, we see that (5.96) is then valid. 

On the other hand, suppose that P/v(0) = and that f^^'-^\z) exists. Then by 
(5.107), we have that 



A^(^)s[^'^i - s^(^)A[^'^i(^) = 0(^2^+1; 



and so it is zero since the left side of this is of degree 2A^. 

Since Pn{0) = 0, Mn+i{z) = -PN+i{0)PNiz) = cMn{z) where c = -P7v+i(0) 
Pn-i{0)~^, and similarly, N]\[^i{z) = cNn{z). It follows that 

An{z) ^ czAn-i{z), Bn{z) = czBn-i{z). 
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Thus, if PAr(O) = but f^'^\z) exists, then 

that is, the [A^ — 1,A^ — 1] Pade approximant would need to have an error of order 

Since Pn{Q) = 0, A^^' = Ap and so we would have that (5o, {Ap ySo) = 'jj for 
J = 0, 1, ... , 2N. But as noted in the remark following the proof of Proposition 5.11, this 
can never happen. We conclude that if PAr(O) = 0, /t^'-^] must fail to exist. 

To complete the proof, define Ax to be A if the problem is determinate and the 
unique von Neumann extension with eigenvalue zero if the problem is indeterminate (so 
Ai^7r(0) = 0). Then the proof of Proposition 5.3 goes through without change and shows 
that ^^ converges to A^: in strong resolvent sense. This implies (5.97). D 

As for other fi^+^~^'^\(^z) for series of Hamburger, we can only define "f^^' moment 
problems for j an even integer (positive or negative). This means that we can use 
(5.49)/(5.59) if £ is an even integer and (5.50)/(5.60) if £ is an odd integer. The result is 
the following: 

Theorem 5.32. Let {'^n}'^=o be a set of Hamburger moments and let Yl'^=oi~'^)^1nZ'^ 
be the associated series of Hamburger. Then: 

(i) If i is any odd integer, the Pade approximants f\-^+^~^'^'i(^z) exist if and only if 
P^)(0) ^ where fori>0, j =£- 1, and M = N and if £ < 0, j = £ + 1 and 
M — N + £ — 1. The functions that exist converge to a finite limit fe{z). 

(ii) If £ is any even integer, the Pade approximants fl^+^~^'^^(^z) all exist and lie in 
a compact subset in the topology of uniform convergence on compact subsets of 
C+. Any limit point, fe{z), is associated to a von Neumann solution, p^^\ of the 
"f^^' moment problem via 

f,iz) = Y^i-iy^,z^ + i-iY^,z' r f'^'^ 

„_n J —OO 



j=0 



+ xz 



if£>0 and 



^-1 /.OO -7 -{£)(-\ ^ -1 



if £ < 0. In particular, if the 7*^^^ problem is determinate, f\-^+^~^'^'i(^z) is 
convergent. 
(iii) Let £ be an even integer. If the 'j^^' problem is determinate {as a Hamburger 
problem), then for £ > 0, f£+i{z) = fe{z) = ■ ■ ■ = fo{z) = f-i{z), and for £ < 0, 
fe-i{z) = fiiz) = ■ ■ ■ = fo{z) = fi{z). In particular, if \'jj\ < C^^'^jl, then all 
j![N+e-i,N]^^^ converges to the same £-independent limits as N ^ 00. 
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§6. Solutions of Finite Order 

In this section, we will continue the analysis of the indeterminate Hamburger moment 
problem. Throughout, {7n}5^o ^^^^ t)e a set of indeterminate Hamburger moments and 
p will be a solution of (1.1), the moment problem for 7. 

We have picked out the von Neumann solutions as coming from self-adjoint extensions 
of the Jacobi matrix, A. But in a certain sense, every p does come from a self-adjoint 
extension of A\ For let Hp = L'^{R,dp), let D{B) = {f e Hp \ J x^lfix)]"^ dp{x) < 00}, 
and let 

{Bf){x)=xf{x). 

Let C[X] be the set of polynomials in x which lies in Tip since all moments are finite, 
and let TYq = C[X]. Corollary 4.15 says that Ho = Ti-p if and only if p is a von Neumann 
solution. Let Ap = B \ <C[X]. Then to say p is a solution of the Hamburger moment 
problem is precisely to say that Ap is unitarily equivalent to the Jacobi matrix. A, 
associated to 7 under the natural map (so we will drop the subscript p). In that sense, 
S is a self-adjoint extension of A, but not in the von Neumann sense. 

B is minimal in the sense that there is no subspace of Tip containing TYq and left 
invariant by all bounded functions of B ( equivalent ly, left invariant by all {B — z)~^ , 
z G C\R or by all e*^'^, s G M). So (if p is not a von Neumann extension) we are in the 
strange situation where D{B) fl Hq is dense in Hq, B[D{B) fl Hq] C Hq, but Hq is not 
invariant for the resolvents of Bl 

It is not hard to see that the set of all solutions of the moment problem is precisely the 
set of all self-adjoint "extensions" of A in this extended sense, which are minimal (TYq is 
cyclic for B) and modulo a natural unitary equivalence. This is a point of view originally 
exposed by Naimark [24, 25, 26, 27] and developed by Livsic [21], Gil de Lamadrid [7], 
and Langer [20]. See the discussion of Naimark's theory in Appendix 1 of [2]. 

In the language of Cayley transforms, the extensions associated to p's that we will 
call of order at most n below are parametrized by unitary maps, U, from /C_|_ © C"^ to 
/C_ © C"^ with Ui, U2 equivalent if and only if there is a unitary map V : C^ -^ C^ so 
Ui{l © y) = (1 © V)U2- From this point of view, these extensions are then parametrized 
by the variety of conjugacy classes of U{n + 1) modulo a U{n) subgroup. This has 
dimension (n -(- 1)^ — n^ = 2n + 1. Our parametrization below in terms of a ratio of two 
real polynomials of degree at most n also is a variety of dimension 2?i + 1, showing the 
consistency of the parametrization. (The p's of exact order n will be a manifold.) 

Given any p, we define the order of p, 

ord(p) = dimCHp/Ho) 

so the von Neumann solutions are exactly the solutions of order 0. Our main result in this 
section will describe all solutions of finite order which turn out to correspond precisely 
to those p's whose Nevanlinna function $p are ratios of real polynomials. 

They will also be p's for which (ip/ nr=i 1-^ ~ ^*l^ ^^ ^^^ measure of a determinate 
moment problem. We will therefore need some preliminaries about such problems. Given 
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zi, . . . ,Zn G C+, by using partial fractions we can write for x real, 



nn I 
,^1 \x - Z^ 



12 i^^ra — 1n\Xi ^Ij • • • i ^n) ~r / 



" " ^7^ 1 ^f 



Z i "^ Z 2 2 2 



(6.1) 



2 



where Pm-2n is a polynomial in x of degree m — 2n. To see (6.1), write ni=i \^ ~ ^i\ 
Y[7=ii^ — Zi){x — Zi) and analytically continue. The sum in (6.1) comes from computing 
the residues of the poles of these functions. Define 

T^^\zi,...,Zn)=E2{Pm-2n{x;Zi,...,Zn)), (6.2) 

where, as usual, E2 is the expectation with respect to any solution of the moment 
problem. For (i, . . . Xn G C+, let 



Tm{zi,...,Zn;Cl^---Xn) = T^^ (^1 , . . . , ^n) + ^1 



1=1 






(6.3) 



Theorem 6.1. Fix z\, . . . ,Zn and Ci, • • • , Cn in C+. There is a one-one correspondence 
between v G M.^ {T{z\^ . . . , Zn\ Ci, • • • , Cn)) and those p in Ai^ {'^) which obey 

^^ ^- = C^, z = l,...,n (6.4) 



X 



under the association 



n 



dp{x) *-^ diy{x) = Y]\x — Zi\ '^ dp{x). (6.5) 

That is, p has moments {7m}m=o '"^^^^ subsidiary conditions (6.4) z/ and on/y z/z/ ^fiven 
by (6.5) /ias t/ie moments {r^(2i, . . . , ^n; Ci, • • • , Cn)}5^=o- 

Proof. This is pure algebra. The functions {xm}m=o^ i x-z- }?=i ^^^ linearly independent 
as analytic functions (finite sums only). (6.1) says that each function x'^/YYi=i{x — 
Zi){x — Zi) = Qm{x) is in the span, 5, of these functions. On the other hand, since 
nr=i(^ ~ ^i)i^ ~ ^i) = -^(^) is a polynomial in x, and x^Qm{x) = Qe+rn{x), L{x)Qm{x) 
is in the span of the {Qmix)}. Similarly, since / _^\ is a polynomial, L{x)Qo{x) t^l^.) 
lies in that span. 

Thus, {Qm{x)} is also a basis of S. It follows that there is a one-one correspondence 
between assignments of numbers Qm{x) i— *> F^ and of assignments x"^ i— *> 7^, -^- 1— *> Q, 

— ^ ^-* Xi since each defines a linear functional on S. For the linear functional to be 

X — Zi ' 

real in F^, representation, all F^, must be real. For reality in (7, C, A) representation, 
we must have 7^ real and A^ = Q- (6-3) is just an explicit realization of the map 
from 7, C, ( to F. Thus, a real dp has moments 7^ with Q conditions if and only if 
/ x™" Yli=i 1"^ ~ -2i|~^ dp{x) = Tm for all m. Since (6.5) shows u is positive if and only if 
p is, we see that the claimed equivalence holds. D 

As an aside of the main theme of this section, we can construct determinate Hamburger 
moments quite close to indeterminate ones. The result is related to Corollary 4.21. 
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Theorem 6.2. Let {7n}5^o ^^ ^ ^^^ of Hamburger moments. Then for any a G (0, oo), 
there is a Ca so that the moment problem with moments 

72n+l = 

72n = 72n ' a^l2n-2 + aSn-2 + ' ' " + {-iVa^'^lO + (-l)"+^C„a2-+2 

(6.6) 

is a determinate moment problem. 

Proof. Let pi be any solution for the Hamburger problem. Since ^ [dpi{x)+dpi{—x)] = dp 

has moments 

( m odd I 

X dp{x) = < = 72^, 

I 72m m even 

we can suppose 72m+i = 0. Let p be a von Neumann solution to the 7^ moment 
problem which is invariant under x -^ —x. If the 7^ problem is determinate, then the 
unique solution is invariant. Otherwise, this follows from the proof of Theorem 2.13. 
Let dp{x) = if^^li ■ By Corollary 4.15 and Theorem 6.1, p is the unique solution of 

the 7' moment problem with J _^.^ _i = Ca-i- Thus, the T m{ia~^ , Ca-^) problem is 
determinate. If Ca = a~^ImCa-i, then a calculation using the geometric series with 
remainder shows that r^(— za~^, Ca-O = 7m given by (6.6). D 

Returning to the main theme of this section, we next examine when a Herglotz function 
is a real rational function. Define for z G C_|_, 

^z{x) = — - — (6.7a) 

viewed as a continuous function on M U {00} with 

iPz{oo) = z. (6.7b) 

If da{x) = [ -.^ 2 ] + c5 00 as a finite measure on MU {00}, then (1.19) can be rewritten as 



$(z) =d+ [ 



LPz{x) da{x). (6.8) 



If a{{oo}) = and / \x\ da{x) < 00, define d = d — j x da{x). 

Proposition 6.3. A Herglotz function is a ratio of two real polynomials if and only if 
the representation (6.8) has a a with finite support. If a has exactly N points in its 
support, there are three possibilities for the degrees of P,Q in ^{z) = -J^ with P,Q 
relatively prime polynomials: 

(i) (a({oo}) =0,(i = 0), degP = A^-l, degg = A^ 

(ii) (cT({oo})=0,(i^O), degP = degg = A^ 

(iii) ((7({cx)})^0), degP-iV, degg = iV-l 
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In all cases, we say N = max(deg(P), deg((5)) is the degree of^. 

Proof. Elementary. D 

Remark. The set of cr's with exactly A^ pure points is a manifold of dimension 2N (N 
points, N weights), d G M is another parameter so the set of such $'s is a manifold of 
dimension 2N + 1. 

As a final preliminary, for any z G C\a{B) = C\ supp{dp) and n = 1,2,..., we 
introduce the functions en{z) on supp((i/9), 

en{z){x) 



{x — z)"^ 



thought of as elements of Tip = L^i^, dp). 
The main result of this section is: 

Theorem 6.4. Let •y be a set of indeterminate Hamburger moments and let po ^ 
^4-^(7). Fix N E {0, 1, . . . }. Then the following are equivalent: 

(1) po has order at most N . 

(2) For some set of distinct zi, . . . ,zn G C\M, TYq U {e-i{zj)}^^^ span Tipf^ . 



(3) For any set of distinct zi, . . . ,zn G C\]R, Hq U {ei{zj)}jL-f^ span Ti 

(4) For some zq G C\M, Hq U {6^(20) }Li span H 



Po- 



(5) For any zq G C\R, Hq U {ej(2o)}jLi span H-p^. 

(6) For some set of zq, zi, . . . , z^ G C+ and Q = Gpg{zj), there is no other p G 
M"{-f) with Gp{zj) =Cj. 

(7) For all sets of zq^zi, . . . ,zn G C_(_ andQj = Gp^^{zj), there is no other p G ^^■^(7) 
with Gp{zj) = Qj. 

(8) For some zo,zi,...,zn G C+ and Cj = Gp„{zj), the r(2o, . . . , ^tv; Co, • • • , Cn) 
moment problem is determinate. 

(9) For any zq, ^i, . . . , ^at G C+ and Q = Gp^^ (zj), the T{zo, . . . , ^at; Co, • • • , Cn) mo- 
ment problem is determinate. 

(10) For some zi, . . . ,zn G C+ and Q = Gp^^{zj), po is a von Neumann solution of 
the T(zi, . . . , zn'i Ci, • • • , Cn) moment problem. 

(11) For any zi, . . . , zn G C_|_ and Q — Gp^^{zj), po is a von Neumann solution of the 
r(zi, . . . , zn; Ci, • • • , Ca^) Tnoment problem. 

(12) The Nevanlinna function $ of pq is a rational function of degree at most N . 

In particular, the Nevanlinna function $ of po has degree N if and only if po has 
order N. Moreover, if one and hence all those conditions hold, po is an extreme point in 
A4^ (j) and is also a pure point measure. 

Remark. The measures obeying (1)-(12) are what Akhiezer calls canonical solutions of 
order A^. 

Lemma 6.5. If ^ E Hq and z ^ o'{B), then 

{B-z)-^ipeno + [ei{z)]. 
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Proof. Since {B — z)~^ is bounded and Hq + [ei(t)] is closed, it suffices to prove this for 
if{x) = P{x), a polynomial in x. But then 

{B - z)-\ ={x- z)-^P{x) = P{z){x - z)-^ + (-P(^) - -P(^)) 

X — z 

^ P{z)ei{z) + R{x) 

for a polynomial R, that is, {B — z)~^Lp G Hq + ei{z). D 

Lemma 6.6. Supose that zq.zi, . . . ,zi G C\(j(i?) are distinct. If e^(2o) is in the span 
of Ho U span[{ej(2o)}^"i^] U span[{ei(2A:)}i=i], then so is em+i{zo). 

Remark. £ can be zero. 

Proof. By hypothesis, for some if G Hq and {cij}"^^ and {l3k}i=i in C, 

m-l £ 

em{zo) = 9^+ X^ "jej(2o) +'^Pkei{zk). 
j=i fc=i 

Apply (S — zo)~^ to this using {B — zo)~^ej{zo) = ej+i{zo), Lemma 6.5, and {B — 
zo)~^ei{zk) = {zo — Zk)~^[ei{zo) — ei{zk)] and we see that em+i{zo) is in the requisite 
span, n 

Proposition 6.7. Let i = ord(p) < oo. Then for any set {zi}f^i of i distinct points in 
C\a{B), Ho U {ei(2,)}Li span Up. 

Proof. By the Stone- Weierstrass theorem, linear combinations oi {{x — z)~^ \ z G 
C\(T(i?)} are dense in the continuous functions on (j{B) vanishing at infinity, and so 
they are dense in L^(M, dp). It follows that we can find {zi}f^i so that Hq U {ei{zi)}^^i 
span Tip. 

Pick Zo G C\a{B) and let W = {B - Z(){B - zq)'^. Then W is bounded with 
bounded inverse, so it maps dense subspaces into dense subspaces. By Lemma 6.5, 
W[no\ cno + [eiizo)]. For i = 1, ...,£- 1, Wciiz,) = ei(^,) + {^l^)[e^izo) - ei(^,)] 
and Wei{zi) = ei{zQ). Thus, W maps the span of {ei{zi)Yi^i into {ei(2i)}j~Q. So 
7ioU{ei(2i)}j~Q span Tip. By successive replacement, we can move the ZiS to an arbitrary 
set of distinct points. D 

Remark. Proposition 6.7 shows (l)-(3) of Theorem 6.4 are equivalent. 

Proposition 6.8. Suppose that £ = ord(p) < oo. Then for any zq G C\a{B), Hq U 
{ej(2o)}j=i spanHp. 

Proof. By hypothesis, there must be some dependency relation 

V7 + 5^«jej(2o) =0 (6.10) 
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for if G Tio and some (cti, . . . , a^+i) 7^ 0. Let /c + 1 = max{j | a^ 7^ 0}. Then solving 
(6.10) for efc_(_i(2o), we see that efc_|_i(2o) hes in the span of Tio U {ej(-2o)}j=i; and so by 
induction and Lemma 6.6, dX\em{zQ) he in this span and so in the spanof 7ioU{ej(zo)}^^i. 
61(2;) is an analytic function in C\M with Taylor coefficients at zq equal to en{zo), so ei{z) 
lies in the span of 7YoU{ej(2;o)}j=i for z in the same half plane as zq. By Proposition 6.7, 
these ei's together with TYq span Tip. D 

Remark. Proposition 6.8 shows that (1), (4), (5) of Theorem 6.4 are equivalent. 

Proposition 6.9. Fix distinct points zi, . . . ,zi E C\]R. Let 

i 
d\i{x) = TT \x — Zi\~'^ dp{x). 

i=l 

Then the polynomials are dense in L^(]R, d^) if and only ifHo U {ei{zj)Yj^i span Tip. 
Proof. Let U be the unitary map from Tip to Tip, given by 



{Uf){x) = \[{x-z,)f{x). 



i=l 

Then t/ maps the span of C[X]U{ei(2i)}f^]^ onto C[X] for P(x) is a polynomial of degree 
m if and only if ni=i(^ ^ z)~^P{x) is a linear combination of a polynomial of degree 
max(0,?7i — i) and {{x — -2i)~^}i=i- Thus, the density assertions are equivalent. D 

Remark. Proposition 6.9 and Corollary 4.15 show the equivalence of (2), (3), (10), and 
(11) of Theorem 6.4 (given that we already know that (2) and (3) are equivalent). 

Proposition 6.10. // ord(po) ^■s finite, then po is a pure point measure with discrete 
support. 

Proof. Let i = ord(po)- Pick distinct {zi}^^i in C+. Then by Propositions 6.7 and 6.9, 
the polynomials are dense in L'^(R, dpi) where for /c = 0, 1, . . . , £, we set 

k 

d^kix) = TT \x — Zi\~'^ dp{x). 
j=i 

Let ko be the smallest k for which the polynomials are dense in L'^(M.,dpk)- If ^0 = 0, 
then dp is a von Neumann solution of the 7 problem. If /cq > 0, then dfXko-i is not a 
von Neumann solution of its moment problem. So by Proposition 4.15 and Theorem 6.1, 

the moment problem for Trn = Tm{zi, . . . , -Zfco-i, Zk„'i Ci, • • • , Cfco-i? Cfc) is indeterminate. 
But since the polynomials are dense in L^(R, (i^Ufc), dfiko is a von Neumann solution of 
an indeterminate problem. 

Either way, dfiko is a von Neumann solution of an indeterminate problem. It follows 
Theorem 4.11 that dpko and so dp is a discrete point measure. D 
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As a final preliminary to the proof of Theorem 6.4, we need a known result about inter- 
polation of Herglotz functions. Given zi, . . . , Zn] wi, . . . ,Wn G C+ with the ^'s distinct, 
define the n x n matrix D by 



Dij{zi,...,Zr,;Wi,...,Wn) 



Wi — w 



J 



Zj 



(6.11) 



Theorem 6.11. Pick zi, 
with 



, Zn G C_|_ with n > 1. There exists a Herglotz function $ 
^z,) = w,, z = l,...,n (6.12) 



if and only if D{zi, . . . ,Zn]Wi, . . . ,Wn) is a {not necessarily strictly) positive definite 
matrix. Moreover, the following are equivalent, given such a $: 

(1) del D{zi,...,Zn;wi,...,Wn)=Q. 

(2) There is a unique Herglotz function $ obeying (6.12). 

(3) ^ is a real rational polynomial of degree at most n — 1 . 

Remark. Since (3) is independent of the choice of {-2i}"=i, so are (1), (2). 

We will sketch a proof of this result (fleshing out some arguments in [1]) below. 

Proof of Theorem 6.4- As already noted. Propositions 6.7, 6.8, and 6.9 show that (1), 
(2), (3), (4), (5), (10), and (11) are equivalent. Theorem 6.1 and Corollary 4.15 prove 
the equivalence of (10) and (11) with (6), (7), (8), and (9). By Theorems 4.14 and 6.11, 
(6) and (12) are equivalent. 

Since ord(p) = A^ is equivalent to ord(p) > N and the negation of ord(p) < A — 1, we 
see that $pg has degree precisely equal to ord(po)- That po is then an extreme point is 
proven in Appendix B. D 

We conclude this section by proving Theorem 6.11. 

Proposition 6.12. Let ^ be a Herglotz function and let Wi = ^{zi), i = l,...,n for 
distinct zi, . . . ,Zn G C+. Let D be given by (6.11). Then D{zi, . . . , Zn', wi, . . . , Wn) is 
{not necessarily strictly) positive definite and det D = if and only if ^ is a real ratio 
of polynomials of degree n — 1 or less. 

Proof. Using the Herglotz representation (1.19), we get a representation of Dij, 

Wi — Wi f , , . 1 1 



c + / d^{x) 



X 



— 1 



from which we obtain for a G C 



N. 



J2 cx^aJD 



*J 



E 

i=l 



a. 



+ 



dn{x] 



E 

i=l 



ai 



Jyi 



proving the positivity. 



(6.13) 
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Suppose det{D) = 0. Then there is a non-zero a G C"^ so that the right side of (6.13) 
is 0. Note that ^^^^ (x-\) ~ YYi=i(^ ~ Zi)~^Q{x), where Q is a polynomial of degree 
n — 1 if X]r=i <^i 7^ ^^^ of degree at most n — 2 if X]r=i '^i ~ ^- Thus, for the right side 
of (6.13) to vanish for some non-zero a in C", either c 7^ and fi is supported at n — 2 
or fewer points, or else c = and fx is supported at n — 1 or fewer points. Either way, by 
Proposition 6.3, $ is a real rational function of degree at most n — 1. 

Conversely, suppose $ is a real rational function of degree precisely n — 1. Then either 
H is supported at n — 1 points (say, xi, . . . , Xn-i) or c 7^ and // is supported at n — 2 
points (say, Xi, . . . , Xn-2)- The map i/j : C^ ^' C^~'^ given by 



t/jjia) 



n 



Xj Zi 



(with 'i/;„_i(a) = ^^^^ a^ if c 7^ 0) has a non-zero kernel by dimension counting. Thus, 
the right side of (6.13) is zero for some a G C^. If degree $ is smaller than n — 1. 



we can find (ai,...,^^) G C^ so that X^i 1=1 <^j<^j-^ij = O5 which still implies that 
det(L>) = 0. D 

Proposition 6.13. Let ^i, . . . , 2„; wi, . . . , w„ lie in C_)_ with the Zi distinct. Suppose 
the matrix Dij given by (6.11) is {not necessarily strictly) positive definite. Then there 
exists a Herglotz function with ^{zi) = Wi, i = 1, . . . ^n. 

Proof. This is a fairly standard use of the Hahn-Banach theorem. We will use the 
representation (6.8) for $ and the function tp^ of (6.7). Begin by noting that without 
loss of generality, we can suppose Zn = i (by mapping C+ -^ C+ with a linear map that 
takes the original z^ to i) and that Rety^i = (by adjusting the constant d in (6.8)). 
Indeed, since 

^^{x) = z, (6.14) 

this choice Re$(z) = is equivalent to (i = in the representation (6.8). 

Let V be the 2n — 1-dimensional subspace of C{M. U {00}), the real- valued continuous 

functions on M U {00} spanned by 1, {Re(pzj{-)}]Zi ^^'^ (I'^'/'^i ( ' )}j'=i • ^^^y / ^ ^ 
can be written: 

n-l 

f = A^ + J2 ^jVj(a^) + Aj ^j{x) (6.15) 

with Ai, . . . , An-i G C and A^ G M. Define a linear functional / : 1/ ^ M by 

n-l 

i(f) = Anlmwn + ^ AjVUj + AjWj (6.16) 

if / has the form (6.15). 
We will prove shortly that 

fix) > for all X =^ £(/) > 0. (6.17) 
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Assuming this, we let X = {f e C{R U {00} | f{x) > for all x} and Y = {f e V \ 
£{f) = 0}. Since X is open and Y closed, by the separating hyperplane version of the 
Hahn-Banach theorem (see [33]), there is a linear functional L on C(M U {00}) so L > 
on X and L < on y. If L is normalized so -^(1) = Imwn, it is easy to see that L 
extends i, and so defines a measure a on M U {cxd} with 

^{zi) = / if^^ (x) da{x) =Wi, i = 1, . . . , n. 

Thus, if (6.17) holds, we have the required $. 

Since f{x) > on the compact set M U {00} implies f = e + g with g > and £ > 
(by hypothesis, 1 E V) and ImWn > 0, we need only show / > implies £{f) > 0. 

If / has the form (6.15), we can write 



nn — l I 19 

, = 1 \x-Zj\^ 



where Q is a polynomial. / > implies Q > 0. Any real non-negative polynomial has 
roots either in complex conjugate pairs or double real roots, so Q must have even degree 
and have the form 



£ 



Q{x) = c^W \x-yk\'^ 



k=i 
for suitable yi, . . . , y^ G M U C+. Thus, with R{x) = cY[k=i{x ~ yk)^ 

f=\h\' with h{x) = f/^^ . 

Since / is bounded on M U {00}, deg{R) < n — 1, and thus h has the form 

n-l 

i=i 
Since {x — Zj)~^ = {zj — '>')~'^[ (^Iz ) ~ ■'■]' ^^ '^^^ rewrite (6.18) as 



■{X-Zj) 

h{x) = y^ a. 
for real x. Thus, 

'^zj{x) -^zi{x) 



X — I 

J 
X — z^ 
i=i ^ 



X +1 

fix) = \h{x)\^ = ^a^a^. ______ = ^a,a, 



Zj Zi 
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SO that 

n 

^(/) = X^ aiajDij, 

which is non-negative by hypothesis. D 

Proof of Theorem 6.11. The first assertion is a direct consequence of Propositions 6.12 
and 6.13. Proposition 6.12 shows the equivalence of (1) and (3). 

To prove (1) ^ (2), suppose (1) holds and, as in the last proof, we can suppose that 
Zn = i and Re Wn = 0. Then by our proof of Proposition 6.12, the set of points where the 
measure a of (6.8) is supported is determined by the a with ^^ Dijaj = as an n — 1 
point set xi, . . . , Xn-i (with Xn-i = cxo allowed). Thus, any such $ is of the form P/Q 
where deg(P) < n — 1 and Q{z) = Ylj=i (^ ~ ^j) (if ^n-i = oo, the product only goes 
from 1 to n — 2), If $i and $2 are two solutions of (6.12), they have the same Q but they 
could be distinct P's, say Pi and P2- But by (6.12), Pi{z) — P2{z) vanishes at 21, ... , Zn- 
Since Pi — P2 is a polynomial of degree n — 1, this is only possible if Pi — P2 = 0, that 
is, $1 = $2. Thus, (1) => (2). 

For the converse, suppose the determinant 

D{zi,...,Zn;wi,...,Wn) >0. (6.19) 

If 71 = 1, since z,w E C+, it is easy to see there are multiple $'s obeying (6.12). (For 
example, as usual we can consider the case z = w = i and then that ^{z) = z or 
$(^) = —z~^.) So we suppose n > 2. 

Consider the function g{w) = det[-D(2;i, . . . , ^^j,; t^i, . . . , Wn-i, w)] where we vary Wn- 
g{w) is of the form c\w\'^ + dw + dw + e with c < 0, e real, and d complex, (c is strictly 
negative since it is — det{D{zi, . . . , Zn-i; t^i, . . . , Wn-i) and we are supposing (6.19) and 
n > 2.) Thus, the set g{w) > is a disk and since g{wn) > 0, Wn is in its interior. Let wq 
be the center of this disk, R its radius, and let w{6) — wq + Pe*^. Since we have proven 
(1) =^ (2), there is a unique $ with 

^g{z,)=w,, i = l,...,n-l; ^0{zn) = w{e) (6.20) 

and it is a rational function of degree at most n — 1. 

As usual, we can suppose zi = wi = i so each ^0{z) has the form 

M^)= f^^-^daeix), (6.21) 

J X - z 

where dag is a probability measure on M U {00}. If 6*^ ^ 6*00, and dag^ -^ dp, then the 
Herglotz function associated to dp obeys (6.20) for 6*00 and so it must be ^e^. It follows 
(since the probability measures are compact) that dp = dag^ and thus {dag} is closed 
and dag 1— > ^g is continuous. Since a continuous bijection between compact sets has a 
continuous inverse, 9 1— > dag is continuous. 
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Since $e is unique, dae is a pure point measure with at most n — 1 pure points. Since 
the function is not determined by (zi, . . . , Zn-i), there must be exactly n — 1 points. It 
follows that the points in the support must vary continuously in 9. 

Note next that for any 9, we can find 9' and tg G (0, 1) so that 



Wr 



^tew{9) + {l-te)w{9'). 



Suppose that there is a unique $ obeying (6.12). It follows that for any 6*, 

^{z)=te<^0{z) + {l-tB)^B'{z) (6.22) 

and thus ^{z) has a representation of the form (6.21) with da a point measure with at 
most 2n — 2 pure points. By (6.22) again, each da^ must be supported in that 2n — 2 
point set, and then by continuity, in a fixed n — 1 point set. But if dag is supported in 
a fixed 6'-independent n — 1 point set, so is da, and thus (6.19) fails. We conclude that 
there must be multiple $'s obeying (6.12). D 



Appendix A: The Theory of Moments and Determinantal Formulae 



The theory of moments has a variety of distinct objects constructed in principle from 
the moments {7n}5^o- ^^^ orthogonal polynomials {Pn{x)}'^^Q, the associated polynomi- 
als {Qnix)}'^^Q, the sums X]?=o-^j(0)^ ^^^ X]?=o^j(0)^' ^^^ Jacobi matrix coefficients. 



and the approximations 



.(^) 



and 



N„{x) 



It turns out most of these objects can 



P„(x) "^^^ M„(x)- 

be expressed as determinants. These formulae are compact and elegant, but for some 
numerical applications, they suffer from numerical round-off errors in large determinants. 
We have already seen two sets of determinants in Theorem 1 . Namely, let Ti^ be the 
N X N matrix, 

/ 7o 7i • • • iN-i \ 
7i 72 • • • In 



Hn 



V 



(A.la) 



7iv-i In 



727V-2 



/ 



and Sn, the matrix 



S 



1 71 72 

72 73 



TV 



IN \ 

7iv+i 



V 



In 7^+1 



72^- 



J 



and define 



hjq = det(7i 



N 



sn = det {Sr, 



(A.lb) 



(A.2) 



(for comparison, D^ from [1] is our hN^i). We will use /iAr(7) if we need to emphasize 
what moments {7n}5^o ^^^ involved. Thus, Theorem 1 says Hn > for all A^ is equivalent 
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to solubility of the Hamburger problem and hpj > 0, sn > for all A^ is equivalent to 
solubility of the Stieltjes problem. 

There is an interesting use of determinants to rewrite Hn and sn in terms of the 
moment problem that makes their positivity properties evident. Suppose dp obeys (1.1). 
Then 



7V-1 



N-1 



hN = det{{x'^+^)o<a,b<N-i) n ^/^(" 

■^ a=0 

/- n -| 7V-j_ 

Y[x'^ det{{xl))o<a,b<N-iY[dp{xa). 
-a=l -I a=0 

Permuting over indices, we see that in na=i('^a) det{x'^), we can replace Xa by x-n-(^a) for 
any permutation, n. Since det(a;^. n) = (— 1)^ det(x|^), we see that 

„ N-l 

hN = (Nl)-' [det{{x%<a,b<N-i)f n ^'^(^«)- 

J „ — n 



Recognizing the Vandermonde determinant, we have 



hN = {N\) 



n-i 



n (^ 



a=0 



N-l 



0<a<b<N-l a=0 



W dp{xc 



(A.3a) 



Similarly, 



N-l 



Sn = {N\) ^ I Xq.-.Xn-I JJ {Xa~Xb] 



l[dp{ 



Xa) 



(A.3b) 



0<a<b<Af-l 



a=0 



The most basic formula is: 
Theorem A.l. Pn{x) is given by 

Pn{x) 



\/hnh 



det 



n+l 



/ 70 7i 

71 72 



In-l In 
\ 1 X 



In \ 

ln+1 



72n-l 

x"^ 



(A.4) 



Proof. Let Sn{x) be the determinant on the right side of (A.4). Then for any solution p 
of the moment problem, we have that J Sn{x)x^ dp{x) is given by the same determinant, 
but with the last row replaced by 7j 'jj+i • • • 7j+n- It follows that 



x^ Sn{x) dp{x) = 0, 

— hn+1 



j = 0, l,...,n- 1 
j = n. 



(A.5) 
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In particular, since 



we see that 



Sn{x) = hnX^ + lower order, 

Sn{xY dp{x) = hnhr 



From (A.5)-(A.7), it follows that P„(x) 



Snjx) 



Y "Ti'ln 



i'n+1- 



D 



(A.6) 
(A.7) 



+1 



From (A. 4), we can deduce formulae for the coefficients an,bn of a Jacobi matrix 
associated to the moments {7n}5^o- ^^^ ^n be the n x n determinant obtained by 
changing the last column in Hn- 



( 



hn = det 



7o 7i 



\7n !■ 



n+l 



7n-2 In \ 

72n-3 72n-l/ 



Thus (A. 4) implies that 



Pnix) 



Theorem A. 2. For n > 0, 



hr 



h 



n+l 



h 
x" — -^ x^~^ + lower order 



(i) 
(ii) 



3=0 



hnhn+2 
'''n+l 
hn+1 



1/2 



h 



n+l 



Proof. By the definition of an and 6^, we have that 

xPn{x) = ttnPn+lix) + h^Pni^x) + tt^-lPn-lix). 

Identifying the a;""*"^ and x"^ terms in (A. 9) using (A. 8), we see that 



hr 



h 



ttr 



n+l 



hn 



+ 1 



hr 
hr 



'h 
hr 



hn+1 
hn+2 

hn+1 
n+2 \ hn+1 



n+l 



+ br. 



hr 



h 



n+l 



(A. 



(A.9) 



(A.IO) 
(A.ll) 
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(A. 10) implies (i) immediately, and given (A. 10), (A. 11) becomes 

^" = J. IT ' 

which implies (ii) by induction if we note the starting point comes from looking at the 
constant term in xPq{x) = x = ao-Pi(x) + 60-^0(2^)7 which implies that bo = ^ = —. □ 

Remark, (ii) has an alternate interpretation. (A. 8) says that jf"^^ is the sum of the n + 1 
roots of Pn+i^x). But Pn+i{x) is a multiple of the determinant of the Jacobi matrix 
Ap . So the sum of the roots is just the trace of Ap , that is, Y^^=o ^j- 

From (A. 4) and Theorem 4.2, (i.e., Qn{x) = Ex{ " x-y" ))' ^^ immediately get 
Theorem A. 3. 

/ 70 71 ••• In \ 

71 72 • • • 7n+l 



Qnix) 



1 



\/Kh 



det 



n+l 



7n-l 
\Rn,o{x) 



In 
Rn,l{^) 



72n-l 
J^n,n\Xj / 



(A.12) 



where 



i-i 

k=0 



J>1 



Proof. 



X-' —y 
x-y 



Rn,j=o{x) = 0. 

Ei=I)l/^"^"^2;'= SO Ex{ 



(A.13a) 
(A.13b) 



2n-l 



Remark. By (A.4), Pn{x)Y.kLo Ik 



X 



x-y 
k-1 



Rn,j{x). D 



has the same form as (A.4) but with the bot- 



tom row replaced by Sn,o{x) ■ ■ ■ Sn,nix) where 

'~^n,j \X 



2n-l 
k=0 



J-1 

E -tj-l-lX^ 

£=-{2n-j) 



= Rn,j{x)+X \j+X ^-fj + i + --- + X ""-fj+n-l+Oix "" ^). 

Recognizing x~'^'^j+x~'^'^j+i + - ■ ■ + x~"~^7j+n-i as x~^ times the first row of the matrix 
in (A.4) plus x~'^ times the second row plus . . . , we conclude that 

Pn{x) I J2 7fc^"'"' ) = Qn{x) + 0(x-"-l) 



fc=0 



and thus 



Qnjx) 
' Pn{x) 



2n-l 
fc=0 



-k-1 



+ 0{x 



-2T^-l^ 



consistent with the f^^ ^'^1 Fade formula (5.28). 

We saw the quantity L of (5.29) is important. Here is a formula for it. 
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Theorem A. 4. L — lim^ 



?«(o) 



P„(0) 



and 



<5n(0) tn 



where Sn = det{SN) and 



det 



Pn{0) Sr 



/ 70 71 

7o 71 72 



V 



In 1n+l 7n+2 



7n \ 

72n+l 



72r 



-J 



(A.14a) 



(A.14b) 



Proof. Follows by putting x = in our formula for Pn{x) and Qn{x). D 

Next we have explicit formulae for M„,{x) and Nn{x), the polynomials introduced for 
Section 5 (see (5.15) and (5.16)). 

Theorem A. 5. 



M^{x) 



1 



hr 



( 7i 72 

72 73 



Sn-l V ^n+1 



7n \ 

7n+l 



7n-l 7n-2 
2 



72n-2 
X^ I 



(A.15) 



iVn(x) 



1 



hr 



( 71 

72 



■Sn-l V "-n+l 



72 
73 



7n-l 7n-2 

\i?„,i(x) i?n,2(a;) 



7n \ 

7n+l 



72n-2 
^n,n\X) ' 



(A.16) 



where Rij{x) is given by (A. 13). 

Remarks. 1. Pn is given by a (n + 1) x (n + 1) matrix. To get M^, which has an n x n 
matrix, we drop the first column and next to last row. 

2. It is interesting that (A.15) does not obviously follow from the basic definition 

(5.15). 

Proof. The right side of (A.15) (call it M^(x)) has the following properties: 
(i) It is a polynomial of degree n. 

(2) It obeys M^{0)=^ 0. 

(3) It obeys E^i^x^ M^^x)) = for j = 0, 1, . . . , n — 2 since the corresponding matrix 
has two equal rows. 

(4) It obeys Mn{x) = ^/hnJK^+l x"^ + lower order so it has the same highest degree 
term as Pnix). 
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These properties uniquely determine Mn{x) so (A. 15) is proven. (A. 16) then foUows 
from (5.24). D 

Remark. By mimicking the argument foUowing Theorem A. 3, one finds that 

2n-2 

SO that 



A;=0 



(X 



Mr,{x) 



2n-2 



^^,^-^-i+0(x-2"), 



fc=0 



consistent with the jW-^^^-M Pade formula used in the proof of Theorem 6 (at the end 
of Section 5). 

Remarkably, there are also single determinant formulae for X]?=o-^j(0)^ ^^'^ X]?=o 

Theorem A. 6. We have that 



Y.PM 



2 '^ri 



J=0 



hn+l 



Y.QM 

J=0 



2 _ _Wn+2_ 
hn+1 



where Vn is given by the n x n determinant, 



v„ = det 



/ 72 

73 



Wn = det 



\7n-1 



7n+l\ 

7n+2 



\ln+l ■■■ 72n / 

and Wn+2 by the {n + 2) x {n + 2) determinant, 



/ 70 71 

70 71 72 
70 71 72 73 



7n-l In 7n+l 7n+2 



7n-l\ 

7n 
7n+l 



72r 



/ 



(A.17) 
(A.18) 



(A.19) 



(A.20) 



Remark. Thus by Theorems 3 and 7, indeterminacy for both the Hamburger and Stieltjes 
problems can be expressed in terms of limits of ratios of determinants. In the Hamburger 
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case, we need (A. 17) and (A. 18) to have finite limits in order that the problem be indeter- 
minate. In the Stieltjes case, (A. 17) and (A. 14) must have finite limits for indeterminacy 
to hold. 



Proof. We actually prove a stronger formula. Let 

N 



Bi,{x,y) = Y,Pn{^)Pn{y)- 



(A.21) 



n=0 



We will show that 



Bn{x,y) = -h^Ldet 



^01 X 

1 7o 7i 

y 7i 72 



x^ ^^ 

In 

1N+1 



\ 



y^ 7n 7n+i 



(A.22) 



72 JV 



/ 



(A. 17) then follows by setting x = y = and noting that if C = {Cij}i<i<j<N is an 
N X N matrix with Cn = 0, Cij = 52j, and Cn = 5^2, then 

det{{Cij)i<ij<]s[) = — det{{Cij)3<ij<N). 
(A.18) then follows from Qj{z) = ^^(I^^M^^PiM' 



(x-z) 



')^ 



since 



N 



Y, Qn{0)^ = E,Ey{[BN{x, y) - Bn{x, 0) - S^(0, y) + S^(0, 0)]x-^y-^). 

n=0 

Thus we need only prove (A.22). To do this, let B^^x^y) be the right side of (A.22). 
Consider Ej;{B]\i{x,y)x^) for < j < N. This replaces the top row in the determi- 
nant by (0 7j 7j-|-i ... 7jv+j). The determinant is unchanged if we subtract the row 
{y^ 7j ... 7Ar+j) from this row. That is, E^{BN{x,y)x^) is given by the determinant 
with top row {—y^ ... 0). Thus, 

E^{Bn{x, y)x^) = -h~l^{~y^)hn+i = y^ ■ 

So, Bi^{x,y) is a reproducing kernel. For any polynomial P{x) of degree A^ or smaller, 
E,{BN{x,y)P{x)) = P{y). But 

N 

Bn{x, y) = J2 EABNix, y)Pj{x))Pj{x) = Bn{x, y) 
j=o 

since {Pj{x)}^^q is an orthonormal basis in the polynomials of degree A^ or less. Thus 
(A.22) is proven. D 
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In terms of the 7*^^-' moment problems (with moments 7- = -^^^^), we recognize v 



as (72)""^/in-i(7^^^)- By (A.17) for the 7^^) problem, 



It 



n-l 



ET^(o) 



Vr 



-1(7^^)) 



But /in(7^'^) = ^^n(7)(72)"" and t;n-i(7^'^) = ll'^Vn-ih) where 

/ 74 ••• 7n+2\ 

75 • • • 7n+3 

l/n-1 = 

\ln+l ■■■ 72n 

Thus, using Proposition 5.13: 
Theorem A. 7. 



(A.23) 



■'^iE/Ao?){E/FW)--^-;^, 



In particular, the Hamburger problem is determinate if and only z/lim^^oo jr~^ = 00. 
So we have a simple ratio of determinants to determine determinacy. 



Appendix B: The Set of Solutions of the Moment Problem as a Compact 
Convex Set 

In this appendix, we will prove (following [1]) that A4^ {'^) and M.^{'^) are compact 
convex sets whose extreme points are dense. Each set is a subset of A^_|_(M U {00}), the 
set of measures on the compact set MUJoo}. This set, with the condition f dp{x) = 70, is 
a compact space in the topology of weak convergence (i.e., / f{x) dpn{x) -^ f f{x) dp{x) 
for continuous functions / on M U {c)o}). 

Theorem B.l. Ai^ {•y) and M.^{'y) are closed in the weak topology and so are compact 
convex sets. 

Remark. Since the x"^'s are unbounded, this does not follow from the definition of the 
topology without some additional argument. 

Proof. By Propositions 4.4 and 4.13, p G Ai^ {•y) if and only if 

f {x - zo)-^ dp{x) e D{zo) (B.l) 

for all zq G C. The set of /u's that obey (B.l) for a fixed zq is closed since D{zo) is 
closed and {x — zq)~^ is in C{M. U {00}). Thus, the intersection over all zq is closed. 
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•^"^(7) = M^i'j) n {/U I j f{x)dii{x) = if supp/ C (—00,0)} is an intersection of 
closed sets. D 

Remark. To get compact sets of measures, we need to consider M U {00} rather than M. 
Thus a priori, the integral in (B.l) could have a point at infinity giving a constant term 
as zq = iy with y ^ 00. Since D{zo) -^ {0} as zq = iy with y ^ 00, this term is absent. 

Theorem B.2 (Naimark). fi E A4^{'y) (resp. ^^"^(7)) is an extreme point if and only 
if the polynomials are dense in L^CR, dp). 

Remark. Compare with density in L^(M, dp) which picks out the von Neumann solutions. 

Proof. This is a simple use of duality theory. The polynomials fail to be dense if and 
only if there exists a non-zero F G L°° (M, dp) so that 

f x''F{x)dp{x)=0 (B.2) 

for all n. We can suppose that ||-F||oo = 1, in which case dp± = (1 ± F) dp both lie in 
■M.^ il) with p = ^{p+ + p_), so p is not extreme. 

Conversely, suppose p = ^{p+ + p-) with p+ 7^ p_ and both in ^^''^(7). Then 
P± ^ 2/9 so /9_|_ is p-absolutely continuous, and the Radon-Nikodyn derivative, -^, obeys 
ll^lloo < 2. Let F = 1 - ^, so ||F||oo < 1 and F ^ since p+ ^ p_. Then 



F(x)x" dp= x'^dp- x"" dp+ = 0. 

Thus, we have proven the result for Ai^{'~f). Since TW (7) = {p E A4^ {"f) \ 
J f{x)dp{x) = for / in C(R U {00}) with support in (— oo,0)}, M.^{'y) is a face of 
Ai^ {'-f), so extreme points of J^^{'y) are exactly those extreme points of A4^ {'j) that 
lie in A^'^(7). D 

Theorem B.3. Let p E Ai^ {•y) have ord(p) < 00. Then p is an extreme point. 

Proof. By Theorem 6.4, for some zi,...,Zn E C_(_, we have that the polynomials are 
dense in L^(R, H^+i 1^ ~ ^j\~'^dp). For any polynomially bounded continuous function 



j \f{x)\dp{x)< l^j \f{x)\^l[\x-z,\-^dp{x)j l^j l[{x-z,)^dp{x)j 



by the Schwarz inequality. It follows that the identification map is continuous from 
L^(M, n?=i(^ ~ '^)~^ dp) i'^to -L^(]R, dp), so the polynomials are dense in L^-norm in the 
continuous functions, and so in L^{M.,dp). By Theorem B.2, p is an extreme point. D 
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Theorem B.4. The extreme points are dense in A4^ {"f) {and in A4^{'~f)). 

Proof. The finite point measures are dense in the finite measures on R U {oo}. Thus, by 
the Herglotz representation theorem in form (6.8), if $ is a Herglotz function, there exist 
real rational Herglotz functions $„ so that ^n{z) -^ ^(-2) for each z G C-|-. 

Now let p G Ai^ {'^) and let $p be the Nevanlinna function of p. Let ^n be as above 
and pn = p$„. Then by (4.35), Gp^{z) — > Gp{z) for each z G C_|_ and so p^ ^> p weakly. 
By Theorems 6.4 and B.3, each such p^ is an extreme point. 

For the Stieltjes case, we need only note that by Theorem 4.18 and the remark after 
it, if p G 7Vf'^(7), then the approximating ^nS can be chosen so that /9$„ G 7^*^(7). D 

Theorem B.5. For any indeterminate set of Hamburger moments {7n}5^o^ M.^ {'^) 
has extreme points p with ord{p) = 00. 

Proof. We first pick positive aj strictly decreasing so that for any p G M.^ {'^) we have 
that 

/n 
Y[{l + a^jx^fdp{x) <2. (B.3) 

We can certainly do this as follows: Since the integral only depends on the moments 7^, 
we need only do it for some fixed po ^ ^^^(7). Since 

lim / (1 + a^x^)^ dp{x) = 1, 

we can pick cti > so 

{1 + ajx^f dp{x) < 2. 

We then pick a2,a3, ■ ■ ■ inductively so the integral is strictly less than 2. Then the sup 
in (B.3) is bounded by 2. The product must be finite for a.e. x w.r.t. dp, so 

00 
J^af<oo. (B.4) 



So if 

then for any real x, 
exists, and by (B.3), 

for any p G M.^ {'^). 



9n{x) = JJ(l + a 
i=i 



2 2\ 
^X ) 



g{x) = lim gn{x) 
g{x)'^ dp{x) < 00 
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Now pick Ci, C2, . . . in C_|_ and pi, p2, ■ • ■ in ^^-^(7) inductively so pk obeys 

{x-ia~^)~^dpk{x) =Cj, j = l,...,k (B.5) 

and then Ck+i is the middle of the disk of allowed values for J{x — ia~^i)~^ dp{x) for 
those p in ^^^(7) which obey (B.5). Pick a subsequence of the pj's which converges to 
some Poo G A4^ {'^). Then poo obeys (B.5) for all k. 
Let 

Mk = {pe M"{-f) I p obeys (B.5)} 



and 



Define for all n, m, 



Moo = nMk. 



/IL 
x^\[{l + ay)-Up^{x). 

By Theorem 6.1, p G 7Vl^(r*^'^^) if and only if dp = Qn dp lies irv M.k- Tm is decreasing 

r^^= fx^g{x)-Up{x) 



to Tm and 



by a simple use of the monotone convergence theorem. 

We claim that p G Ai^ (T^^^) if and only if dp = g dp lies in A^oo- For 



7n = lim x"^g{x) ^gn{x)dp{x) 



and the right side only depends on the moments T^ and similarly for calculation of 
f{x — iaj)~^ dp{x). 

Now let p be a von Neumann solution of the T^°^' moment problem and dp = g dp. 
Since J g^ dp < 00 by construction, the proof of Theorem B.3 shows that p is an extreme 
point. On the other hand, since p obeys (B.5) and Q is not in the boundary of its allowed 
circle, ord(p) = 00. D 



Appendix C: Summary of Notation and Constructions 

Since there are so many objects and constructions associated to the moment problem, 
I am providing the reader with this summary of them and where they are discussed in 
this paper. 
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CI Structure of the Set of Moments. 

{7n}^=Q is called a set of Hamburger moments if there is a positive measure p on 
(—00, 00) with 'jn = f x'^ dp{x) and a set of Stieltjes moments if there is a p with support 
in [0,00). We take 70 = 1 and demand supp(p) is not a finite set. Theorem 1 gives 
necessary and sufficient conditions for existence. If there is a unique p, the problem is 
called determinate. If there are multiple p's, the problem is called indeterminate. If a 
set of Stieltjes moments is Hamburger determinate, it is a fortiori Stieltjes determinate. 
But the converse can be false (see the end of Section 3). 

M.^ {'^) (resp. Ai^{'y)) denotes the set of all solutions of the Hamburger problem 
(resp. all p G A4^ {'^) supported on [0, 00)). They are compact convex sets (Theorem B.l) 
whose extreme points are dense (Theorem B.4). Indeterminate Hamburger problems 
have a distinguished class of solutions we have called von Neumann solutions (called A^- 
extremal in [1] and extremal in [35]). They are characterized as those p G A4^ ("f) with 
the polynomials dense in L^(R, dp). They are associated to self-adjoint extensions Bt of 
the Jacobi matrix associated to 7 (see C4 below for the definition of this Jacobi matrix). 
The parameter t lies in R U {00} and is related to Bt and the solution pt £ M.^ {•y) by 

t={do,Bi^do) = fx-^dptix). (C.l) 



In the indeterminate Stieltjes case, there are two distinguished solutions: pp, the 
Friedrichs solution, and px, the Krein solution. Both are von Neumann solutions of 
the associated Hamburger problem and are characterized by inf[supp(/Ui?)] > inf[suppp] 
for any other p G A^^(7) and by /Uk({0}) > and /Uk({0}) > /o({0}) for any other 
p G M.^ {""y) (see Theorems 3.2 and 4.11 and 4.17). All solutions of the Stieltjes problem 
lie between pp and px in the sense of (4.42) (see Theorem 5.10). 

We define 

C+ = {^ G C I Im^ > 0}. 

For any probability measure p, we define its Stieltjes transform as a function of 

C\supp(p) by 



X — z 
These functions map C-|- to C-(- and have an asymptotic series 



CXD 



G,{z)r^-z-'\^^,z-n (C.2) 

as \z\ -^ 00 with min(|Arg(2)|, |Arg(— 2)|) > e. (C.2) characterizes A^^(7) and the 
asymptotics are uniform over M.^ {^) (Proposition 4.13). 

For an indeterminate moment problem, we defined ord(p), the order of p, to be the 
codimension of the closure of the polynomials in L^(]R, dp). The solutions of finite order 
are dense in ^^^(7) and in ^^"^(7), and each is an extreme point. Solutions of finite 
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order (and, in particular, the von Neumann solutions) are discrete pure point measures; 
equivalently, their Stieltjes transforms are meromorphic functions (Theorems 4.11 and 

6.4). 

C2 Nevanlinna Parametrization. 

We used T to denote the analytic maps of C_|_ to C_|_ where the closure is in the 
Riemann sphere. The open mapping theorem implies that if $ G JF, either $(2) = t for 
some t G MU {00} or else $ is a Herglotz function which has a representation of the form 
(1.19). 

Associated to any indeterminate moment problem are four entire functions, A(^z\ 
B(z), C(z), D(z), obeying growth condition of the form \f{z)\ < Cg exp(e\z\). We defined 
them via the transfer matrix relation (4.16), but they also have explicit formulae in 
terms of the orthogonal polynomials, P and Q (defined in C4 below) — these are given 
by Theorem 4.9. 

We define a fractional linear transformation F{z) : C U {00} ^ C U {00} by 

^/ N/ N C(z)w + A(z) ,^ , 

F(z)(w) = --r-i Z-V4- C.3 

For indeterminate Hamburger problems, there is a one-one correspondence between 
p G M.^{'^) and $ G JF given by 

G (z) = F(z)mz)) = _.^M^M±^M (c 4) 

Wl^J ^[^)['^[z)) Diz)^iz)+B{z) ^ ' 

for z G C-i- (Theorem 4.14). $ is called the Nevanlinna function of p, denoted by $p. 
The von Neumann solutions correspond precisely to ^{z) = t where t is given by (C.l) 
(Theorem 4.10). The solutions of finite order correspond precisely to the case where $ 
is a ratio of real polynomials (Theorem 6.4). 

Given a set of Stieltjes moments, the p G ^^'^(7) are precisely those p G A4^ {'j) whose 
Nevanlinna function has the form (Theorem 4.18) 



r 

Jo 



dpj{x) 

X — z 



with J dp{x) < cxD and do>tp = j dppix) = {6Q,Ap^6o). 

For z G C_|_, the image of C_|_ U {c)o} under F{z) is a closed disk denoted by ^{z). 
Von Neumann solutions p have Gp{z) G dV{z) for all z while other solutions have Gp{z) G 
'D(z)™* for all z (Theorem 4.3, Proposition 4.4, and Theorem 4.14). 
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C3 Derived Moment Problems. 

Any set of moments has many families of associated moments. For each real c, 7(c) 
is defined by 

^r.{c) = J2(%'^r.-j. (C.5) 

There is a simple map (p ^-^ p( ■ — c)) that sets up a bijection between Ai^ {'^) and 
A4^ {'y(c)) and, in particular, 7 is Hamburger determinate if and only if 7(c) is Ham- 
burger determinate. But the analog is not true for the Stieljes problem. Indeed, (end of 
Section 3), if 7 is a set of indeterminate Stieltjes moments and cp = — inf supp((i/Ui?), 
then 7(c) is a set of Stieltjes moments if and only if c > cp, and it is Stieltjes determi- 
nate if and only if c = c^?. The orthogonal polynomials for 7(c) and 7 are related via 
P^/'^\z) = Pm{z-c). 

Given a set of Stieltjes moments {7n}5^0' o^^ defines Hamburger moments {F^jJ^o 
by 

F2m = 7m, F2m+1 = 0. (C.6) 

There is a simple map {dp i-^ diJ,{x) = |x[o,oo)(2;) dp{x^) -\- \x{-oo,o]{x) dp{x^)) that sets 
up a bijection between M.^ {^) and those p G M.^ {V) with dp{—x) = dp{x). Then 7 
is Stieltjes determinate if and only if the F problem is Hamburger determinate (The- 
orem 2.13). There are simple relations between the orthogonal polynomials for 7 and 
those for F ((5.89)-(5.92)). 
For £ = 1,2, ..., define 

7f = ^ . (C.7) 

If 7 is a set of Stieltjes moments, so are "f^^' for all £, and if 7 is a set of Hamburger 
moments, so are 7*^^^ for £ = 2, 4, ... . The map p 1-^ dp^^^ = x^ ^^^' is a map from A^^ (7) 

to M.^ {'^'^'^') which is injective. But in the indeterminate case, it is not^ in general, 
surjective. If 7 is indeterminate, so is 7*^^-', but the converse may be false for either the 
Hamburger or Stieltjes problems (Corollary 4.21). For general £, the connection between 
the orthogonal polynomials for 7 and 7*^^-' is complicated (see the proof of Proposition 5.13 
for £ = 2), but (see the remarks following Theorem 5.14) 

p(l). . _ dN[PN{Q)PN+l{z) - Pn+MPn{z)] 

Given a set of Hamburger moments {7j}°^0' ^"^^ defines a new set {7- }fLo by the 
following formal power series relation: 

00 -1 —1 ^00 \ 

Y.i-^TlnZ'' =1-11Z + (72 - ll)z' r£^-lT^^^^Z- . (C.8) 



■n=0 



n=Q 



p(0) 


{z) = 


■ aoQN+i{z). 


^(^) = 


. (^(o; 


i)(-^) 
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This complicated formula is simple at the level of Jacobi matrices. If {a^j^gi {^njj^o 

(resp. {an }5^0' {^" }5^o) ^^^ ^^^ Jacobi matrix coefficients associated to 7 (resp. 7*-°''), 
then 

a^°) = a^+i, 6^°^ = ^n+i- 

The map p ^^ p given by (5.53) sets up a bijection between M.^ {'~^) and A^''^(7'^'^^), and, 
in particular, 7 is a Hamburger determinate if and only if 7*^*^^ is (Proposition 5.15). This 
is not, in general, true in the Stieltjes case (see Proposition 5.17). There is a simple 
relation between the orthogonal polynomials for 7 and for 7*^^^^; see Proposition 5.16. In 
particular, 

(C.9) 
For I = —1, —2, . . . , we defined 

(CIO) 



C4 Associated Polynomials. 

Given a set of Hamburger moments, the fundamental orthogonal polynomials Pn{z) 
are defined by requiring for any p G M.^ (cy) that 

P,{x)Pt{x)dp{x)=5,t (C.ll) 

and Pn{x) = cnnx^ + lower order terms with cnn > 0. They obey a three-term 
recursion relation 

xPn{x) = aNPN+i{x) + bNPNix) + aN-iPN-i{x). (C.12) 

The Jacobi matrix A associated to 7 is the tridiagonal matrix given by (1.16). 

The associated polynomials Qn{x) are of degree A^ — 1 and defined to obey (C.12) 
but with the starting conditions, Qo{x) = 0, Qi{x) = -^. They are related to P/v by 
(Theorem 4.2) 

J x-y 

for any p G ^^^(7). As noted in (C.9), they are up to constants, orthogonal polynomials 
for another moment problem, namely, 7'-'^''. 

In the Stieltjes case, we defined polynomials Mn by 

M,(.) = P,(.) - ^^^P^^ . (C,14) 

They vanish at x = so — ^^^-^ is also a polynomial. It is of degree A^ — 1 and is up 
to a constant, the principal orthogonal polynomial of the moment problem 7*^^-*. N]\[{x) 
is defined analogously to (C.13) with Pn replaced by Mn. In Section 5, it was useful 
to change the normalizations and define multiples of Pn{—x), Qn{—x), Mjv(— x), and 
Nn{—x) as functions Un{x), Vn{x), Gn{x), Hn{x) given by (5.71)-(5.74). 
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C5 Pade Approximants and Finite Matrix Approximations. 

If {7n}5^o i^ ^ ^^^ of Hamburger (resp. Stieltjes) moments, the formal power series 

Yl'^=oi^^)^lriZ^ is called a series of Hamburger (resp. a series of Stieltjes). Formally, it 

sums to 

dp{x) 

1 + xz 

if p e A4^{'j). The Fade approximants f^^'^\z) to these series, if they exist, are defined 
by (5.25)-(5.27) as the rational function which is a ratio of a polynomial of degree A^ to 
a polynomial of degree M , whose first N + M + 1 Taylor coefficients are {(— l)"7n}n=o • 
For a series of Stieltjes, for each fixed £ = 0, ±1, ... , 



N 



lim /[^+^-^'^l(-)-/.W 



exists for all z G C\(— oo,0) and defines a function analytic there. Indeed, for x G 
[0, oo), {—ly f'-^~^^~^'^\x) is monotone increasing (Theorem 6, Theorem 5.14, and The- 
orem 5.18). Moreover, 



f dprix) f dpK{x 

foiz) = / -r— , fi[z) = / -— — 

J 1 + xz J 1 + xz 



xz 

where pp, Pk are the Friedrichs and Krein solutions of the moment problem. In partic- 
ular, the moment problem is Stieltjes determinate if and only if /o = /i- 

There is a connection between y[-^+^-i.-^] (2;) and the moment problems "f^^' and '~f^^~^> 
(if £ < 0, "f^^' and 'j^^^^') given in Theorem 5.14 and Theorem 5.18. In particular, if 
i > and 7*^^^ is Stieltjes determinate, then f-i{z) = fo{z) = ■ ■ ■ = f£j^i{z), and if £ < 
and 7*^^'' is Stieltjes determinate, then fi{z) = fo{z) = ■ ■ ■ = fe{z) = fe-i{z). 

The situation for series of Hamburger is more complicated, but, in general, 
j![N+£-i,N]^^^ converges for z G C+ and i — ±1, ±3, . . . (Theorem 5.31 and Theo- 
rem 5.32; see Theorem 5.31 for issues of existence of the Fade approximant). 

There are connections between the Fade approximants and the polynomials P, Q, M, A^ 
as well as two finite matrix approximations Ap ' and A'^ ' to the Jacobi matrix, A. Ap ' 
is the upper right N x N piece of A. A'^ ' differs by adjusting the A^A^ matrix element 
of A^p^ so det(A^^) = 0. Then 






lN,N](^^_/x^ n , ^a[^+^]\-1x.\- z^Nn+i{-^) 



The connection to the continued fractions of Stieltjes is discussed after Theorem 5.28. 
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C6 Criteria for Determinacy. 

Criteria for wtien a Hamburger problem is determinate can be found in Proposition 1.5, 
Theorem 3, Corollary 4.5, Proposition 5.13, Theorem A. 6, and Theorem A. 7. 

Criteria for when a Stieltjes problem is determinate can be found in Proposition 1.5, 
Theorem 7, Corollary 4.5, Theorem 5.21, Corollary 5.24, Theorem A. 4, and Theorem A. 6. 
In particular, if one defines 

C2,.=£^. = _[a,_iP,(0)P,_i(0)]-i 
C2J-1 =mj = |P,_i(0)p 

(the c's are coefficients of Stieltjes continued fractions; the m's and ts are natural pa- 
rameters in Krein's theory), then the Stieltjes problem is determinate if and only if 
(Theorem 5.21 and Proposition 5.23) YlTLo ^j ~ ^^• 
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